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A  HYPERPSECTRAL  GAS  ANALYSIS  SYSTEM  (HyGAS) 

FINAL  REPORT 


1.0  Summary 

The  Hyperspectral  Gas  Analysis  System  (HyGAS)  Phase  II  project  objective  was  to 
design  and  develop  a  software  system  for  use  with  the  Hyperspectral  Image  Processing 
System  (HIPS)'  to  do  standoff  (remote)  gas  detection,  identification  and  analysis. 

To  define  practical  specifications  for  HyGAS,  we  reviewed  existing  techniques  for 
standoff  gas  analysis.  Techniques  developed  by  the  U.  S.  Army  for  application  in  field 
scenarios  (as  opposed  to  laboratory  settings)  were  given  precedence  for  incorporation  into 
HyGAS. 

The  following  techniques  from  U.  S.  Army  researchers  were  initially  identified  as 
operationally  functional  for  limited  standoff  gas  applications  under  actual  field  situations. 

•  gaussian  spatial  edge-detection,  background  clutter  suppression,  and  multiple 
focal  plane  filters  (M.  Althouse,  ERDEC1 2) 

•  linear  discriminant  spectral  detection  (R.  Kroutil,  ERDEC) 

•  narrow-band  infrared  video  camera  (M.  Norton,  NVEOD3) 

•  neural  net  spectral  identification  (P.  Rauss,  ARL4) 

The  two  spectral  techniques  (Kroutil  and  Rauss)  were  demonstrated  using  point  FOV 
FTIR  spectrometers.  The  remaining  techniques  (Althouse  and  Norton)  were 
demonstrated  using  infrared  spatial  imaging  arrays  (FLIRs  or  IR  video  cameras),  in 
conjunction  with  one  or  more  associated  narrow-band  filters.  Although  each  of  these 
techniques  holds  in  common  the  potential  for  real-time  application,  each  is  limited  by 
currently  available  hardware  to  producing  data  in  only  two  of  the  three-dimensions 
accommodated  by  a  standard  HIPS  image  cube5. 

We  conceptualized  HyGAS  as  an  interactive  software  package  for  HIPS  that  would 
enable  the  user  to  interactively  display  and  analyze  both  multiple  images  and  multiple 
spectra.  Spectral  techniques  could  be  applied  interactively  to  selected  spatial  images; 
similarly,  spatial  techniques  could  be  applied  interactively  to  selected  spectral  bands. 

After  investigating  each  of  the  above  techniques,  only  the  two  spectral  techniques  proved 
feasible  for  incorporation  into  the  existing  HIPS  architecture  of  performing  spectral 
operations  on  a  traditional  2-D  spatial  scene.  To  simulate  the  use  of  these  techniques 
with  next-generation  hardware  systems  capable  of  generating  full  3-D  hyperspectral 
image  cubes,  we  created  simulated  hyperspectral  image  cubes  based  on  the  actual 
Government  Furnished  Information  (GFI)  field-acquired  spectral  data  for  each  technique. 


1  Hyperspectral  Image  Processing  System,  developed  under  previous  SBIR  funding  for  the  U.S.  Army 
Topographic  Engineering  Center,  Alexandria,  Virginia. 

2  U.  S.  Army  Edgewood  Research  Development  and  Engineering  Center,  Edgewood,  Maryland. 

3  U.  S.  Army  Night  Vision  and  Electro-Optical  Directorate,  Ft.  Belvoir,  Virginia. 

4  U.S.  Army  Research  Laboratory,  Ft.  Belvoir,  Virginia 

5  a  HIPS  image  cube  is  usually  represented  as  a  two-dimensional  spatial  image  with  full  spectral  data 
defined  for  each  image  pixel. 
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We  then  designed  the  interactive  HyGAS  user  interface  to  incorporate  each  function 
necessary  to  apply  these  techniques.  Details  on  supplied  GFI  datasets,  gas  databases,  and 
creation  of  simulated  cubes  are  given  in  the  appendices  of  this  report.  Details  on  the 
functional  use  of  HyGAS  are  given  separately  in  the  HyGAS  User's  Manual/T utorial. 

2.0  Technical  Objectives 

HyGAS  is  a  software  module  prototype  designed  to  meet  the  scientific  requirements  of 
standoff  (remote)  gas  detection  and  exploitation  for  an  imaging  spectrometer  system. 
HyGAS  applies  the  methodologies  used  to  process/analyze  data  from  point  spectrometers 
to  process  a  hyperspectral  image  cube  data  set,  i.e.,  a  2-D  image  with  full  spectral  data  for 
each  image  pixel. 

HyGAS  is  designed  to  function  within  the  framework  of  the  existing  HIPS  software 
package. 

The  HyGAS  Phase  II  effort  was  divided  chronologically  into  two  major  parts: 

(1)  define  the  practical  specifications  required  for  HyGAS  development,  as  built 
upon  the  Phase  I  background  research  and  system  requirements 

(2)  implement  these  specifications  into  the  HyGAS  prototype  software  system. 

HyGAS  development  comprised  four  tasks: 

Task  1 — System  Development 
Task  2 — Software  Development 

Task  3 — Spectrum  Database/Measurements  Development 
Task  4 — Documentation 

2.1  Task  1 — System  Development 

The  objectives  of  the  system  development  task  were: 

(1)  prepare  a  HyGAS  conceptual  design 

(2)  define  the  specifications  necessary  to  implement  this  design 

(3)  incorporate  these  specifications  into  preliminary  and  final  designs  for  the  HyGAS 
system. 

This  task  included  reviewing  the  research  conducted  during  Phase  I  and  conducting  new 
research,  selecting  suitable  applications  and  requirements  from  that  research,  and  defining 
the  specifications  required  to  implement  software  for  these  applications  and  requirements. 

The  HyGAS  software  system  was  developed  for  hardware  systems  that  did  not  currently 
exist,  but  were  expected  to  become  available  in  the  near  future.  The  design  of  HyGAS 
was  biased  towards  hardware  systems  that  were  currently  under  development,  and 
systems  that  were  likely  to  adapt  well  to  existing  methods  of  gas  study.  Such  systems 
include  imaging  fourier  transform  infrared  (FTIR)  systems  (FTS),  and  hyperspectral 
imaging  (HSI)  cameras.  To  adapt  best  to  gas  studies,  such  systems  should  operate  in  the 
8-12  pm  spectral  region  with  a  spectral  resolution  of  at  least  8  cm'1.  Systems  such  as 
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these  will  allow  the  application  of  novel  techniques  of  gas  detection  and  integration 
currently  in  use  for  point  FTIR  spectrometers  and  multispectral  imaging  (MSI)  cameras. 

HyGAS  was  designed  to  incorporate  into  HIPS  current  methods  of  standoff  gas  analysis 
using  both  spectral  and  spatial  (image)  processing  techniques.  Primary  emphasis  was 
gas/vapor  cloud  detection,  i.e.,  “Is  there  a  gas  cloud  present,  and  if  so,  where  is  it?” 
Secondary  emphasis  was  gas  identification,  i.e.,  “If  there  is  a  gas  contaminant  present, 
what  gas(es)  is  it?”  Tertiary  emphasis  was  gas  analysis,  i.e.,  “How  much  gas  is  there?”, 
“Where,  in  detail,  is  the  gas?”  and/or  “Where  is  the  gas  cloud  likely  to  go?”  Gas  cloud 
detection  was  currently  being  done  by  both  spectral  (Kroutil)  and  spatial  (Althouse) 
processing  techniques.  Gas  identification  was  only  being  done  by  spectral  (Rauss) 
techniques. 

Gas  analysis  usually  requires  both  spectral  and  spatial  techniques  as  well  as  some  sort  of 
depth  profiling  technique  (such  as  LIDAR),  or  some  sort  of  time  profiling  technique 
(such  as  real  time  video  synchronization,  e.g.  Norton).  As  delivered,  only  spectral 
techniques  were  implemented  into  HyGAS.  Spatial  techniques  and  depth  and  time 
profiling  techniques  should  be  considered  for  future  implementation. 

2.2  Task  2 — Software  Development 

The  objective  of  the  software  development  task  was  to  implement  the  software  processes 
required  to  build  the  HyGAS  system  defined  in  Task  1.  This  task  included  acquisition  of 
necessary  hardware,  acquisition  and  modification  of  necessary  existing  software,  and 
coding  of  additional  new  software. 

A  Gas  Analysis  option  was  added  to  the  HIPS  main  menu  selection  window.  This  Gas 
Analysis  menu  contains  the  individual  gas  analysis  functions,  which  can  be  run  non- 
interactively.  It  also  includes  the  Interactive  HyGAS  module  which  incorporates  all  of 
these  individual  functions  (and  more)  into  the  interactive  gas  analysis  package. 

Individual  HyGAS  functions  include  Linear  Discriminant  Test  (Kroutil  gas  detection 
algorithm),  Background  Removal  (subtraction  or  division  method),  Interferogram  to 
Spectrum  and  Spectrum  to  Interferogram  (FFT  and  inverse  FFT),  Library  Search  and 
Neural  Net  search  (gas  identification/analysis),  and  image  statistics.  Most  of  these  are 
spectral  functions  which  can  be  applied  either  to  an  image  cube  or  an  individual  spectrum 
(or  interferogram).  [Additional  functions  originally  considered  for  incorporation,  such  as 
the  spatial  edge  detection  functions  from  M.  Althouse,  turned  out  not  to  be  feasible  for 
this  implementation  of  HyGAS.] 

2.3  Task  3 — Spectral  Database/Measurements  Development 

The  objectives  of  the  spectral  database/measurements  development  task  were  to  provide 
suitable  data  for  testing  and  demonstrating  the  HyGAS  system  developed  under  Task  1, 
and  for  demonstrating  and  validating  the  performance  of  each  software  process  developed 
under  Task  2.  This  task  included  acquisition  of  existing  datasets  and/or  data  libraries, 
measurement  of  new  datasets  and/or  data  libraries,  and/or  modeling  of  theoretical  datasets 
and/or  data  libraries.  Datasets  and  gas  libraries  supplied  (GFI)  for  use  with  HyGAS  are 
described  in  detail  in  Appendices  A  and  B. 
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Hyperspectral  image  data  suitable  for  gas  applications  was  not  available  for  HyGAS. 
Therefore,  we  synthesized  datasets  to  demonstrate  and  test  the  functionality  of  all 
HyGAS  routines.  These  datasets  are  based  on  actual  spectral  data  synthesized  into  a 
hyperspectral  image  cube  format  and  are  designed  to  approximate  realistic  acquisition 
applications.  Details  are  given  in  Appendix  C. 

2.4  Task  4 — Documentation 

The  objective  of  the  documentation  task  was  to  manage,  plan,  track  and  provide  proper 
closure  of  the  project.  This  task  included  preparation  and  delivery  of  the  items  on  the 
contract  data  requirements  list  (CDRL).  This  task  also  included  the  test  and  validation 
plans,  training  and  demonstration,  the  Factory  Acceptance  Test  (FAT)  and  the  Site 
Acceptance  Test  (SAT). 

Documentation  included  technical  progress  reports,  work  plans,  software  users  manual 
and  final  technical  report. 

3.0  Review  of  Work  Performed 

3.1  Task  1 — System  Development 

HyGAS  was  designed  to  process  hyperspectral  image  data  for  standoff  (remote) 
gas/vapor  cloud  detection,  identification  and  analysis.  During  HyGAS  development,  no 
suitable  hyperspectral  imaging  sensor  systems  were  operational  for  these  applications. 
Hyperspectral  imaging  systems  did  exist,  but  they  had  neither  the  spectral  coverage  nor 
the  spectral  resolution  necessary  for  gas-specific  applications.  However,  gas  sensor 
systems  with  the  required  spectral  coverage  and  resolution  are  in  design  and  development 
stages,  and  the  HyGAS  project  intent  was  to  provide  software  for  processing  the  data 
from  such  systems  as  soon  as  they  become  available. 

Software  routines  developed  for  HyGAS  were  based  on  current  state-of-the-art 
techniques  developed  for  gas  analysis  using  existing  sensor  systems.  Existing  sensor 
systems  designed  for  standoff  (remote)  gas  applications  fall  into  one  of  the  following 
system  types:  Type  I — hyperspectral  point  spectrometers,  or  Type  II — multispectral 
imaging  cameras.  Under  development  are  hybrid  systems  of  Type  III — hyperspectral 
imaging  systems. 

All  of  these  systems  are  designed  to  operate  in  the  thermal  IR  spectral  region  (2-14/tm), 
often  referred  to  as  the  “fingerprint”  spectral  region  for  gas  studies.  Type  I  systems, 
hyperspectral  point  spectrometers,  are  usually  Fourier  Transform  Spectrometer  (FTS) 
systems.  These  systems  produce  sufficient  (high  resolution)  spectral  information  for  gas 
studies,  but  only  for  a  single,  usually  wide,  instantaneous  spatial  field  of  view  (IFOV). 
Type  II  systems,  multispectral  imaging  cameras,  are  usually  Forward  Looking  Infrared 
Radiometers  (FLIRs)  or  IR  video  cameras  converted  to  “multispectral”  spatial  imagers. 
These  systems  produce  sufficient  spatial  information,  but  only  at  a  very  few  (2-3) 
relatively  narrow  spectral  bands.  Type  III  systems,  hyperspectral  imaging  systems,  are 
currently  under  design  and  development  only  as  system  prototypes.  These  systems  will 
produce  both  sufficient  spectral  and  spatial  information  for  gas  specific  applications,  but 
at  the  expense  of  time  dependence  in  one  dimension.  Hyperspectral  imaging  hardware 
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systems,  by  the  inherent  restriction  of  a  2-D  focal  plane,  must  compromise  one  of  the 
three  dimensions  (spatial-spatial-spectral)  in  a  hyperspectral  image  cube,  and  collect  that 
dimension  as  a  function  of  time.  Spatial  analysis  techniques  are  most  useful  for  hardware 
that  generates  a  spatial-spatial  focal  plane  (i.e.  imaging  cameras).  Spectral  analysis 
techniques  are  most  useful  for  hardware  that  generates  a  spatial-spectral  focal  plane 
(imaging  spectrometers).  The  remaining  dimension  (spectral  for  imaging  cameras,  spatial 
for  imaging  spectrometers)  must  be  collected  (scanned)  as  a  function  of  time.  This 
naturally  creates  difficulties  with  time  variable  phenomena,  which  necessitates  that 
integration  times  and  time  differentials  be  kept  as  short  as  possible. 

In  the  HyGAS  project,  we  integrated  gas  detection  algorithms  and  techniques  already  in 
practice  with  Type  I  and  II  systems,  into  the  existing  HIPS  software  system  designed  to 
handle  data  from  Type  III  systems.  Thus  HyGAS  is  an  integrated  software  system  ready 
to  process  data  from  hyperspectral  spatial  imaging  sensor  systems  with  the  spectral 
resolution  and  range  of  coverage  required  for  standoff  gas  applications. 

During  Phase  II,  we  selected  the  techniques  in  use  for  Type  I  and  II  systems  as  most 
likely  to  be  successfully  implemented  into  HyGAS.  A  certain  bias  was  assumed  to 
support  systems  of  Type  III  most  likely  to  be  available  in  the  near  future. 

3.1.1.  Type  I  systems — Hyperspectral  point  spectrometers 

Spectral  analysis  of  gases  and  vapors  is  traditionally  done  in  the  thermal  infrared  spectral 
region  (700-5000  cm-',  or  2-14  /mi),  at  high  spectral  resolution  (0.1-16  cm-',  or  AMX  = 
0.0001-0.01),  using  integrated  point  field  of  view  laboratory  Fourier  Transform 
Spectrometer  (FTS)  systems.  The  infrared  spectral  region  is  often  referred  to  as  the 
“fingerprint”  region  for  gas  identification:  each  molecule  exhibits  specific  characteristic 
absorption  (or  emission)  lines  due  to  specific  changes  in  energy  associated  with  the 
electronic,  vibrational  or  rotational  state  of  the  molecule.  For  gases,  changes  in 
vibrational  and  rotational  states  are  most  common.  Changes  in  vibrational  state  exhibit 
broader  bands  than  changes  in  rotational  state  (rotational  fine  structure),  and  the  two  are 
often  superimposed,  creating  a  vibration-rotation  band.  Rotational  fine  structure  (lines  of 
bandwidth  -0.1  cm-')  is  important  primarily  for  the  detailed  study  of  small,  symmetrical 
molecules.  For  the  study  of  larger,  more  complicated  molecules,  so  many  broad, 
overlapping  lines  are  generated  due  to  vibrational  changes,  that  high  resolution  gains  very 
little  over  coarser  resolution.  For  studies  of  these  agents,  resolutions  on  the  order  of  16 
cm-1  can  suffice. 

Laboratory  observations  are  generally  conducted  under  carefully  controlled  atmospheres, 
with  the  major  absorbing  atmospheric  constituents,  water  and  carbon-dioxide,  removed. 
For  standoff  field  applications,  the  strength  of  the  ambient  atmospheric  water  vapor 
absorptions  typically  leave  only  the  “window”  spectral  regions  of  800-1200  cm-1  (~8— 12 
fim)  and  2000-3000  cm-1  (~3-5/im)  useful.  Additionally,  the  reduced  intensity  of  the 
background  source  radiation  and  the  resultant  reduced  temperature  differential  between 
that  source  and  the  gas  under  study  (ambient  thermal  background  instead  of  a  hot  thermal 
laboratory  source)  reduces  the  signal-to-noise  ratio  to  near  detectability  limits  for  ambient 
atmosphere  field  applications. 
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FTS  systems,  with  their  inherent  Fellgett  and  Jacquinot  advantages,  are  generally 
preferable  to  dispersive  (grating)  systems  for  low  signal-to-noise  scenarios.  The  Fellgett 
advantage  refers  to  the  inherent  benefit  achieved  by  an  FTS  system  which  can  collect  the 
energy  from  all  spectral  frequencies  simultaneously,  over  that  of  a  comparable  single 
detector  dispersive  system,  which  could  only  collect  the  energy  from  each  spectral  band 
sequentially  as  the  spectrum  is  scanned  over  the  detector.  This  advantage  is  achieved 
with  the  FTS  system  by  frequency-division  multiplexing,  providing  wide  wavelength 
coverage  as  well  as  greatly  improved  SNR  (perhaps  several  orders  of  magnitude)  over 
that  of  a  comparable  range  scanning  spectrometer.  The  Jacquinot  advantage  refers  to  the 
greater  optical  throughput  (etendue,  defined  as  the  product  of  the  beam  area  A  and  the 
solid  angle  Ci)  achieved  by  an  FTS  system  over  a  comparable  resolution  dispersive 
system.  This  advantage,  shared  also  by  Fabry-Perot  interferometers,  is  realized  mostly 
on  practical  considerations  of  the  physical  size  of  a  reasonable  beamsplitter  element  for 
an  FTS  system,  as  compared  to  the  much  larger  grating  required  for  a  similar  resolution 
dispersive  system. 

In  the  laboratory,  FTS  systems  produce  sufficiently  high  signal-to-noise  data  at  high 
enough  spectral  resolution  to  make  gas  identification  and  analysis  easy  and  practicable. 
Data  collected  in  time  (Fourier)  domain  are  transformed  into  spectral  (frequency  or 
wavenumber)  domain  using  a  Fast  Fourier  Transform  (FFT)  technique.  Data  of  both  the 
target  and  a  suitable  background  reference  are  required  in  order  to  remove  instrumental  or 
observational  effects  and  isolate  the  spectral  features  of  the  target  gas  itself.  This  is 
referred  to  as  background  removal,  and  is  typically  performed  on  the  target  measurement 
by  dividing  by,  or  in  some  cases  subtracting,  the  background  measurement.  By  choosing 
an  appropriately  calibrated  background  reference,  the  resultant  spectrum  of  the  target  can 
be  calibrated  to  relative  absorbance,  and  searched  against  a  standard  library  of  likely 
candidate  materials  to  be  identified  based  on  traditional  search  and  match  techniques. 

In  the  field,  where  the  signal-to-noise  ratio  is  much  lower  than  in  the  laboratory,  the 
traditional  library  search  and  match  techniques  fail.  Two  promising  methods  of  dealing 
with  this  difficulty  were  selected  for  implementation  into  HyGAS:  (1)  a  technique 
developed  by  R.T.  Kroutil  of  the  U.S.  Army  Edgewood  Research  and  Development 
Engineering  Center  (ERDEC),  which  is  applied  to  raw  FTS  interferograms  to  detect  the 
presence  or  absence  of  a  particular  gas;  and  (2)  a  technique  developed  by  Patrick  J.  Rauss 
of  the  U.S.  Army  Research  Laboratory  (ARL),  which  is  applied  to  the  transformed  FTS 
spectra  to  identify  an  unknown  gas. 

3. 1.1.1  Gas  Detection 

Kroutil  developed  a  unique  technique  suitable  for  gas  detection  in  actual  real-time  field 
scenarios,  using  raw  FTS  interferogram  data.  Based  on  a  predetermined  linear 
discriminant  vector,  a  simple  mathematical  test  is  performed  in  the  field  on  each 
interferogram  acquired,  detecting  the  presence  or  absence  of  one  particular  gas.  The  key 
to  this  technique  is  the  advance  determination  of  the  linear  discriminant  vector.  The 
accuracy  of  this  method  depends  on  the  quality  and  universality  of  this  vector.  In  theory , 
if  one  can  measure  a  sufficiently  diverse  set  of  background  materials  in  advance,  both 
with  and  without  gas  presence,  to  define  this  vector  with  sufficient  accuracy,  one  can 
attain  a  very  high  (-98%)  positive  detection  ratio  with  very  few  false  alarms.  It  is  then  a 
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simple  mathematical  task  to  use  this  discriminant  in  the  field  to  detect  the  presence  of  the 
desired  gas.  The  disadvantage  of  this  technique  is  that  a  different  discriminant  vector 
must  be  determined  for  each  particular  gas  one  wants  to  search  for. 

This  technique  is  very  good  at  dealing  with  the  difficulties  of  low  signal-to-noise  ratio 
and  changing  backgrounds  common  to  typical  field  data  collection  scenarios.  This 
technique  also  is  fast  because  it  operates  on  the  raw  FTS  interferogram  data,  obviating 
any  need  to  perform  an  FFT,  and  it  is  space  efficient  because,  in  general,  it  operates  on 
only  a  small  subset  of  the  interferogram  (i.e.  100  instead  of  1000  channels).  While  not  of 
particular  importance  to  a  point  FTS  system,  this  consideration  becomes  much  more 
important  for  an  imaging  FTS  system.  The  trade-off  to  obtain  this  speed  of  application  is 
twofold:  (1)  extensive  measurement  and  computation  is  required  in  advance  of  the  actual 
field  application  in  order  to  determine  an  accurate  discriminant  vector  for  any  agent;  and 
(2)  the  “search”  is  actually  the  detection  likelihood  of  one  single  specific  agent — 
“searching”  against  a  “library”  of  agents  would  be  very  fast  if  such  a  library  were  to  exist, 
but  developing  such  a  library  would  require  the  extensive  advance  measurement  and 
computation  effort  for  each  individual  agent  included  in  the  library.  Although  this 
technique  is  currently  impracticable  for  more  than  a  very  few  individual  agents,  future 
prospects  may  be  improved  by:  (1)  the  use  of  suitable  hyperspectral  imaging 
spectrometers  to  collect  a  sufficient  advance  measurement  database  for  each  desired 
agent;  and  (2)  the  use  of  a  super  computer  to  speed  up  the  advance  computation  time. 
Also,  in  practicality,  the  number  of  agents  one  would  be  searching  for  should  be 
relatively  small  for  most  applications,  and  development  of  a  suitable  “library”  of 
discriminant  vectors  would  be  for  only  a  few  agents. 

For  HyGAS,  we  designed  the  user  interface  for  HIPS  to  perform  the  discriminant  test 
using  Kroutil's  code  and  discriminants.  We  provide  this  “discriminant  test”  function  for 
application  to  a  raw  interferogram  “image”  dataset.  We  generated  a  suitable  input  test 
dataset  of  interferograms  from  data  we  acquired  from  Kroutil.  We  provided  for  a  small 
“library”  of  discriminants,  which  we  populated  with  two  discriminants  for  test  purposes. 
We  obtained  these  discriminants  from  Kroutil.  [Details  are  given  in  Appendix  E.] 
Standard  HIPS  routines  for  display  can  be  used  to  identify  and  locate  the  detected  agent 
in  the  image  based  on  the  results  of  the  discriminant  test. 

For  HyGAS,  it  was  not  practicable  to  implement  the  actual  advance  pattern  recognition 
techniques  developed  by  Kroutil  to  define  a  discriminant  vector.  This  situation  might  be 
feasible  if  one  acquired  a  Silicon  Graphics  machine  of  suitable  compatibility  to  Kroutil's. 

3. 1.1.2  Gas  Identification 

Rauss  developed  a  unique  technique  for  gas  identification  using  a  trained 
backpropagation  neural  network  to  search  a  processed  spectrum  from  an  FTS  system 
against  a  small  but  traditional  spectral  library.  This  technique  was  developed  by  training 
a  pattern  recognition  technique  to  recognize  noisy  processed  spectra  based  on  a  training 
dataset  of  library  spectra  corrupted  with  random  noise. 

For  HyGAS,  we  incorporated  this  technique  into  the  HyGAS  package  using  a 
commercially  available  neural  net  package.  We  initially  intended  to  acquire  the  trained 
net  and  identification  code  from  Rauss.  We  acquired  the  FFT  routine  as  part  of  the  initial 
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code  received  from  R.  Kroutil  to  process  any  interferogram  to  yield  a  “spectrum."  Data 
for  both  the  target  and  a  suitable  background  must  be  acquired,  transformed,  and  then 
processed  through  a  background  removal  technique.  Given  a  suitable  known 
background,  spectral  calibration  can  be  performed,  and  the  resultant  spectrum  can  be 
searched  using  these  special  techniques  “trained”  on  noisy  library  data  to  identify  the 
unknown  gas. 

Additional  details  on  this  technique  are  given  in  Appendix  D. 

3.1.2  Type  II  systems — Multispectral  imaging  cameras 

Image  analysis  of  gases  and  vapors  based  on  thermal  IR  spectral  signatures  is  a  relatively 
new  developmental  technique.  Traditionally,  image  analysis  uses  various  edge  detection 
algorithms  to  detect,  identify  and  analyze  discrete  targets:  buildings,  tanks,  missiles,  etc. 
These  detection  algorithms  can  sometimes  use  multispectral  information  to  enhance  their 
detection  accuracy.  Vapor  clouds  are  nebulous  entities,  not  discrete  targets,  and 
detection/identification  algorithms  that  work  for  discrete  targets  will  not  work  for  vapor 
clouds.  Also,  with  current  imaging  hardware,  identification  of  unknown  gases  is  not 
feasible  because  the  spectral  resolution  is  too  low  and  the  range  of  spectral  coverage  is 
too  small;  however,  detection  of  specific  gases  is  possible  in  cases  where  specific  filters 
can  be  selected  a  priori  to  enhance  the  spectral  contrast  between  a  specific  gas  and  most 
common  backgrounds. 

Current  hardware  systems  researched  for  potential  HyGAS  applications  were  designed  to 
be  incorporated  into  existing  forward  looking  infrared  radiometer  (FLIR)  video  systems, 
for  use  in  the  field  in  real-time  (i.e.  at  the  30Hz  frame  rate  of  standard  video  systems), 
and  with  the  intent  to  work  eventually  as  part  of  an  autonomous  unmanned  automatic 
sensor  system.  Hardware  approaches  were  designed  to  convert  existing  raster  scan  FLIR 
thermal  imaging  systems  into  a  sort  of  narrow  bandpass  spectrometer  (multispectral 
imaging  camera)  by  means  of  filtering  the  8-12  pm  input  radiation  before  it  reaches  the 
FLIR  detector.  Typically,  two  or  three  specific  filters  are  selected  to  “match”  the  spectral 
signatures  of  the  known  specific  gas(es)  to  be  studied.  One  filter,  the  “target”  filter,  is 
selected  centered  on  a  spectral  feature  contributed  purely  by  the  gas;  and  the  other,  the 
“background”  filter,  is  selected  just  to  one  side  of  this  spectral  feature.  The  intent  is  to 
image  the  “target”  gas  in  one  filter  and  the  “background”  in  the  other  filter(s),  such  that 
the  image  of  the  “background”  filter  can  be  subtracted  from  (or  divided  into)  the  image  of 
the  “target”  filter,  eliminating  all  background  clutter  except  for  the  gas  cloud. 

Three  different  hardware  systems  were  considered  for  potential  HyGAS  application. 

Two  of  these  are  pupil  plane  filtering  techniques;  the  other  is  a  focal  plane  filtering 
technique.  Of  the  pupil  plane  filtering  techniques,  one  uses  traditional  discrete 
interference  filters;  the  other  uses  a  tunable  Fabry-Perot  Etalon  filter. 

The  focal  plane  filtering  technique  (developed  by  M.  Althouse,  ERDEC)  deposits  thin 
film  coatings  in  discrete  physical  spectral  filters  directly  on  top  of  a  linear  detector  array. 
Resulting  video  images  are  “multispectral”  in  that  within  each  video  frame,  stripes  across 
the  image  represent  data  collected  at  different  spectral  bands.  In  theory,  a  cloud  with  the 
suspect  gas  contamination  would  appear  distinctly  “striped”  in  the  video  image,  while 
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everything  else  would  appear  nearly  normal.  Analysis  of  actual  data  of  this  type  is 
somewhat  difficult  for  humans,  but  almost  trivial  for  computers. 

The  discrete  filtered  pupil  plane  technique  (also  developed  by  M.  Althouse)  uses  one  or 
more  filters  in  a  stepped  or  spinning  filter  wheel  in  front  of  the  FLIR  lens.  The  spinning 
of  the  filter  wheel  is  synchronized  to  the  video  frame  rate.  Discrete  filters  are  selected 
based  on  the  gas(es)  one  expects  to  search  for.  Again,  image  analysis  by  computer  is 
easier  than  image  analysis  by  a  human  operator. 

The  tunable  Fabry-Perot  Etalon  filtered  pupil  plane  technique  (developed  by  M.  Norton 
of  NVEOD)  uses  two  circular  variable  filters  (CVFs),  one  broad  band  and  one  narrow 
band,  both  placed  in  front  of  the  FLIR  lens.  Again,  the  filter  selection  is  synchronized  to 
the  video  frame  rate.  However,  in  this  case,  the  filters  are  not  discrete  filters,  but  are 
tunable  to  any  position  across  the  8-12  /nn  spectral  range.  Theoretically,  this  means  that 
tuning  can  be  matched  quickly  and  easily  to  any  specific  gas;  this  gives  this  type  of 
system  much  more  flexibility  over  discrete  filtered  systems.  With  this  system,  the 
minimum  number  of  filters  required  to  detect  the  gas  are  selected.  This  saves  time  and 
processing.  The  system  is  currently  designed  as  an  aid  to  visual  detection  of  a  vapor 
cloud  in  real  time.  One  aid  to  the  operator  is  an  inverse  video  technique,  such  that  a 
warm  gas  detected  against  a  cool  sky  background  may  look  bright,  and  against  a  warmer 
terrain  background  may  look  dark.  With  inverse  video,  the  operator  can  make  the  gas 
look  bright  in  both  cases. 

The  practical  difficulties  of  all  of  these  systems  are  low  signal  to  noise,  and  high 
background  clutter.  Different  approaches  have  been  selected  to  deal  with  these 
difficulties.  Techniques  that  are  designed  to  function  with  the  standard  FLIR  video 
output  for  immediate  analysis  by  an  operator  must  operate  fast,  generally  at  the  standard 
30Hz  video  rate.  For  automatic  applications,  operations  can  conceivably  be  somewhat 
slower  and  significantly  computer  enhanced. 

One  simple  technique  for  image  enhancement  is  multiple  frame  averaging.  More 
complicated  techniques  generally  involve  statistical  methods  for  edge  detection  and 
background  clutter  suppression.  M.  Althouse  has  experimented  with  several  edge 
detection  and  background  clutter  suppression  techniques  suitable  for  application  to  vapor 
clouds  using  multispectral  images. 

3. 1.2.1  Target  detection  based  on  statistical  concentration  models 

Althouse  is  currently  developing  statistical  methods  suitable  for  real-time  autonomous 
target  detection  of  chemical  vapor  clouds,  using  data  from  focal  plane  spectrally  filtered 
FLIRs.  Althouse  has  demonstrated  that  statistical  Gaussian-distribution  vapor 
concentration  models  developed  by  R.E.  Warren  can  be  successfully  applied  to  low 
signal-to-noise  bi-spectral  (2-spectral  bands)  images  to  detect  vapor  clouds  (i.e.  targets 
with  neither  distinct  edges  nor  distinct  dimensions);  however,  these  methods  are  rather 
time  consuming  to  compute.  Althouse  has  derived  a  recursive  Kalman  filtering  technique 
based  on  the  Warren  model  that  computes  much  faster.  He  also  has  developed  adaptive 
filtering  preprocessing  algorithms  for  suppression  of  background  clutter  to  further 
enhance  the  usefulness  of  this  technique  in  low  signal-to-noise  situations.  Although  also 
computationally  time  consumptive,  the  application  of  systolic  array  algorithms  with 
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appropriate  parallel  processing  hardware  can  reduce  this  computation  time  to  the  point  of 
being  real-time  possible  in  a  deployable  FLIR  package. 

For  HyGAS  we  decided  that  these  techniques,  although  quite  promising,  were  sufficiently 
developmental,  and  until  they  are  further  refined  for  application  on  realistic  data,  it  was 
not  practicable  to  incorporate  them  into  the  HyGAS  system.  We  also  decided  that  the 
concept  of  a  focal  plane  multispectral  filter  was  not  compatible  with  the  initial  design 
concept  of  HyGAS.  Further  research  of  this  technique  would  be  required  before  it  would 
be  suitable  for  incorporation  into  HyGAS.  However,  this  approach  may  eventually  be 
applicable  to  pupil  plane  filtering  techniques. 

3. 1.2.2.  Other  techniques  developed  by  Althouse 

Althouse  also  experimented  with  edge  detection  using  multistage  predictive  coding 
(MPC),  and  also  with  image  thresholding  using  relative  entropy  theory.  These  techniques 
were  not  selected  for  integration  into  HyGAS. 

3. 1.2.3.  Target  detection  based  on  tunable  Fabry-Perot  filter  FLIRs 

Norton  (NVEOD)  experimented  with  tunable  Fabry-Perot  Etalon  filters  for  vapor  cloud 
detection  using  current  FLIR  systems.  This  approach  used  only  very  simple  image 
processing  techniques,  such  as  image  subtraction  (or  division)  to  do  background  removal, 
inverse  video,  and  image  coadding.  These  techniques  already  existed  within  the  HIPS 
system,  so  applying  them  to  HyGAS  is  simple.  We  developed  these  techniques  into 
Interactive  HyGAS  to  make  them  easier  to  use. 

3.1.3  Type  III  systems:  Hyperspectral  imaging  systems 

Thermal  infrared  hyperspectral  imaging  systems  suitable  for  gas  studies  were  not  yet 
available  during  HyGAS  development.  However,  such  systems  are  under  design  and 
development,  and  may  be  suitable  for  gas  studies  in  the  not  too  distant  future.  One 
potential  system  is  a  3-5/vm  FTS  system  under  modification  at  the  University  of  Hawaii. 
Another  potential  FTS  system  is  under  development  by  Bomem. 

For  HyGAS,  we  unsuccessfully  tried  to  acquire  actual  test  data  of  a  spatial  scene 
collected  in  interferogram  space  to  use  as  input  for  HyGAS  functional  testing. 

3.1.4  Task  1  Summary 

HyGAS  incorporates  methods  of  gas/vapor  cloud  detection  and  gas  identification  using 
primarily  spectral  rather  than  image  (spatial)  analysis  techniques. 

Gas  detection  in  HyGAS  can  only  be  done  spectrally  (on  raw  interferograms  using 
Kroutil’s  linear  discriminant  test).  Gas  detection  using  spatial  techniques  (on 
multispectral  bands,  using  Althouse’ s  statistical  concentration  model  and  a  recursive 
Kalman  filtering  technique)  should  be  kept  under  consideration  for  future 
implementations. 

Gas  identification  can  only  be  done  spectrally  (on  transformed  spectra,  using  Rauss’s 
trained  search  technique).  This  technique  is  incorporated  into  HyGAS. 
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All  of  these  methods  are  based  heavily  on  statistical  probabilities,  and  meaningful  image 
display  and  interpretation  may  require  additional  information.  One  suggestion  of  the 
COTR  is  to  allow  for  low  spatial  resolution  analysis  results  to  be  superimposed  on  higher 
resolution  image  data.  This  is  a  technique  that  should  be  considered  for  future 
implementation. 

3.2  Task  2 — Software  Development 

HyGAS  was  designed  to  be  part  of  the  existing  HIPS  software  package.  This  enables 
HyGAS  to  access  all  of  the  existing  HIPS  math  functions,  display  routines,  libraries  and 
user  interfaces.  Hardware  and  software  issues  required  compatiblity  with  the  existing 
platform. 

3.2.1  Hardware  acquisition  and  basic  system  platform  preparation 

A  SUN  SPARCstation  10  was  selected  and  purchased  (GFE)  as  the  HyGAS  software 
system  platform.  This  hardware  was  selected  for  compatibility  with  the  existing  HIPS 
platform,  but  with  significantly  increased  speed  to  enable  faster  processing.  This 
hardware  initially  arrived  with  the  incorrect  operating  system,  which  caused  some  delays 
in  getting  the  hardware  operational.  This  difficulty  was  resolved  when  the  correct 
software  was  supplied  by  the  manufacturer. 

HyGAS  was  designed  to  function  as  a  module  of  the  existing  HIPS  software  system, 
which  was  developed  on  a  SPARCstation  1  under  the  Sun  O/S  4. 1  operating  system, 
using  the  original  Kernighan  and  Ritchie  C-compiler.  The  SPARCstation  10  runs  the 
Solaris  2.3  operating  system  and  uses  the  more  universal  standard  ANSI  C-compiler. 

This  required  all  existing  HIPS  code  to  be  modified  to  conform  to  ANSI  C  standards 
before  it  would  compile  properly  on  the  SPARCstation  10.  This  process  involved 
changing  all  function  declarations  from  the  Kernighan  and  Ritchie  standard  to  ANSI  C 
standard.  Once  this  was  completed,  all  HIPS  routines  had  to  be  tested  for  any  run-time 
errors  caused  by  the  new  compiler. 

The  existing  HIPS  source  code  was  transferred  from  the  SPARCstation  1  onto  the 
SPARCstation  10,  and  was  converted  to  ANSI  C.  Most,  but  not  all  HyGAS  HIPS 
(version  1.7)  routines  have  been  tested  for  run-time  errors.  Some  HIPS  modules  may  still 
contain  run-time  errors. 

3.2.2  HyGAS  software  development 

HyGAS  software  development  concentrated  on  obtaining  existing  algorithms  and  code 
from  other  parties  and  manufacturers  for  incorporation  into  Hy GAS/HIPS.  Three 
primary  efforts  were  made:  (1)  acquisition  of  suitable  contour  plotting  and  3-D  display 
software;  (2)  acquisition  of  suitable  library  search  routines  and/or  detailed  gas  analysis 
software;  and  (3)  acquisition  of  specialized  source  code,  such  as  discriminant  testing,  etc. 

3.2.2. 1  Contour  plotting  and  3-D  display  software 

Three  commercially  available  software  packages  were  acquired  for  testing  on  a  trial 
basis:  (1)  S-Plus/Stat  Sci;  (2)  NCAR;  and  (3)  XRT/3d.  None  of  these  packages  proved 
suitable  for  incorporation  into  HyG AS/HIPS.  The  S-Plus  package  did  not  provide 
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suitable  results  (we  want  labeled/identified  contour  levels),  and  the  separate  software 
library  modules  required  to  create  a  run-time  package  summed  up  to  $5K.  The  NCAR 
package  produced  nice  output,  but  was  very  slow.  Also,  this  package  produces  a 
metafile,  i.e.  a  computer  generated  graphics  code  file  that  with  an  appropriate  driver  can 
be  plotted  on  the  monitor  screen.  However,  there  is  no  way  to  overlay  this  file  with  a 
HIPS  image  file,  and  no  way  to  access  this  file  with  any  HIPS  function.  The  XRT/3d 
package  had  severe  limitations  on  the  number  of  pixels  it  could  handle,  and  it  was  very 
slow. 

Given  these  difficulties,  we  elected  not  to  pursue  3-D  plotting  for  HyGAS. 

Two  options  seem  feasible  for  future  consideration:  the  Interactive  Data  Analysis 
Language  (IDL);  and  the  Advanced  Visualization  System  (AVS).  Both  of  these  packages 
are  combined  image  processing  and  plotting  packages,  and  ought  to  be  able  to  handle 
contour  plots  and  image  overlays.  (The  three  prior  packages  were  just  plotting  packages, 
not  imaging  packages.)  We  anticipate  that  one  or  both  of  these  packages  should  be  able 
to  do  both  image  and  contour  overlays,  but  we  would  need  to  find  out  if  the  data  formats 
are  accessible  to  HIPS.  They  both  have  accessible  software/display  libraries,  so  they  may 
prove  suitable  on  both  counts. 

3.2.2.2  Library  search  routines  and  detailed  gas  analysis  software 

We  purchased  a  PC-compatible  software  package  for  laboratory  gas  analysis  called 
GRAMS/386  (previously  referred  to  as  SpectraCalc).  Although  this  package  is  not  UNIX 
compatible,  we  decided  that  since  it  was  strongly  suggested  by  both  ERDEC  and 
NVEOD,  as  well  as  other  contacts  (and  referred  to  as  the  unofficial  “Army  standard”  for 
spectral  analysis),  we  would  acquire  it  for  possible  interfacing  with  HyGAS.  SpectraCalc 
is  designed  for  use  with  laboratory  spectrometers,  primarily  for  FTS  instruments.  As 
such,  it  is  designed  to  deal  with  high  spectral  resolution  data.  Its  search  routines  are 
based  on  the  premise  that  laboratory  collected  data  will  generally  be  at  higher  spectral 
resolution  than  library  data,  and  matching  the  laboratory  spectral  resolution  of  the 
unknown  agent  to  the  standard  resolution  of  the  spectral  library  will  be  a  process  of  high- 
to-low  convolution  (such  that  convolution  can  be  done  on  the  unknown  spectrum  to 
match  the  lower  resolution  library),  or,  if  low-to-high  convolution,  not  at  significantly 
different  resolutions  (such  that  either  linear  or  spline  interpolation  would  not  introduce 
significant  discrepancies).  This  situation  is  not  likely  to  be  attained  for  imaging 
spectrometers  in  the  near  future,  but  the  capability  should  not  be  dismissed.  HyGAS  was 
not  designed  to  interface  with  SpectraCalc  (as  originally  intended),  to  be  able  to  extract 
spectra  from  a  HIPS  image  and  store  them  in  standard  SpectraCalc  and/or  JCAMP 
format.  However,  HyGAS  was  designed  to  export  ASCII_XY  data,  so  that  any  PC 
networked  to  the  SPARCstation  should  be  able  to  import  these  spectra  directly  to 
SpectraCalc  for  traditional  gas  analysis. 

We  researched  the  availability  of  traditional  gas  phase  spectral  libraries  for  use  with 
Hy GAS/HIPS.  Traditional  gas  phase  spectral  libraries  are  available  in  PC-compatible 
format  from  several  different  vendors  (Aldrich  Chemical,  Sadtler  Research,  Sprouse 
Scientific,  Infrared  Analysis,  etc.).  Many  of  these  vendors  supply  copyrighted  and/or 
encrypted  libraries  that  can  only  be  accessed  with  specific  PC  software  packages.  This 
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means  HIPS  cannot  access  these  libraries,  nor  can  we  write  translation  routines. 
SpectraCalc  can  access  all  of  these  libraries  (as  well  as  any  other  library  in  “Galactic” 
format).  [Galactic  Industries  is  the  manufacturer  of  the  GRAMS/386  and  the  related 
SpectraCalc  and  LabCalc  programs.]  Very  few  of  these  libraries  are  available  in 
SPARC/UNIX  compatible  format;  a  few  are  available  in  JCAMP  format,  for  an 
additional  fee.  We  approached  Aldrich  Chemical  about  obtaining  their  library(ies)  in 
SPARC/UNIX  compatible  format,  but  were  unable  to  reach  a  solution.  Hence,  for 
HyGAS  we  were  unable  to  provide  either  search  routines  or  access  routines  for  such 
libraries. 

We  approached  Aldrich  about  acquiring  data  access  rights,  perhaps  by  signing  a  non¬ 
disclosure  agreement  or  by  providing  encryption  protection  within  Hy GAS/HIPS.  The 
difficulty  is  that  for  HyGAS,  we  must  provide  for  data  collected  at  lower  spectral 
resolution  than  that  of  the  libraries  (perhaps  at  substantially  lower  resolution,  such  that 
linear  or  spline  interpolation  methods  would  not  be  applicable).  For  encrypted  and/or 
copyrighted  libraries,  this  will  only  be  possible  if  we  can  obtain  the  approval  of  the 
library  supplier  to  not  only  access  the  library,  but  also  to  convolve  it  to  lower  resolutions 
and  store  the  convolved  library  as  a  HIPS  library.  In  the  future,  we  should  consider 
encryption  of  HIPS  libraries,  so  that  appropriate  convolution  of  the  library  to  the 
HyG AS/HIPS  sensor  can  be  performed  without  infringing  on  the  manufacturer’s 
restrictions.  In  the  unlikely  situation  of  HyGAS/HIPS  resolution  higher  than  that  of  a 
library,  it  is  a  simple  matter  to  convolve  the  higher  resolution  HyGAS  spectrum  to  the 
lower  resolution  of  the  library  before  performing  the  spectral  search. 

We  acquired  the  novel  search  method  using  backpropagation  neural  nets  developed  by 
Rauss.  Details  of  this  approach  are  described  in  Appendix  D. 

3.2.23  Specialized  source  code 

Specialized  source  code  to  perform  the  linear  discriminant  test  was  acquired  from  Kroutil 
in  PC-compatible  code.  This  code  was  transferred  to  the  SPARCstation  10,  and  rewritten 
to  conform  to  standard  ANSI  C.  Graphics  libraries  accessed  by  the  PC  code  were  not 
available  for  the  SPARC  hardware,  and  a  few  required  PC  routines  were  rewritten  to  call 
standard  X-window  libraries. 

Specialized  code  to  actually  determine  the  linear  discriminant  is  available  from  Kroutil, 
but  only  in  SiliconGraphics  (SGI)  C-code.  This  code  accesses  many  of  the  specialized 
graphics  display  features  available  only  on  an  SGI  machine.  For  the  future,  we  should 
consider  acquiring  this  code  and  a  low-end  SGI  machine  to  run  it  on.  (Estimated  GFE 
price  for  such  a  machine  is  ~$5K.)  The  intent  would  be  to  allow  the  user  to  select  several 
image  cube  datasets,  and  run  these  datasets  through  the  special  pattern  recognition 
algorithms  developed  by  Kroutil  to  determine  a  less  robust  discriminant  vector,  but  one 
that  would  likely  be  usable  for  subsequent  datasets.  [For  example,  Kroutil  has  spent 
months  of  CPU  time  determining  a  robust  discriminant  vector  for  one  specific  gas  based 
on  several  hundred  thousand  spectra  acquired  with  a  point  spectrometer  (MID AC  FTS 
sensor).  With  a  set  of  hyperspectral  images,  perhaps  as  few  as  5  images  carefully 
selected  and  defined  would  produce  enough  data  to  define  a  reasonably  efficient 
discriminant.  This  discriminant  could  be  determined  off-line,  using  the  SGI  machine,  and 
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the  resultant  discriminant  could  then  be  added  to  the  library  on  the  HyGAS  machine  for 
consideration  in  subsequent  datasets.] 

Specialized  code  to  do  image  enhancement  (“edge  detection”  for  vapor  clouds,  and 
background  clutter  suppression)  could  be  obtained  from  M.  Althouse.  We  intended  to 
incorporate  these  techniques  into  HyGAS,  but  were  unable  to  do  so. 

3.2.3  Task  2  Summary 

HyGAS  appears  under  a  Gas  Analysis  menu  within  HIPS.  This  menu  contains  the 
individual  gas  analysis  functions,  which  can  be  run  non-interactively,  as  well  as  the 
HyGAS  module  which  incorporates  them  into  an  interactive  package.  Current  functions 
include  Linear  Discriminant  Test  (gas  detection),  Background  Removal  (subtraction  or 
division),  Fast  Fourier  Transform  (interferogram  to  spectrum)  and  reverse,  and  Library 
Search  and  Neural  Net  search  (gas  identification/analysis). 

Analysis  can  be  performed  on  raw  instrumental  interferograms  or  on  transformed  (or  raw 
instrumental)  spectra.  Techniques  developed  by  Kroutil  to  do  gas  detection  by  linear 
discriminant  tests  are  performed  in  interferogram  space.  Spectral  search  techniques 
developed  by  Rauss  are  performed  in  spectral  space.  Data  from  an  FTS  may  be 
processed  through  either  or  both  techniques. 

3.3  Task  3 — Spectral  Database/Measurements  Development 

Simulation  of  a  test  dataset  required  a  significant  effort  in  HyGAS.  We  acquired  some 
data  from  R.  Kroutil  and  from  P.  Rauss,  but  without  much  descriptive  information. 

We  considered  J.  Salisbury  (of  Johns  Hopkins  University)  as  a  source  for  alternative 
thermal  IR  data,  along  with  the  possibility  of  renting  a  MID  AC  sensor  to  measure  actual 
interferograms  of  realistic  backgrounds  to  match  Kroutil’s  data  for  use  with  his 
predetermined  discriminants.  Neither  of  these  options  was  successful. 

Full  details  on  GFI  databases  and  libraries,  and  a  description  of  the  use  of  these  datasets 
to  develop  simulated  test  image  cubes,  are  found  in  the  appendices. 

3.4  Task  4 — Documentation 

Since  HyGAS  was  developed  as  a  functional  part  of  HIPS,  the  software  documentation, 
or  user’s  manual,  for  HyGAS  was  expected  to  become  part  of  the  HIPS  User’s  Manual. 
This  was  the  plan  when  the  interim  software  documentation  was  set  up  as  Chapter  12- 
Gas  Analysis,  in  the  standard  format  of  the  HIPS  User’s  Manual,  with  preliminary 
function  and  parameter  specification  descriptions.  However,  as  HyGAS  progressed,  we 
elected  to  create  the  User's  Manual  in  the  form  of  a  tutorial,  and  we  delivered  the  HyGAS 
User's  Manual/Tutorial  instead  of  the  reference  manual  format  of  the  HIPS  User  s 
Manual. 

Complete  details  on  the  use  of  HyGAS  with  the  provided  and  the  simulated  datasets  is 
given  in  the  HyGAS  User's  Manual/Tutorial. 
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4.0  Additional  Information 

Additional  details  on  GFI  datasets,  GFI  databases,  creation  of  simulated  image  cubes, 
neural  net  and  linear  discriminants,  are  provided  in  the  appendices. 
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HyGAS  Field  Instrument  Datasets  (GFI) 


HyGAS  Field  Instrument  Datasets  (GFI) 

Two  types  of  field  instrument  datasets  were  supplied  (GFI)  for  use  with  HyGAS. 

•  RaussJHowden  spectral  datasets 

•  Kroutil_MIDAC  interferogram  datasets 


1.  Rauss_Howden  spectral  datasets 

The  Rauss_Howden  spectral  datasets  came  from  a  field  experiment  by  E.  Howden  of 
NVEOD.  These  measurements  (11  and  12  August  1992)  were  collected  using  a  point 
sensor  FTIR  spectrometer.  Individual  diskette  listings  are  included  below.  Figures  1  and 
2  show  associated  logbook  notes.  Table  1  lists  the  measurements  of  the  better  dataset  (12 
August)  in  chronological  order;  corresponding  spectra  are  shown  in  Figures  3-5. 


Rauss_Howden  Field  Data,  Disk  1  Listing. 

Volume  in  drive  A  is  FIELDDATA 
Volume  Serial  Number  is  2C63-13F9 
Directory  of  A:\ 


ACEl 

SPC 

365,824 

08-12-92 

2  : 30p 

ACE2 

SPC 

365,824 

08-12-92 

2  : 38p 

ACE3 

SPC 

365, 824 

08-12-92 

2  : 21p 

BG1 

SPC 

2,052 

08-12-92 

3 : 39p 

BG2 

SPC 

36,784 

08-12-92 

2 : 08p 

BG3 

SPC 

36,784 

08-12-92 

3  : 40p 

BG4 

SPC 

36,784 

08-12-92 

4 : 3  lp 

README 

README 

BAR 

192 

1,028 

04-15-94 

04-15-94 

1 : 04p 
2  :  54p 

9  f ile (s )  1,211,096  bytes 


244,224  bytes  free 


Rauss_Howden  Field  Data,  Disk  2  Listing. 

Volume  in  drive  A  is  FIELDDATA2 
Volume  Serial  Number  is  283A-14D0 
Directory  of  A:\ 

MEK1  SPC  914,224  08-12-92  4 : 57p 

MEK2  SPC  365,824  08-12-92  5:05p 

2  file (s)  1,280,048  bytes 

_ 177,152  bytes  free _ 


Rauss_Howden  Field  Data,  Disk  3  Listing. 

Volume  in  drive  A  is  FIELDDATA3 
Volume  Serial  Number  is  0E2D-14E9 


Directory  of 

A:  \ 

81292 

<DIR> 

04-15-94 

1 : 40p 

81192 

<DIR> 

04-15-94 

1 :44p 

2 

!  file(s) 

0  bytes 

Directory  of 

A:\81192 

<DIR> 

04-15-94 

1 : 44p 

<DIR> 

04-15-94 

1 : 44p 

BG10 

SPC 

36,784 

08-11-92 

1 : 45p 

BG11 

SPC 

36,784 

08-11-92 

2  : 03p 

BG12 

SPC 

36,784 

08-11-92 

2  : 50p 

BG13 

SPC 

91,624 

08-11-92 

3  : 21p 

BG14 

SPC 

91,624 

08-11-92 

6  : 14p 

BG15 

SPC 

91,624 

08-11-92 

6 : 22p 

BGROUND4 

SPC 

5,708 

08-11-92 

10:22a 

BGROUND5 

SPC 

36,784 

08-11-92 

10:49a 

BGROUND6 

SPC 

36,784 

08-11-92 

11:00a 

BGROUND7 

SPC 

36,784 

08-11-92 

11:02a 

BGROUND8 

SPC 

36,784 

08-11-92 

11:06a 

BG ROUND 9 

SPC 

36,784 

08-11-92 

11:19a 

SF62 

SPC 

67,860 

08-11-92 

11:34a 

15 

filets)  642,712  bytes 

Directory 

'  of 

A:\81292 

<DIR> 

04-15-94 

1 : 40p 

<DIR> 

04-15-94 

1 : 4  Op 

SF61 

SPC 

365,824 

08-12-92 

3  :  59p 

SF62 

SPC 

365,824 

08-12-92 

4  :  07p 

README 

63 

04-15-94 

1 : 42p 

5  file(s)  731,711  bytes 


Total  files  listed: 

22  filets)  1,374,423  bytes 
_ 79,872  bytes  free 


Rauss_Howden  Field  Data,  Disk  4  Listing. 

Volume  in  drive  A  is  FIELDDATA 4 
Volume  Serial  Number  is  0C3C-14F0 
Directory  of  A:\ 

ACE1  SPC  914,224  08-11-92  3:52p 

ACE2  SPC  365,824  08-11-92  4:00p 

2  filets)  1,280,048  bytes 
_ 177,152  bytes  free _ 


RaussJHowden  Field  Data,  Disk  5  Listing. 

Volume  in  drive  A  is  FIELDDATA5 
Volume  Serial  Number  is  0410-15CA 


Directory  of  A: 

:  \ 

ACE3 

SPC 

365,824 

08-11-92 

4  : 08p 

ACE4 

SPC 

365,824 

08-11-92 

4  : 16p 

ACE5 

SPC 

365,824 

08-11-92 

4  : 29p 

SF61 

SPC 

183,024 

08-11-92 

11:32a 

4  filets)  1,280,496  bytes 
_ 176,128  bytes  free 


RaussJHowden  Field  Data,  Disk  6  Listing. 

Volume  in  drive  A  is  FIELDDATA6 
Volume  Serial  Number  is  203A-15CE 
Directory  of  A:\ 

SF63  SPC  914,224  08-11-92  2:35p 

1  filets)  914,224  bytes 

543,232  bytes  free 


Volume  in  drive  A  is  FIELDDATA 7 
Volume  Serial  Number  is  2F18-15D1 
Directory  of  A:\ 

MEK1  SPC  914,224  08-11-92  7:00p 

1  filets)  914,224  bytes 

_ 543,232  bytes  free 


Table  1.  Summary  of  Rauss_Howden  field 
data  from  12  August  1992.  The  datasets  listed 
here  (in  chronological  order)  were  identified  as 
preferable  to  those  from  1 1  August  1992. 
Corresponding  spectral  plots  for  these  datasets 
are  shown  in  Figures  3-5. 


BG2 

SPC 

36,784 

08-12-92 

2 : 08p 

ACE3 

SPC 

365,824 

08-12-92 

2 : 21p 

ACEl 

SPC 

365,824 

08-12-92 

2  : 30p 

ACE2 

SPC 

365,824 

08-12-92 

2  : 38p 

BG1 

SPC 

2,052 

08-12-92 

3  : 39p 

BG3 

SPC 

36,784 

08-12-92 

3  : 40p 

SF61 

SPC 

365,824 

08-12-92 

3  : 59p 

SF62 

SPC 

365,824 

08-12-92 

4  : 07p 

BG4 

SPC 

36,784 

08-12-92 

4  : 3  lp 

MEK1 

SPC 

914,224 

08-12-92 

4  :  57p 

MEK2 

SPC 

365,824 

08-12-92 

5  : 05p 
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Figure  1.  Logbook  notes  for  Rauss_Howden  experimental  gas  release  on  12  August  1992. 

Annotations  indicate  this  dataset  is  preferable  to  that  of  1 1  August  1992. 
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Figure  2.  Logbook  notes  for  Rauss_Howden  experimental  gas  release  on  11  August  1992. 

These  field  measurements  cover  the  spectral  range  833-1265  cm'1  at  2  cm'1  resolution 
(449  datapoints  each,  Nyquist  sampling  at  1  cm'1  datapoint  spacing).  [Presumably,  these 
spectra  were  originally  collected  as  interferograms  by  an  FTIR  spectrometer,  were 
subsequently  transformed  into  the  spectral  domain,  and  then  were  trimmed  to  exclude 
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datapoints  outside  the  given  spectral  range.]  Laboratory  absorbance  spectra  (at  a  coarser 
resolution  of  12  cm'1  resolution)  for  these  same  gases  are  shown  in  Figure  6. 

RaussJHowden  Field  Data 


Wavenumber  (inverse  cm) 

Figure  3.  Spectral  measurements  for  acetone  and  associated  field  backgrounds  from  12  August  1992. 

Absorptions  due  to  acetone  are  clearly  visible  in  sample  acei. 

Rauss_Howden  Field  Data 


Wavenumber  (inverse  cm) 

Figure  4.  Spectral  measurements  for  sulfur  hexafluoride  (SF<0  and  associated  field  backgrounds 
from  12  August  1992.  Absorptions  due  to  SF6  are  barely  discernible  in  sample  sf  62.  [Refer  to  Figure  4 
for  corresponding  laboratory  absorbance  spectra.] 
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RaussJHowden  Field  Data 

12  August  1992,  MEK  test  release 


Wavenumber  (inverse  cm) 

Figure  5.  Spectral  measurements  for  methyl-ethyl-ketone  (MEK)  and  associated  field  backgrounds 
from  12  August  1992.  Absorptions  due  to  MEK  are  clearly  visible  in  sample  mek2.  However,  there  is  a 
significant  difference. in  the  overall  radiance  for  this  sample  and  the  corresponding  background,  which  may 
be  caused  by  possible  alterations  in  the  instrumental  gain  settings. 


Spectral  Signature  Comparison 


(from  NIST/EPA  Gas  Phase  Database) 


Wavenumber  (inverse  cm) 

Figure  6.  Laboratory  absorbance  spectra  for  acetone,  sulfur  hexafluoride  and  methyl-ethyl-ketone. 

These  absorbance  spectra  came  from  the  NIST/EPA  Database  35,  Gas  Phase  Infrared,  and  are  measured  at 
a  significantly  coarser  spectral  resolution  (12  cm"1  spectral  resolution)  than  the  field  data  (at  2  cm" ). 
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2.  Kroutil_MIDAC  Interferogram  Datasets 

The  Kroutil_MIDAC  interferogram  datasets  came  from  both  field  and  laboratory 
experiments  by  R.  Kroutil  of  ERDEC.  All  measurements  [presumably]  were  collected 
using  a  MIDAC  point  sensor  FTIR  spectrometer,  operated  at  8  cur1  spectral  resolution. 
This  mode  of  operation  generates  interferograms  with  1024  datapoints,  which  correspond 
to  spectra  of  512  datapoints  over  the  spectral  range  0  to  1971  cm'1  wavenumbers. 

Three  interferogram  datasets  were  supplied. 

•  mdtOl 

•  kft0002.dat 

•  kft0048.dat 

Dataset  mdtOl  is  [presumably]  a  series  of  field  interferograms  collected  from  a  helicopter 
overflight  during  a  ground  test  release  of  sulfur  hexafluoride  (SF6).  No  additional  details 
were  supplied  with  this  dataset.  Unfortunately,  since  no  discriminant  coefficients  to 
detect  the  presence  of  SF6  were  supplied  with  this  dataset,  we  were  unable  to  determine 
which  records  within  this  file  included  SF6  and  which  were  field  background  only. 

Dataset  kft0002.dat  is  [presumably]  a  collection  of  interferograms  measured  in  a 
laboratory  under  controlled  conditions,  using  a  MIDAC  spectrometer,  a  gas  cell,  and  a 
laboratory  thermal  blackbody  source.  Various  conditions  (background  source 
temperature,  gas  temperature,  gas  concentration)  are  presumably  represented  for 
situations  with,  and  without,  the  gas  MEK  in  the  optical  path.  Discriminant  coefficients 
for  both  acetone  and  MEK  were  supplied  with  this  dataset.  [However,  no  measurements 
of  acetone  were  detected  within  this  dataset.]  No  additional  details  were  supplied  with 
this  dataset. 

Dataset  kft0048.dat  is  [presumably]  a  collection  of  laboratory  background  interferograms 
measured  under  various  conditions,  but  all  without  the  gas  MEK.  No  additional  details 
were  supplied  with  this  dataset. 

Example  interferograms  from  each  dataset  are  shown  superposed  in  Figure  7.  Examples 
expanded  about  the  centerburst  region  are  shown  in  Figure  8.  All  interferograms  from 
dataset  mdtOl  appear  nearly  identical,  suggesting  a  uniform  field  background  differing 
only  in  the  presence  or  absence  of  SF6  in  the  optical  path.  Interferograms  from  the  other 
two  datasets  exhibit  comparable  variability  within  distinct  groups  of  contiguous 
measurements,  but  vary  between  groups.  This  suggests  differences  in  the  conditions 
under  which  each  group  of  interferograms  was  collected,  such  as,  perhaps,  different 
temperatures  for  the  background  thermal  source  in  each  group.  Details  describing  the 
different  groups  of  measurements  in  each  dataset  were  not  available. 

Individual  diskette  listings  of  the  datasets  (and  programs)  supplied  by  Kroutil  are 
included  below.  Listings  of  the  ERDEC  dataset  diskettes  also  are  included. 
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Interferogram  Datasets  from  Kroutil 
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Figure  7.  Example  interferograms  from  datasets  mtdoi,  kft0002.dat,  and  kftoo48.dat.  Overall 
intensity  is  comparable  to  the  example  shown  for  all  records  within  mtdoi,  but  varies  between  groups  of 
comparable  data  for  the  other  datasets. 


(example  of  first  10  records  in  each  file) 
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Interferogram  Datasets  from  Kroutil 


Interferogram  Channel 

Figure  8.  Example  interferograms  from  datasets  mtdoi,  kft0002.dat,  and  kft0048.dat,  region  near 
centerburst,  only.  The  centerburst  occurs  at  channel  1 16  for  interferograms  in  mtdoi,  at  channel  48  for 
kf  tooo2 .  dat,  and  at  channel  84  (a  negative  burst)  for  kf  too48 .  dat.  Example  interfrograms  from  the  first 
two  groups  of  interferograms  in  kft0002  .dat  are  shown. 
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Kroutil  Disk  1,  Collect  Programs. 

Volume  in  drive  A  has  no  label 


Directory  of  A:\ 


COLLECT 

<DIR> 

07-15-94 

12 : 32p 

1  file ( s) 

0  bytes 

Directory  of 

A : \ COLLECT 

<DIR> 

07-15-94 

12  : 32p 

<DIR> 

07-15-94 

12 : 32p 

MIDCOL 

EXE 

54,173 

12-29-93 

12 : 20p 

MIDCOL 

VI 

51,025 

09-25-92 

3  :  lOp 

REPLAY 

EXE 

42,394 

05-03-93 

8:51a 

HEADERST 

DEF 

2,505 

05-04-92 

9  : 34p 

REPLAY 

C 

29,299 

05-01-93 

3  : 53p 

OSCCONV 

C 

9.113 

04-26-93 

8:54a 

ASM 

BAT 

745 

05-16-92 

7 : 57p 

MIDCOL 

BAT 

135 

08-15-92 

9 :  07p 

MIDDEF 

H 

2,816 

07-24-92 

7:13a 

MIDCOLD 

C 

31,783 

08-25-92 

3  : 44p 

MIDCOLD 

BAT 

138 

08-25-92 

3  : 46p 

MIDCOLD 

EXE 

46,863 

10-04-92 

3  : 41p 

HEADERS 

DEF 

1,578 

09-27-92 

2 : 36p 

REPLAY 

VI 

23,775 

11-05-92 

11:46a 

REPLAY 

BAT 

132 

05-01-93 

11:47a 

INSTALL 

BAT 

267 

10-13-92 

2  : 50p 

MIDLAB 

C 

68,692 

02-17-93 

8:26a 

MIDLAB 

EXE 

54,019 

02-17-93 

8:22a 

MIDCOL 

V2 

61,271 

11-05-92 

11:51a 

MIDLAB 

BAT 

135 

01-26-93 

2 : 56p 

MIDLAB 

VI 

61,752 

02-05-93 

9:51a 

F5 

FSP 

2,072 

02-13-93 

7  :  51p 

CONVINTF 

BAT 

137 

04-27-93 

9:01a 

F6 

FSP 

2,072 

02-13-93 

7 : 51p 

F7 

FSP 

2,072 

02-13-93 

7  : 51p 

F8 

FSP 

2,072 

02-13-93 

7 : 52p 

REPLAY 

V2 

26,068 

02-12-93 

8 : 48p 

SCCONV 

C 

14,719 

04-30-93 

9:22a 

EXSCCONV 

DEF 

381 

04-27-93 

10:01a 

SCCONV 

BAT 

99 

04-26-93 

8:54a 

SCALC 

DEF 

970 

04-02-93 

6 : 51p 

CONVINTF 

EXE 

239,818 

04-27-93 

1 : 37p 

SCCONV 

EXE 

175,742 

04-30-93 

9:22a 

MENU 

H 

2,084 

01-18-90 

6  : 46p 

EXCONV 

DEF 

402 

04-27-93 

8:50a 

EXREPLAY 

DEF 

478 

05-01-93 

11:56a 

CONVINTF 

C 

11,915 

05-03-93 

9:07a 

MIDCOLV 

BAT 

180 

05-04-93 

6:08a 

MENU 

C 

11,208 

01-18-90 

6 : 44p 

EXMIDCOL 

DEF 

359 

05-03-93 

8 : 22p 

MIDCOLV 

EXE 

190,247 

12-29-93 

12 : 32p 

MIDCOL 

C 

68,434 

12-29-93 

12 : 19p 

MIDCOLV 

c 

72,357 

12-29-93 

12 : 31p 

45  file(s)  1, 366, 496  bytes 


Total  files  listed: 

46  file ( s)  1,366,496  bytes 
_ 77,824  bytes  free 


Kroutil  Disk  2,  Pattern  Recognition 
Programs. 


1  Volume  in  drive  A 

has  no 

label 

1  Volume 

Serial  Number  is 

2A4A-0FFC 

Directory  of 

A :  \ 

MTRX 

<DIR> 

07-15-94 

12 : 3  5p 

1  filets) 

0  bytes 

Directory  of 

A:\MTRX 

<DIR> 

07-15-94 

12 : 3  5p 

<DIR> 

07-15-94 

12 : 35p 

MTRX 

DEF 

318 

08-17-92 

12:08a 

MIDDEF 

H 

2,816 

07-24-92 

7:13a 

AMTRX 

EXE 

52,843 

04-15-93 

9:52a 

MTRXD 

BAT 

132 

08-17-92 

1 : 44p 

AMTRXD 

EXE 

49,403 

04-15-93 

9:53a 

MTRX 

BAT 

129 

08-17-92 

1 : 46p 

HEADERS 

DEF 

1,578 

09-27-92 

2:36p 

FILTER1 

INC 

28,285 

04-02-93 

2  : 01p 

MATGRAPH 

EXE 

45,034 

10-02-92 

3  : 52p 

MATGRAPH 

BAT 

112 

10-01-92 

10:45a 

MTRXD 

C 

16,243 

02-17-93 

10:19a 

DISCRIM1 

INC 

5,122 

04-03-93 

9  :  07a 

MTRX 

EXE 

56,035 

12-27-93 

7 : 03p 

MTRX2 

EXE 

65,443 

12-27-93 

7 : 02p 

MTRXD 

EXE 

52,595 

04-15-93 

9:55a 

MTRX2 

BAT 

132 

02-09-93 

3  : 37p 

MATGRAPH 

C 

5,516 

02-17-93 

10:18a 

FILTER2 

INC 

18,421 

04-02-93 

2  :21p 

DISCRIM2 

INC 

5,  122 

04-03-93 

9:06a 

MFILTER 

INC 

28,285 

04-02-93 

2 :  Olp 

AFILTERl 

INC 

18,421 

04-02-93 

2 : 20p 

AFILTER2 

INC 

18,421 

04-02-93 

2:21p 

ADIS1 

INC 

5,122 

04-03-93 

9:04a 

ADIS2 

INC 

5,122 

04-03-93 

9:06a 

MDIS1 

INC 

5,122 

04-03-93 

9:07a 
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MTRX  2 

C 

42,995  12-27-93  7:00p 

MTRX 

C 

40,406  12-27-93  7 : OOp 

29  file (s) 

569,173  bytes  ; 

Total 

files  listed: 

30  filets) 

569,173  bytes 

880,128  bytes  free 

Kroutil  Disk  3,  Datafile  #1,  kft0002.dat. 


Volume  in  drive  A  has  no  label 
Volume  Serial  Number  is  192A-15F0 

Directory  of  A:\ 

KFT0002  DAT  587,648  07-15-94  12:25p 

1  file(s)  587,648  bytes 

Total  files  listed: 

1  file (s)  587,648  bytes 

_ 869,888  bytes  free 


Kroutil  Disk  4,  Datafile  #2,  kft0048.dat. 


Volume  in  drive  A  has  no  label 
Volume  Serial  Number  is  3A15-1AEC 

Directory  of  A:\ 

KFT0048  DAT  1,162,112  07-15-94  8:21a 

1  file (s )  1,162,112  bytes 

Total  files  listed: 

1  file(s)  1,162,112  bytes 

_ 274,432  bytes  free 


Kroutil  Disk  5,  FTIR  data  collection 
software,  version  2,  update  1, 3/17/93. 


Volume 

in  drive 

has  no  label 

Volume 

Serial  Number  is 

2A3C-12E7 

Directory  of 

A:  \ 

COLLECT 

<DIR> 

07-06-88 

4:50a 

PC  I  DA 

<DIR> 

07-06-88 

4:50a 

DATA 

<DIR> 

07-06-88 

4:51a 

INSTALL 

BAT 

257 

02-28-93 

6 : 22p 

READ 

ME 

1,306 

02-28-93 

6 : 25p 

5  filets) 

1,563  bytes 

1  Directory  of 

A :  \ COLLECT 

<DIR> 

07-06-88 

4:50a 

<DIR> 

07-06-88 

4:50a 

MIDCOL 

EXE 

54,091 

02-17-93 

5  : 26p 

REPLAY 

C 

29,155 

02-28-93 

5  : 26p 

REPLAY 

EXE 

43,057 

02-28-93 

5  : 17p 

HEADERST 

DEF 

2,505 

05-04-92 

9  : 34p 

REPLAY 

BAT 

135 

08-25-92 

3  : 48p 

CONVINTF 

42,852 

04-26-93 

8:56a 

MIDCOL 

BAT 

135 

08-15-92 

9 :  07p 

MIDDEF 

H 

2,816 

07-24-92 

7:13a 

READ 

ME 

463 

02-25-93 

2  : 45p 

MIDCOL 

C 

68,614 

04-06-93 

9:02a 

HEADERS 

DEF 

1,578 

09-27-92 

2  : 36p 

MIDDOC 

TXT 

21,479 

10-13-92 

2  : 3  8p 

MIDDOC 

WP 

21,600 

10-13-92 

2  : 3  8p 

F5 

FSP 

2,072 

02-13-93 

7  :  51p 

F6 

FSP 

2,072 

02-13-93 

7 : 51p 

F7 

FSP 

2,072 

02-13-93 

7 : 51p 

F8 

FSP 

2,072 

02-13-93 

7 : 52p 

SCALC 

DEF 

970 

04-02-93 

6  :  51p 

20  filets) 

297,738  bytes 

Directory  of 

A : \ DATA 

<DIR> 

07-06-88 

„  4 ; 51a 

<DIR> 

07-06-88 

4:51a 

MDT01 

361,664 

02-20-93 

5  : 29p 

READ 

ME 

288 

02-28-93 

5:  lip 

4  filets) 

361,952  bytes  j 

Directory  of 

A : \ PC IDA 

<DIR> 

07-06-88 

4:50a 

<DIR> 

07-06-88 

4:50a 

PC  I  DA 

C 

9,793 

03-24-93 

10:08a 

TEMPS 

C 

14,472 

07-22-92 

2  : 18p 

RMATH 

c 

18,771 

07-23-92 

10 : 44p 

DEFINES 

IDA 

375 

04-20-90 

12 : 38p 

EXTERNS 

IDA 

525 

04-03-90 

8:31a 

HEADERS 

IDA 

1,504 

07-20-90 

9:58a 

PC  IDA 

DOC 

3,061 

12-06-90 

10:02a 

INTF_TMP 

IDA 

40,960 

07-24-92 

4  :  51p 

SPEC_TMP 

IDA 

20,480 

07-24-92 

4  : 51p 

MICPRINT 

SMS 

1,629 

10-16-89 

12:  OOp 
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ROMAN 3 

SMS 

13,789 

10-16-89 

12 : OOp 

SBO 

SMS 

4,669 

10-16-89 

12 : OOp 

MATH 

BAT 

46 

07-23-92 

12 : 49p 

I  DAI 

BAT 

41 

07-25-92 

12 : 29p 

PCIDA 

EXE 

128,754 

04-19-93 

1 : 55p 

IDALINK 

BAT 

47 

04-15-93 

4  : 39p 

BOB 

FSP 

17,408 

07-22-92 

2  : 56p 

MATH 

C 

18,594 

04-19-93 

10  : 33a 

HELP 

IDA 

15,809 

03-24-93 

9:39a 

PCIDA 

TXT 

161 

10-05-92 

9 : 02a 

GRAFIX 

BAT 

44 

07-22-92 

4  : 12p 

IDAPRINT 

BAT 

24 

07-25-92 

1 : 03p 

READ 

ME 

340 

02-25-93 

2  : 47p 

DISK 

C 

21,600 

02-12-93 

3  : 56p 

GRAFIX 

C 

15,692 

02-12-93 

3  : 49p 

OUTPUT 

C 

24,085 

03-24-93 

9:39a 

IDACOMP 

BAT 

205 

04-19-93 

1 : 43p 

Total  files  listed: 
58  file(s) 


1,034,131  bytes 
407,040  bytes  free 


Kroutil  Disk  6,  Version  3  "C"  code 
MIDCOL. 


Volume  in  drive  A  has  no  label 
Volume  Serial  Number  is  411A-19D5 

Directory  of  A:\ 


MIDCOL 

EXE 

54,091 

02-17-93 

5  : 26p 

EXREPLAY 

DEF 

478 

05-01-93 

11:56a 

REPLAY 

BAT 

132 

05-01-93 

11:47a 

HEADERST 

DEF 

2,505 

05-04-92 

9  : 34p 

REPLAY 

C 

29,299 

05-01-93 

3  :  53p 

CONVINTF 

BAT 

137 

04-26-93 

8:53a 

MIDCOL 

BAT 

135 

08-15-92 

9  : 07p 

MIDDEF 

H 

2,816 

07-24-92 

7:13a 

READ 

ME 

463 

02-25-93 

2  : 45p 

CONVINTF 

C 

11,915 

05-03-93 

9:07a 

HEADERS 

DEF 

1,578 

09-27-92 

2 : 36p 

MIDDOC 

TXT 

21,479 

10-13-92 

2  : 38p 

MIDDOC 

WP 

21,600 

10-13-92 

2  : 38p 

F5 

FSP 

2,072 

02-13-93 

7  : 51p 

F6 

FSP 

2,072 

02-13-93 

7  : 51p 

F7 

FSP 

2,072 

02-13-93 

7  : 51p 

F8 

FSP 

2,072 

02-13-93 

7  : 52p 

SCALC 

DEF 

970 

04-02-93 

6  : 51p 

MIDCOL 

C 

68,614 

04-06-93 

9:02a 

CONVINTF 

EXE 

239,818 

04-27-93 

5  : 12p 

EXCONV 

DEF 

402 

04-27-93 

8:50a 

REPLAY 

EXE 

42,394 

05-03-93 

8:52a 

MENU 

H 

2,084 

01-18-90 

6 : 46p 

MIDCOLV 

BAT 

180 

05-04-93 

6:08a 

MIDCOLV 

EXE 

190,141 

05-07-93 

6:48a 

MENU 

C 

11,208 

01-18-90 

6 : 44p 

EXMIDCOL 

DEF 

359 

05-03-93 

8  : 22p 

MIDCOLV 

C 

71,963 

05-07-93 

6:46a 

SCCONV 

C 

14,719 

04-30-93 

9:22a 

SCCONV 

BAT 

99 

04-26-93 

8:54a 

SCCONV 

EXE 

175,742 

04-30-93 

9:22a 

Total  files  listed: 
31  file(s) 


973,609  bytes 
475,136  bytes  free 


Kroutil  Disk  7,  Version  3  PCIDA. 


Total  files  listed: 
27  file(s) 


581,210  bytes 
869,376  bytes  free 


Kroutil  Disk  8,  Version  3  MTRX  source. 


Volume  in  drive  A  has  no  label 
Volume  Serial  Number  is  0C42-19D4 


Directory  of  A:\ 


MTRX 

DEF 

318 

08-17-92 

12:08a 

MIDDEF 

H 

2,816 

07-24-92 

7:13a 

AMTRX 

EXE 

52,843 

04-15-93 

9:52a 

MTRXD 

BAT 

132 

08-17-92 

1 : 44p 

AMTRXD 

EXE 

49,403 

04-15-93 

9:53a 

MTRX 

BAT 

129 

08-17-92 

1 : 46p 

HEADERS 

DEF 

1,578 

09-27-92 

2  :36p 

FILTERl 

INC 

28,285 

04-02-93 

2  :  Olp 

MATGRAPH 

EXE 

45,034 

10-02-92 

3  : 52p 

MATGRAPH 

BAT 

112 

10-01-92 

10:45a 

MTRXD 

C 

16,243 

02-17-93 

10:19a 

DISCRIM1 

INC 

5,122 

04-03-93 

9:07a 

MTRX  2 

C 

42,997 

04-15-93 

9:56a 

MTRX 

EXE 

56,035 

04-15-93 

9:56a 

MTRX  2 

EXE 

65,443 

04-15-93 

9:57a 

MTRXD 

EXE 

52,595 

04-15-93 

9:55a 

MTRX2 

BAT 

132 

02-09-93 

3  : 37p 

MATGRAPH 

C 

5,516 

02-17-93 

10:18a 

FILTER2 

INC 

18,421 

04-02-93 

2  : 21p 

DISCRIM2 

INC 

5,122 

04-03-93 

9:06a 

MFILTER 

INC 

28,285 

04-02-93 

2:  Olp 

AFILTER1 

INC 

18,421 

04-02-93 

2 : 20p 

AFILTER2 

INC 

18,421 

04-02-93 

2  :  21p 

ADIS1 

INC 

5,122 

04-03-93 

9:04a 

ADIS2 

INC 

5,122 

04-03-93 

9:06a 

MDIS1 

INC 

5,122 

04-03-93 

9:07a 

MTRX 

C 

40,409 

04-15-93 

9:51a 

Total  files  listed: 
27  f ile {s ) 


569,178  bytes 
881,152  bytes  free 


ERDEC  Disk  1,  Edgewood  Chemical 
Database,  Disk  1  of  2. 


Volume  in  drive  A  has  no  label 


Directory  of  A:\ 


Volume  in  drive 

A  has  no  label 

Volume  Serial  Number  is 

382A-19D7 

Directory  of  A:\ 

DISK 

C 

22,870 

05-04-93 

10 : 08p 

RMATH 

C 

18,771 

07-23-92 

10 : 4  4p 

OUTPUT 

c 

24,067 

05-04-93 

7  : 42p 

HEADERS 

IDA 

1,504 

07-20-90 

9:58a 

PCIDA 

DOC 

3,061 

12-06-90 

10:02a 

INTF_TMP 

IDA 

40,960 

07-28-93 

9:22a 

SPEC_TMP 

IDA 

20,480 

07-28-93 

9:22a 

MICPRINT 

SMS 

1,629 

10-16-89 

12:  OOp 

ROMAN 3 

SMS 

13,789 

10-16-89 

12:  OOp 

SBO 

SMS 

4,669 

10-16-89 

12:  OOp 

MATH 

BAT 

46 

07-23-92 

12  : 49p 

IDA1 

BAT 

41 

07-25-92 

12  : 29p 

PCIDA 

EXE 

352,244 

07-27-93 

9  : 23p 

IDALINK 

BAT 

47 

04-15-93 

4  : 3  9p 

GRAFIX 

C 

15,908 

05-04-93 

10:44p 

PCIDA 

TXT 

161 

10-05-92 

9:02a 

GRAFIX 

BAT 

44 

07-22-92 

4  : 12p 

IDAPRINT 

BAT 

24 

07-25-92 

1 : 03p 

READ 

ME 

340 

02-25-93 

2  : 47p 

IDACOMP 

BAT 

205 

05-02-93 

8  : 51p 

EXTERNS 

IDA 

646 

05-04-93 

7  :38p 

HELP 

IDA 

15,809 

05-04-93 

7  : 58p 

MATH 

C 

18,556 

07-27-93 

9 : 22p 

DISK 

BAT 

34 

05-04-93 

5 : 38p 

DEFINES 

IDA 

376 

05-04-93 

10  : 47p 

TEMPS 

C 

14,697 

05-04-93 

9  : 05p 

PCIDA 

C 

10,232 

07-27-93 

8  : 49p 

27  file ( s ) 

581.210  bvtes 

IRSPEC 

EXE 

72,610 

02-02-88 

8:17a 

IRSPEC 

LIR 

610,100 

10-21-93 

9:35a 

IRSPEC01 

SPC 

8,301 

12-11-90 

1 : 06p 

IRSPEC02 

SPC 

8,225 

12-11-90 

1 : 08p 

IRSPEC03 

SPC 

8,234 

12-11-90 

1 : 12p 

IRSPEC05 

SPC 

8,226 

12-11-90 

2  :  21p 

IRSPEC04 

SPC 

8,235 

12-11-90 

2  : 23p 

IRSPEC06 

SPC 

8,279 

12-11-90 

2  :23p 

IRSPEC07 

SPC 

8,256 

12-11-90 

2  :24p 

IRSPEC08 

SPC 

8,481 

12-11-90 

2  :24p 

IRSPEC09 

SPC 

8,362 

12-11-90 

2  : 24p 

IRSPEC10 

SPC 

8,478 

12-11-90 

2  :25p 

IRSPEC11 

SPC 

8,533 

12-11-90 

2  : 25p 

IRSPEC12 

SPC 

8,556 

12-11-90 

2 : 26p 

IRSPEC13 

SPC 

8,318 

12-11-90 

2 : 26p 

IRSPEC14 

SPC 

9,232 

12-11-90 

2 : 27p 

IRSPEC15 

SPC 

8,475 

12-11-90 

2 : 28p 

IRSPEC16 

SPC 

9,647 

12-12-90 

9:34a 

IRSPEC 

BAK 

675 

12-12-90 

11:36a 

IRSPEC50 

ASP 

21,887 

12-14-90 

2 : 49p 

IRSPEC51 

ASP 

21,887 

12-14-90 

2 : 49p 

IRSPEC52 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC53 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC54 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC55 

ASP 

21,887 

12-14-90 

2  : 49p 

IRS PEC 5 6 

ASP 

21,887 

12-14-90 

2 : 49p 

IRSPEC57 

ASP 

21,890 

12-14-90 

2  :49p 

IRSPEC58 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC 5 9 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC60 

ASP 

21,887 

12-14-90 

2 : 49p 

IRSPEC 

CNT 

634 

02-14-91 

2  : 19p 

README! ! 

DOC 

1,008 

02-14-91 

2  :27p 

Total  files  listed: 
32  file (s) 


1,061,625  bytes 

388,096  bytes  free 


ERDEC  Disk  2,  Edgewood  Chemical 
Database,  Disk  2  of  2. 


Volume  in  drive 

A  has  no 

label 

Directory  of  A;\ 

IRSPEC49  ASP 

21,887 

12-14-90 

2 : 49p 

IRSPEC48  ASP 

21,887 

12-14-90 

2 : 49p 

IRSPEC47  ASP 

21,887 

12-14-90 

2  : 49p 
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IRSPEC20 

ASP 

21,956 

12-14-90 

2  : 47p 

IRSPEC19 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC18 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC17 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC16 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPECX5 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC14 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC13 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC12 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC11 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC10 

ASP 

21,887 

12-14-90 

2 : 47p 

IRS PEC 09 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC08 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC07 

ASP 

21,887 

12-14-90 

2 : 47p 

IRSPEC06 

ASP 

21,898 

12-14-90 

2:47p 

IRSPEC05 

ASP 

21,888 

12-14-90 

2 : 46p 

IRSPEC04 

ASP 

21,904 

12-14-90 

2 : 46p 

IRSPEC03 

ASP 

21,899 

12-14-90 

2 : 46p 

IRS PEC 02 

ASP 

21,887 

12-14-90 

2 : 46p 

IRS PEC 01 

ASP 

21,887 

12-14-90 

2 : 46p 

49  file(s) 

1,072,950  bytes 

Total  files  listed: 

49  file(s) 

1,072,950  bytes 

378,368  bytes  free 

IRSPEC46 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC45 

ASP 

21,887 

12-14-90 

2  : 49p 

IRSPEC44 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC43 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC42 

ASP 

21,899 

12-14-90 

2  : 48p 

IRSPEC41 

ASP 

21,925 

12-14-90 

2  : 48p 

IRSPEC40 

ASP 

21,903 

12-14-90 

2  : 48p 

IRSPEC3  9 

ASP 

21,890 

12-14-90 

2  : 48p 

IRSPEC3  8 

ASP 

21,911 

12-14-90 

2  : 48p 

IRSPEC37 

ASP 

21,896 

12-14-90 

2  : 48p 

IRSPEC36 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC35 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC34 

ASP 

21,887 

12-14-90 

2:48p 

IRSPEC33 

ASP 

21,887 

12-14-90 

2 : 48p 

IRSPEC32 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC3 1 

ASP 

21,889 

12-14-90 

2  : 48p 

IRSPEC30 

ASP 

22,156 

12-14-90 

2  : 48p 

IRSPEC29 

ASP 

21, 888 

12-14-90 

2  : 48p 

IRSPEC28 

ASP 

21, 887 

12-14-90 

2  : 48p 

IRSPEC27 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC26 

ASP 

21,889 

12-14-90 

2  : 48p 

IRSPEC25 

ASP 

21,887 

12-14-90 

2  : 48p 

IRSPEC24 

ASP 

21,888 

12-14-90 

2  : 47p 

IRSPEC23 

ASP 

21, 887 

12-14-90 

2  : 47p 

j IRSPEC22 

ASP 

21,887 

12-14-90 

2  : 47p 

IRSPEC21 

ASP 

21, 887 

12-14-90 

2  : 47p 
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Gas  Library  Databases  (GFI) 


Gas  Library  Databases  ( GFI) 

Two  gas  library  databases  were  obtained  (GFI)  for  use  with  HyGAS. 

•  Rauss_Sprouse  database 

•  ERDEC_Emery  database 


1.  Rauss_Sprouse  database 

The  Rauss_S prouse  database  (Table  1)  was  supplied  by  P.  Rauss  of  the  U.S.  Army 
Research  Laboratories  (ARL),  Ft.  Belvoir,  Virginia.  This  database  contains  80  spectra. 
This  includes  the  absorbance  spectra  for  78  different  gas  species  at  2  cm'1  spectral 
resolution.  (Nyquist  sampling  for  2  cm'1  spectral  resolution  requires  a  datapoint  spacing 
of  1  cm'1.)  One  gas  was  measured  at  two  concentrations  (ethylene  dibromide,  at  259.13 
ppm  and  12.4  ppm);  one  gas  was  measured  at  two  spectral  resolutions  (sulfur 
hexafluoride,  at  2  cm'1  and  0.5  cm'1). 

This  database  was  supplied  in  Spectra  Calc™/GRAMS™  (trademarks  of  Galactic 
Industries  Corporation)  binary  format,  which  we  converted  into  standard  JCAMP  format 
and/or  ASCII_XY  format,  for  incorporation  into  a  HIPS  spectral  library.  Each  of  the 
spectra  at  2  cm'1  resolution  contains  4150  datapoints  (the  0.5  cm'1  spectrum  contains 
14730  datapoints),  decreasing  in  x-value  from  4400  cm  1  to  400  cm  1  wavenumbers;  the 
y-value  is  in  absorbance  units.  The  original  measurements  for  this  gas  database  were 
done  by  Sprouse  Scientific  Systems,  to  specifications  defined  by  Rauss. 


Table  1.  HyGAS  GFI  Rauss_Sprouse  database.  The  80  spectra  included  in  the  Sprouse  database  are 
listed  numerically  by  filename  (spec  [nn]  [ r ] ).  Original  filenames  (i.e.  spectti.spc)  included  the  extension 
" .  spc"  (indicating  the  Spectra  Calc/GRAMS  format  of  Galactic  Industries  Corporation),  and  a  ”#” 
character  in  the  title,  which  was  not  recognized  as  a  valid  filename  by  GRAMS.  Descriptions  were 
extracted  from  the  converted  JCAMP  format  header.  Cross-reference  index  numbers  for  each  derived 
HyGAS  library  are  listed.  [Note  indexing  starts  at  0  for  HIPS  spectral  libraries,  and  1  for  the  HyGAS 
Neural  Net.]  Files  spec44a  and  spec44b  refer  to  two  different  species;  spec58  and  spec58r  are  the  same 
species,  as  are  spec7i  and  spec7ir.  [This  explains  the  original  index  numbers  ranging  from  1—77  rather 
than  the  expected  1-80.) 


Title  Description 


spec_01 

spec_02 

spec_03 

spec_04 

spec_05 

spec_06 

spec_07 

spec_08r 

spec_09 

spec_10r 

spec_llr 

spec_12 

spec_13 

spec_14 

spec_15 

spec_16 

spec_17 

spec_18r 

spec_19r 

spec_20 

spec_21 

spec_22r 

spec_23r 

spec_24r 

spec_25r 

spec__26r 

spec_27r 

spec_28r 

spec_29r 

spec_30r 


METHYL  CHLORIDE,  800  ppm 
ETHYL  CHLORIDE,  800  ppm 
ISOPROPYL  CHLORIDE,  49.54  ppm 
DIMETHYL  ETHER,  200  ppm 
DIMETHYL  DISULFIDE, 

PINACOLYL  ALCOHOL,  193.50  ppm 
2 -CHLOROETHANOL , 

FORMAMIDE,  404  ppm,  T=60  deg  C 

TRIMETHYL  PHOSPHITE,  19.00  ppm 

PYRIDINE,  144  ppm,  T=25  deg  C 

ALLYL  CHLORIDE,  352  ppm,  T=25  deg  C 

DIMETHYL  METHYLPHOSPHONATE,  46.52  ppm,  T=45  deg  C 

DIETHYL  METHYLPHOSPHONITE,  26.88  ppm,  T=45  deg  C 

DIETHYL  METHYLPHOSPHONATE,  40.24  ppm,  T=45  deg  C 

TRIISOPROPYL  PHOSPHITE,  14.68  ppm,  T=45  deg  C 

DIISOPROPYL  METHYLPHOSPHONATE,  22.96  ppm,  T=45  deg  C 

DIETHYL  METHYLPHOSPHONOTHIOATE, 

ACETONITRILE,  1696  ppm,  T=25  deg  C 
ACETONE,  146  ppm,  T=25  deg  C 
THIODIGLYCOL,  10,837.67  ppm,  T=250  deg  C 
DIISOPROPYL  AMINOETHANOL,  96.61  ppm,  T=73  deg  C 
BROMOBENZENE,  152  ppm,  T=60  deg  C 
ALLYLBENZENE,  284  ppm,  T=25  deg  C 
BROMOETHANE,  468  ppm,  T=25  deg  C 
NITROMETHANE,  99  ppm,  T=25  deg  C 
CYCLOHEXANONE,  87  ppm,  T=25  deg  C 
ACETIC  ANHYDRIDE,  49  ppm,  t=25  deg  C 
ACETIC  ACID,  129  ppm,  T=25  deg  C 
BENZYL  CHLORIDE,  242  ppm,  T=25  deg  C 
ETHYL  PHENYL  KETONE,  205  ppm,  T=60  deg  C 


Sprouse  Rauss-  HyGAS 

_80  Sprouse  Neural 

HIPS  _65  Net 

Library  HIPS  Index 

Index  Library 

Index 

0 

0 

1 

1 

1 

2 

2 

2 

3 

3 

3 

4 

4 

4 

5 

5 

5 

6 

6 

6 

7 

7 

7 

8 

8 

8 

9 

9 

9 

10 

10 

10 

11 

11 

11 

12 

12 

12 

13 

13 

14 

13 

14 

15 

14 

15 

16 

15 

16 

17 

16 

17 

18 

17 

18 

19 

18 

19 

20 

19 

20 

21 

20 

21 

22 

21 

22 

23 

22 

23 

24 

23 

24 

25 

24 

25 

26 

25 

26 

27 

26 

27 

28 

27 

28 

29 

28 

29 

Appendix  B-l 


HyGAS  Final  Report 


Gas  Library  Databases  (GFI) 


Title 

Description 

Sprouse 

_80 

HIPS 

Library 

Index 

Rauss- 

Sprouse 

__65 

HIPS 

Library 

Index 

HyGAS 

Neural 

Net 

Index 

spec  31 

THIONYL  CHLORIDE,  61.72  ppm 

30 

29 

30 

spec  32 

PHOSPHORUS  TRICHLORIDE,  19.26  ppm,  T=45  deg  C 

31 

30 

31 

spec  33 

PHOSPHOROUS  OXYCHLORIDE,  216.30  ppm,  T=45  deg  C 

32 

31 

32 

spec  34 

DIMETHYLPHOSPHORAMIDIC  DICHLORIDE,  70.60  ppm,  T=73  deg  C= 

33 

spec  35 

DIMETHYLPHOSPHITE,  54.97  ppm,  T=45  deg  C 

34 

spec  36 

DIISOPROPYL  HYDROGEN  PHOSPHITE,  20.22  ppm,  T=45  deg  C 

35 

32 

33 

spec  37r 

TETRAHYDROFURAN,  133  ppm,  T=25  deg  C 

36 

33 

34 

spec  38r 

ETHYL  ACETATE,  45  ppm,  T=25  deg  C 

37 

34 

35 

spec  39r 

METHACRYLATE,  80  ppm,  T=25  deg  C 

38 

35 

36 

spec  40 

METHYL PHOSPHONYL  DICHLORIDE,  36.37  ppm,  T=45  deg  C 

39 

36 

37 

spec  41r 

FORMIC  ACID,  119  ppm,  T=25  deg  C 

40 

37 

38 

spec  42 

METHYL  DICHLOROPHOSPHINE,  103.07  ppm,  T=45  deg  C 

41 

38 

39 

spec  43 

HYDROGEN  SULFIDE,  20,000  ppm 

42 

39 

40 

spec  44a 

SULFUR  DICHLORIDE,  53.46  ppm,  T=45  deg  C 

43 

spec  44b 

SULFUR  MONOCHLORIDE,  63.02  ppm,  T=45  deg  C 

44 

spec  45 

ETHYLENE,  55.7  ppm 

45 

40 

41 

spec  46 

METHANOL,  52.8  ppm 

46 

41 

42 

spec  47 

ETHANOL,  69.1  ppm 

47 

42 

43 

spec  48 

ISOPROPANOL,  34.9  ppm 

48 

43 

44 

spec  49 

VINYL  CHLORIDE,  49.20  ppm 

49 

44 

45 

spec  50 

ETHYLENE  OXIDE,  55.8  ppm 

50 

45 

46 

spec_51 

DIMETHYLAMINE,  400  ppm 

51 

46 

47 

spec_52 

TRIETHYLAMINE,  38.59  ppm 

52 

47 

48 

spec_53 

CARBON  TETRACHLORIDE,  11.1  ppm 

53 

spec_54 

CHLOROFORM,  20.1  ppm 

54 

48 

49 

spec_55 

HYDROGEN  CHLORIDE,  lOOppm 

55 

spec_56 

PARATHION,  4,293.97  ppm,  T=200  deg  C 

56 

49 

50 

spec_57 

DIETHYL  ETHER,  172.59  ppm 

57 

50 

51 

spec_58 

SULFUR  HEXAFLUORIDE,  5.01  ppm 

58 

51 

52 

spec_58r 

SULFUR  HEXAFLUORIDE,  5.01  ppm,  T=25  deg  C,  0.5  cm-1  Resolution 

59 

spec_5  9 

METHANE,  50.30  ppm 

60 

52 

53 

spec_60 

ETHANE,  48.50  ppm 

61 

53 

54 

spec_61 

PROPANE,  54.50  ppm 

62 

spec_62 

PROPYLENE,  56.20  ppm 

63 

54 

55 

spec_63 

1,3 -BUTADIENE,  200  ppm 

64 

55 

56 

spec_64 

n-BUTANE,  52.5  ppm 

65 

56 

57 

spec_65 

C i s - 2 - PENTENE ,  83.09  ppm 

66 

57 

58 

spec_66 

ETHYL  MERCAPTAN,  363.21  ppm 

67 

58 

59 

spec_67 

BUTYL  MERCAPTAN,  133.92  ppm 

68 

59 

60 

spec_68 

BENZENE,  20.2  ppm 

69 

60 

61 

spec_69 

SULFUR  DIOXIDE,  49.8  ppm 

70 

61 

62 

spec_70 

NITROUS  OXIDE,  47.2  ppm 

71 

62 

63 

spec_71 

ETHYLENE  DIBROMIDE,  259.13 

72 

63 

64 

spec_71r 

ETHYLENE  DIBROMIDE,  12.4  ppm,  T=25  deg  C 

73 

spec_72r 

NITROETHANE,  62  ppm,  T=25  deg  C 

74 

64 

65 

spec_73 

CARBON  DIOXIDE,  40.4  ppm 

75 

spec_74 

CARBON  MONOXIDE,  50.40  ppm 

76 

spec_75 

WATER,  1,492.30  ppm,  T=45  deg  C 

77 

spec_7  6 

n-HEXANE,  89  ppm,  T=25  deg  C 

78 

spec_77 

n-HEPTANE,  79  ppm,  T=25  deg  C 

79 

2.  ERDECJEmery  database 

The  ERDECJEmery  database  (Table  2)  was  supplied  by  R.  Kroutil  of  the  U.S.  Army 
Edgewood  Research  Development  and  Engineering  Center  (ERDEC),  Edgewood, 
Maryland.  This  database  contains  60  spectra.  This  includes  absorbance  spectra  for 
selected  chemical  warfare  (CW)  agents,  CW  agent  simulators,  and  potential  contaminants 
such  as  engine  vapors,  smoke  and  dust,  all  at  2  cm'1  spectral  resolution. 

This  database  was  supplied  in  a  generic  ASCII  format,  which  we  converted  into 
ASCII_XY  format  for  incorporation  into  a  HIPS  spectral  library.  Each  of  the  original 
spectra  contains  2184  datapoints,  increasing  in  x-value  from  388  cm'1  to  2500  cm'1 
wavenumbers;  y-values  are  in  optical  cross-section  units.  Not  all  spectra  include  valid 
datavalues  over  the  full  spectral  range.  Ranges  of  datapoints  with  y-values  of  zero  were 
omitted,  defining  the  valid  subset  spectral  range  for  each  spectrum;  minimum  and 
maximum  wavenumbers  for  each  spectrum  are  indicated  in  the  appropriate  columns  of 
Table  2.  The  original  measurements  for  this  database  were  [presumably]  done  by 
ERDEC,  to  specifications  defined  by  the  Army,  and  were  compiled  into  this  database  by 
Silvio  Emery,  at  ERDEC. 
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Table  2.  HyGAS  GFI  Emery  database.  The  60  spectra  included  in  the  Emery  database  are  listed  by 
original  filename  (irspectnnj).  Descriptions  were  extracted  from  the  original  ASCII  header.  Index 
numbers  for  the  derived  HyGAS  library  also  are  included,  for  cross-reference  by  library  index  number. 
Minimum  and  maximum  wavenumber  values  are  included  for  each  spectrum,  indicating  the  effective  range 
of  valid  datavalues. 


Title 

Description 

HIPS 

Library 

Index 

min 

wavenumber 

cm"1 

max 

wavenumber 

cm’1 

irspecOl 

EA0001  GA  (LIQUID) 

0 

396 

2500 

irspec02 

EA0002_GA_ (VAPOR) 

1 

427 

2500 

irspec03 

EA0003_GB_ (LIQUID) 

2 

404 

2500 

irspec04 

EA0004_GB_ (VAPOR) 

3 

411 

2500 

irspec05 

EA0005_GD_( LIQUID) 

4 

419 

2500 

irspec06 

EA0006  GD_( VAPOR) 

5 

427 

2500 

irspec07 

EA000 7_VX_ ( LIQUID ) 

6 

411 

2500 

irspec08 

EA0008_CS_ (LIQUID) 

7 

466 

2500 

irspec09 

EA0009  CN  (LIQUID) 

8 

473 

2500 

irspeclO 

EA0010  HD_{ LIQUID) 

9 

481 

2500 

irspecll 

EA0011_HQ_( LIQUID) 

10 

489 

2500 

irspecl2 

EA00 1 2_HD_ ( VAPOR ) 

11 

496 

2500 

irspeclB 

EA0013_DM_( LIQUID) 

12 

396 

2500 

irspecl4 

EA0014  HC_( VAPOR) 

13 

396 

2500 

irspecl5 

EA0015  HC_( PARTICULATES) 

14 

404 

2500 

irspecl6 

EA0016  GREEN  SMOKE_ (VAPOR ) 

15 

411 

2500 

irspecl7 

EA0017  GREEN  SMOKE_ ( PARTICULATES ) 

16 

419 

2500 

irspecl8 

EA0018  RED_SMOKE_( VAPOR) 

17 

427 

2500 

irspecl9 

EA0019  RED  SMOKE_ ( PARTICULATES) 

18 

435 

2500 

irspec20 

EA0020  YELLOW_SMOKE_ (VAPOR) 

19 

401 

2500 

irspec21 

EA0021  YELLOW_SMOKE_ ( PARTICULATES) 

20 

450 

2500 

irspec22 

EA0022  VIOLET  SMOKE_ (VAPOR ) 

21 

458 

2500 

irspec23 

EA0023  VIOLET  SMOKE_( PARTICULATES) 

22 

466 

2500 

irspec24 

EAO  02  4_GRASS_SMOKE_ ( VAPOR ) 

23 

473 

2500 

irspec25 

EA0025  GRASS_SMOKE_( PARTICULATES) 

24 

481 

2500 

irspec26 

EAO  02  6_COTTON_SMOKE_ ( VAPOR ) 

25 

489 

2500 

irspec27 

EA0027  COTTON_SMOKE_( PARTICULATES) 

26 

496 

2500 

irspec28 

EAO  0  2  8„RUBBER_SMOKE_ ( VAPOR ) 

27 

504 

2500 

irspec29 

EA0029  WOOL„SMOKE_( VAPOR) 

28 

512 

2500 

irspec30 

EA003 0_WOOL_SMOKE_ ( PARTICULATES ) 

29 

388 

2500 

irspec31 

EAO  03 1_PAPERS_ ( VAPOR ) 

30 

527 

2500 

irspec32 

EAO  0  3  2_GASOLINE_ENG INE_EXHAUST 

31 

404 

2500 

irspec33 

EAO  0  3  3 _GASOL INE_ ( VAPOR ) 

32 

411 

2500 

irspec34 

EAO  0  3  4_ORANGE_SMOKE 

33 

419 

2500 

irspec35 

EAO  0  3  5_WHITE_SMOKE 

34 

427 

2500 

irspec36 

EA003 6 _BLUE_SMOKE 

35 

435 

2500 

irspec37 

EAO 03 7_METHYLCELL 

36 

442 

2500 

irspec38 

EA0038J3ICDI 

37 

504 

2500 

irspec39 

EAO  0  3  9_TH  IOLAMINE 

38 

512 

2500 

irspec40 

EA0040_DIMP_ (LIQUID) 

39 

700 

1350 

irspec41 

EAO  04 1_DMMP_ ( VAPOR ) 

40 

701 

1350 

irspec42 

EA0042_DMMP_ (LIQUID) 

41 

700 

1350 

irspec43 

EA0043_ETHANOL_ (VAPOR) 

42 

726 

1350 

irspec44 

EAO 04  4_FT . _BENNING_OUST 

43 

435 

1608 

irspec45 

EA0045_FT . _BENNINGJDUST 

44 

450 

1608 

irspec46 

EAO  0  4  6 _  S A  LT_LAK  E_C I TY_DUST 

45 

427 

1613 

irspec47 

EA 0  0  4  7_PANAMA__DUST 

46 

435 

1613 

irspec48 

EAO  04  8_S INGLEEDON_DUST 

47 

401 

1605 

irspec49 

EAO  0  4  9_SOUTH_MIAMI_DUST 

48 

466 

1613 

irspec50 

EAO  0  5  0_PANAMA_DUST 

49 

401 

1605 

irspec51 

EAO  0  5 l_SINGLEEDON JDUST 

50 

401 

1605 

irspec52 

EAO 0 52_FT . _MCCLELLEN_DUST 

51 

401 

1600 

irspec53 

EAO 0 5 3 _FT . _BENNING_DUST 

52 

401 

1601 

irspec54 

EA0054_FT . _BENNING_DUST 

53 

401 

1599 

irspec55 

EAO  0  5  5_WP_ ( VAPOR ) 

54 

741 

1428 

irspec56 

EAO  0  5  6_NEUTRAL 

55 

388 

2500 

irspec57 

EAO 057_UN KNOWN 

56 

747 

1303 

irspec58 

EAO  0  5  8_UNKNQWN 

57 

708 

1400 

irspec59 

EA0059_UNKNOWN 

58 

666 

1400 

irspec60 

EAO  0  6  0_UNKNOWN 

59 

699 

1400 

3.  HyGAS  Spectral  Libraries 

For  use  with  HyGAS,  three  HIPS  spectral  libraries  were  created  from  these  databases. 

•  HyGAS_Sprouse_8 0 

•  HyGAS_Rauss-Sprouse_65 

•  HyGAS_ERDEC_Emery 

The  first  two  of  these  libraries  are  derived  from  the  Sprouse  database;  the  third  library  is 
derived  from  the  ERDEC  Emery  database.  Descriptions  for  these  libraries  are  listed  in 
Table  3. 
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Table  3.  HyGAS  HIPS  Spectral  Library  Descriptions.  This  information  is  extracted  from  the  HIPS 
ASCII  file  $dod_dir/library/ tables /General .  lib.  Included  for  each  library  is  a  list  of  HIPS  spectral 
sensors  to  which  the  library  has  been  resampled. 


Name 

HyGAS__Sprouse_8  0 

Description 

"Gas  Spectra  <2cm-l  resolution)  measured  by  Sprouse  Scientific 
Systems  for  P.  Rauss  of  US  Army  NVEOD,  4150  pts . ,  4400-400  cm- 
1,  DECREASING  wavenumber;  (also  includes  SF6  at  0.5cm-l 
resolution) " 

Directory 

/usr/HIPS  /  LIBRARIES  /  HyGAS_Sprouse__8  0 

Mode 

420 

Sensors 

Rauss_449_spec . sensor  Kroutil_513_spec . sens 

Kroutil__113_spec .  sens 

Name 

HyGAS_Rauss-Sprouse_65 

Description 

"Gas  Spectra  (2cm-l  resolution)  measured  by  Sprouse  Scientific 
Systems  for  P.  Rauss  of  US  Army  NVEOD,  selected  65  items,  449 
pts.,  833-1265  cm-1,  sorted  by  increasing  wavenumber, 

NORMALIZED  0  to  1" 

Directory 

/usr /HI PS /LIBRARIES /HyGAS_Rauss-Sprouse_6 5 

Mode 

420 

Sensors 

Rauss_449_spec. sensor 

Name 

:  HyG  AS_ERDEC_JEme  r  y 

Description : 

:  "Wavenumber  vs.  optical  cross  section  of  chemicals  stored  in 

the  IRSPEC  database  V.1.0  received  from  Edgewood  Arsenal  and 
compiled  by  Dr.  Silvio  L.  Emery" 

Directory 

/ us r / HI P S / L IBRARI ES / HyGAS_ERDEC_Emery 

Mode 

420 

Sensors 

Rauss_449_spec . sensor  Kroutil_113_spec . sens 

The  HIPS  spectral  library  HyGAS_Sprouse_80  includes  all  80  spectra  from  the  Sprouse 
database,  at  original  spectral  resolution  and  at  original  absorbance  units.  Spectra  included 
in  this  library  are  referenced  by  HIPS  library  index  (0-79)  in  Table  1. 

The  HIPS  spectral  library  HyGAS_Rauss-sprouse_65  includes  65  selected  spectra  from 
the  full  Sprouse  database.  For  this  library,  data  values  outside  the  spectral  region  of 
interest  (-8-12  pm,  or  -800-1250  cm'1)  were  eliminated.  Gas  species  with  no  significant 
spectral  features  in  this  region  also  were  eliminated,  as  were  any  replicated  species.  The 
y-axis  values  for  the  remaining  65  spectra  were  then  scaled  (minimum  to  maximum 
within  each  spectrum)  to  lie  within  the  arbitrary  range  0  to  1.  [This  is  a  standard 
technique  for  creating  spectral  search  library  datasets,  and  is  the  method  used  by  Rauss 
for  his  neural  net  training  library.  This  method  also  maximizes  the  variability  in 
absorbance  as  a  function  of  wavelength,  which  is  the  primary  identifying  metric  for  any 
gas  species  (i.e.  its  unique  and  characteristic  "fingerprint"  spectral  signature).]  Spectra 
included  in  this  library  are  referenced  by  HIPS  library  index  (0-64)  in  Table  1.  These  are 
the  same  spectra  included  in  the  HyGAS  Neural  Net  search  library,  also  referenced  in 
Table  1.  [Note  the  Neural  Net  index  values  are  1-65,  while  the  corresponding  HIPS 
library  index  values  are  0-64.] 

The  HIPS  library  HyGAS_ERDEC_Emery  includes  all  60  spectra  in  the  original  Emery 
database,  but  only  within  the  defined  valid  spectral  range.  This  range  is  indicated  in 
Table  2.  Although  this  range  varies,  all  spectra  include  the  primary  spectral  region  of 
interest  (-8-12  pm,  or  -800-1250  cm'1)  for  HyGAS. 
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4.  HIPS  Sensor  Files 

HIPS  can  associate  a  spectral  sensor  characteristics  (or  instrument)  file  with  any  HIPS 
dataset  (Table  4).  Such  a  file  will  define  a  spectral  value  for  each  spectral  channel  (or 
index)  in  a  HIPS  dataset.  This  is  useful  for  spectral  plots,  when  the  HIPS  user  needs  to 
assign  a  physical  meaning  to  the  spectral  axis,  particularly  so  in  instances  where  the 
functional  relationship  between  the  spectral  channels  is  not  a  simple  linear  function. 

HIPS  also  uses  the  spectral  sensor  characteristics,  albeit  not  directly,  when  applying  a 
library  search  algorithm.  Mathematically,  HIPS  compares  the  unknown  spectrum  to  the 
library  spectrum,  channel  by  channel.  In  practice,  this  means  the  spectral  channels  in  a 
spectral  library  should  correlate  directly  with  the  spectral  channels  in  the  unknown.  If  a 
spectral  library  is  created  with  the  same  spectral  sensor  characteristics  as  the  unknown, 
HIPS  can  simply  search  against  the  so-called  Raw  spectrum  in  the  library.  If,  however, 
the  samples  in  a  spectral  library  are  not  consistent  (such  as  is  the  case  with  the  Emery 
library),  or  there  is  not  a  1:1  correlation  between  the  spectral  channels  of  the  unknown 
and  the  spectral  channels  in  the  library,  the  library  must  first  be  "convolved,"  or 
resampled  to  match  the  spectral  sensor  characteristics  of  the  unknown  species.  In  this 
case,  HIPS  creates  an  additional  binary  version(s)  of  the  library,  with  data  values 
resampled  to  match  the  spectral  characteristics  of  the  specified  sensor(s). 

Table  4.  Associating  a  HIPS  spectral  sensor  characteristics  file.  This  is  an  example  abs  file  in  a  HIPS 
dataset  directory,  which  defines  the  parameters  of  the  associated  HIPS  data  file  data.xys.  The  instrument 
field  indicates  an  associated  spectral  sensor  characteristics  file  Rauss_449_spec.  sensor  which  defines  the 
spectral  axis  values  for  each  of  the  449  spectral  channels.  This  sensor  must  exist  in  the  directory 
$ dod_dir/library/ sensors,  for  HIPS  plotting  functions  to  assign  values  to  the  independent  variable,  or 
spectral  axis,  in  a  spectral  plot. 


OWNER  =  "hardy:" 

DATE  =  "Thu 

Sep  8 

11:54:58 

1994" 

DATATYPE  =  float 

X_AXIS_SIZE 

=  240 

Y_AXIS_SIZE 

=  150 

S_AXIS_SIZE 

=  449 

INSTRUMENT  = 

Rauss_ 

_449_spec . 

.  sensor 

The  sensors  supplied  for  use  with  HyGAS  are  listed  in  Table  5.  These  sensors  were 
defined  to  match  specific  HyGAS  test  datasets,  either  for  use  with  HyGAS  plotting 
functions,  or  for  use  with  the  HyGAS  search  libraries.  For  HyGAS,  each  of  the  spectral 
libraries  has  been  resampled  to  match  the  spectral  sensor  characteristics  for  appropriate 
simulated  HyGAS  test  datasets.  Table  3  indicates  the  sensors  for  which  each  of  these 
libraries  has  been  resampled  to  match. 

The  sensor  Kroutii_i024_intf .  sen  was  defined  to  match  the  characteristics  of  the  field 
data  supplied  by  Kroutil  as  GFI  for  use  with  HyGAS.  This  data  was  supplied  as  1024 
data  point  interferograms,  at  unit  channel  spacing,  ranging  from  channel  index  0  to 
channel  index  1023.  This  sensor  was  defined  for  use  during  plotting  only,  and  is  not 
applicable  for  use  with  spectral  libraries  (no  interferogram  domain  spectral  libraries  were 
attempted). 

The  sensor  Kroutil _ 513 _ spec .  sens  was  defined  to  match  the  recommended  Spectra 

Calc  functional  translation  of  the  Kroutil_MIDAC  interferograms.  Spectra  Calc 
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generates  a  spectral  domain  of  513  data  points  from  the  raw  interferograms  of  1024  data 
points.  This  sensor  was  defined  both  for  plotting  and  for  preliminary  testing  with  a 
spectral  library.  Only  the  Sprouse_80  library  was  resampled  to  this  sensor.  With  a 
spectral  range  of  0  to  1974.75  cm'1  wavenumbers  (at  8  cm'1  spectral  resolution,  Nyquist 
sampling  requires  a  datapoint  spacing  of  4  cm1),  much  of  the  spectral  data  for  this  sensor 
falls  outside  the  range  of  valid  library  spectral  data  (typically  400  cm'1  to  4400  cm'1 
wavenumbers);  in  fact,  in  order  to  function  with  this  sensor,  such  a  spectral  library  must 
be  extrapolated,  essentially  "padded"  with  invalid  zero-data-value  filler  points,  on  one 
end  of  its  valid  spectral  range,  to  match  the  spectral  characteristics  of  the  sensor.  This 
approach  is  not  recommended  in  practice,  and  was  applied  for  testing  purposes  only. 

The  sensor  Kroutii_5i2_spec .  sens  was  defined  to  match  an  in-house  developed  FFT 
method,  which  generated  an  identical  spectral  domain  to  that  for  the  Kroutil  FFT  code, 
but  with  one  less  data  point.  The  sensor  Kroutii_5i2 .  spec  is  functionally  identical  to 
this  sensor,  and  is  included  purely  for  backwards  compatibility  with  an  outdated  version 
of  the  HyGAS  code. 

The  sensor  Kroutii_ii3_spec .  sens  is  a  subset  of  the  other  Kroutil  spectral  domain 
sensors,  which  eliminates  those  datapoints  outside  the  primary  spectral  region  of  interest 
(-8-12  pm,  or  -800-1250  cm'1)  to  HyGAS.  Restricting  this  sensor  to  the  primary 
spectral  region  of  interest  avoids  any  difficulties  associated  with  trying  to  extrapolate 
beyond  the  range  of  valid  spectral  data  for  the  library.  The  Emery  library  and  the 
Sprouse_80  library  are  both  resampled  to  this  sensor,  at  8  cm'1  resolution.  (Since  the 
Rauss-Sprouse_65  library  already  is  a  spectral  subset  of  the  Sprouse_80  library,  little  if 
anything  is  to  be  gained  by  duplicating  this  sensor  for  the  subset  Rauss-Sprouse_65 
library.  The  sole  difference  would  be  in  the  values  associated  with  the  re-scaled  y-axis, 
and  this  difference  can  be  tested  using  the  Rauss_449  sensor,  described  below.) 

The  sensor  Rauss_449_spec .  sensor  is  defined  to  match  the  spectral  characteristics  of 
the  field  data  (from  the  Rauss-Howden  spectrometer)  supplied  by  Rauss  as  GFI  for  use 
with  HyGAS.  These  data  were  supplied  as  449  data  points  ranging  from  1265  cm'1  to  833 
cm"1  wavenumbers,  at  2  cm'1  spectral  resolution  (Nyquist  sampling  at  1  cm'1  datapoint 
spacing).  All  three  HyGAS  spectral  libraries  were  resampled  to  this  sensor  resolution. 

A  detailed  listing  of  the  sensor  files  defined  for  HyGAS  is  given  in  Table  6.  This  listing 
details  the  format  of  a  HIPS  sensor  file;  however,  only  the  first  five  and  last  five  channels 
of  each  sensor  are  listed. 

Table  5.  Sensors  supplied  for  use  with  HyGAS.  This  is  a  listing  of  the  HIPS  directory 
$dod_dir/ library/ sensors.  Sensor  names  include  a  reference  to  a  HyGAS  simulation  cube  test  dataset 
(Kroutil  or  Rauss),  the  number  of  spectral  channels,  and  a  reference  to  whether  the  sensor  is  an 
interferogram  domain  sensor,  or  a  spectral  domain  sensor.  The  sensor  Kroutii_5i2 .  spec .  tbi  is 
functionally  identical  to  the  sensor  Kroutii_5i2_spec .  sens .  tbi,  and  is  included  for  backwards 
compatibility  with  an  outdated  version  of  the  HyGAS  code.  Actual  sensor  names  omit  the  suffix  .  tbi. 

Kroutil_1024__intf .  sen .  tbi 
Kroutil_513_spec . sens . tbi 
Kroutil_512_spec . sens . tbi 
[Kroutil_512.spec .tbi] 

Kroutil_113_spec . sens . tbi 
Rauss_449_spec . sensor . tbi _ 
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Table  6.  Detailed  format  listing  of  HIPS  sensor  files  for  HyGAS.  Only  the  first  five  and  last  five 
channels  of  each  sensor  are  included. 


==>  Kroutil_1024_intf . sen. tbl  <== 
Name  :  K 

Description  :  " 

Channels  :  1 

Band  Table  : 

0.000000  0.100000  1.000000  0.00 
1.000000  0.100000  1.000000  0.00 


:  Kroutil_1024_intf . sen 

:  "Kroutil  MIDAC  1024pt  inter ferogram,  0-1023" 
:  1024 

0.000000 

0.000000 


2.000000  0.100000  1.000000  0.000000 
3.000000  0.100000  1.000000  0.000000 

4.000000  0.100000  1.000000  0.000000 


1019.000000  0.100000  1.000000  0.000000 

1020.000000  0.100000  1.000000  0.000000 

1021.000000  0.100000  1.000000  0.000000 

1022.000000  0.100000  1.000000  0.000000 

1023.000000  0.100000  1.000000  0.000000 


==>  Kroutil_113_spec .sens .tbl  <= 

Name 

Description  : 

spectra) ■ 

Channels  : 

Band  Table  : 

833.099976  3.860000  1.000000  0.000000 

836.950012  3.860000  1.000000  0.000000 

840.809998  3.860000  1.000000  0.000000 

844.669983  3.860000  1.000000  0.000000 

848.530029  3.860000  1.000000 


Kroutil_113_spec . sens 

"Kroutil  spectral  sensor,  113pts,  B33-1265cm-l  {channels  217-329  in  512  or  513pts 


3.860000 

3.860000 


3.860000  1.000000 


0.000000 

0.000000 

0.000000 


1249.650024  3.860000  1.000000  0.000000 
1253.500000  3.860000  1.000000  0.000000 
1257.359985  3.860000  1.000000  0.000000 
1261.219971  3.860000  1.000000  0.000000 
1265.069946  3.860000  1.000000  0.000000 


==>  Kroutil_512_spec . sens . tbl  <== 

Name  :  Kroutil_512_spec . sens 

Description  :  "simulated  spectral  sensor  for  Kroutil /MIDAC  FTIR  spectrometer,  for  use  with  HyGAS  FT 

function  that  generates  512  datapoints  from  a  1024  point  inter ferogram" 

Channels  :  512 

Band  Table  : 

0.000000  3.860000  1.000000  0.000000 


3.860000  3.860000  1.000000 

7.710000  3.860000  1.000000 


0.000000 

0.000000 


11.570000  3.860000  1.000000  0.000000 

15.430000  3.860000  1.000000  0.000000 


1955.469971  3.860000  1.000000  0.000000 
1959.319946  3.860000  1.000000  0.000000 
1963.180054  3.860000  1.000000  0.000000 
1967.040039  3.860000  1.000000  0.000000 

1970.890015  3.860000  1.000000  0.000000 


==>  Kroutil_513_spec . sens . tbl  <== 

Name  :  Kroutil_5 

Description  :  "Kroutil  i 

Channels  :  513 

Band  Table  : 

0.000000  3.860000  1.000000  0.000000 

3.860000  3.860000  1.000000  0.000000 

7.710000  3.860000  1.000000  0.000000 

11.570000  3.860000  1.000000  0.000000 

15.430000  3.860000  1.000000  0.000000 

1959.319946  3.860000  1.000000  0.000000 

1963.180054  3.860000  1.000000  0.000000 

1967.040039  3.860000  1.000000  0.000000 

1970.890015  3.860000  1.000000  0.000000 

1974.750000  3.860000  1.000000  0.000000 


Kroutil_513_spec . sens 

"Kroutil  spectral  sensor,  513pts,  0-1974 . 75cm- 1 
513 


==>  Rauss_449_spec . sensor . tbl  <=- 

Name  :  Rauss_449_spec . sensor 

Description  :  "simulated  Howden  field  sensor,  449  spectral  channels,  833-1265  inverse  cm,  from  data 

supplied  by  Rauss  as  GFI,  bandwidth  defined  by  datapoint  spacing" 

Channels  :  449 

Band  Table  : 

833.000000  0.964000  1.000000  0.000000 

833.960022  0.964000  1.000000  0.000000 

834.929993  0.964000  1.000000  0.000000 

835.890015  0.964000  1.000000  0.000000 

836.859985  0.964000  1.000000  0.000000 

1261.130005  0.964000  1.000000  0.000000 

1262.099976  0.964000  1.000000  0.000000 

1263.060059  0.964000  1.000000  0.000000 

1264.020020  0.964000  1.000000  0.000000 

1265.000000  0.964000  1.000000  0.000000 
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5.  Additional  notes 

The  HIPS  spectral  libraries  supplied  with  HyGAS  can  be  plotted,  listed  and  accessed 
using  the  standard  HIPS  functions. 

Sample  HIPS  command  lines  used  to  create  these  libraries  and  add  the  spectra  to  them  are 
given  in  Tables  7  and  8.  These  command  line  files  were  supplied  with  HyGAS.  The 
respective  files  are  in  directory  paths 

$hygas/Rauss/Sprouse_data/ascii_xy/xstuf f 

and 

$hygas / Emery /ascii_xy/xs tuff. 

Subdirectory  references  are  indicated  relative  to  these  path  names. 

Table  7.  Sample  HIPS  command  lines  for  creating  HIPS  spectral  libraries  for  HyGAS.  These 
command  lines  can  be  executed  as  UNIX  shell  scripts  (ex.  sh  do_criib_sprouse_80)  from  the  appropriate 
HyGAS  directory. 


==>  do_crlib_Sprouse_80  <== 

crlib  -t  "Gas  Spectra  (2cm-l  resolution)  measured  by  Sprouse  Scientific  Systems  for  P.  Rauss  of  US  Army  NVEOD,  4150 
pts.,  4400-400  cm-1,  DECREASING  wavenumber;  (also  includes  SF6  at  0.5cm-l  resolution)'  -S  -G  -d  /usr/HIPS/LIBRARIES 
HyGAS_Spr ous e_ 8 0 


==>  do_crlib_Sprouse_65  <== 

crlib  -t  'Gas  Spectra  (2cm-l  resolution)  measured  by  Sprouse  Scientific  Systems  for  P.  Rauss  of  US  Army  NVEOD,  selected 
65  items,  449  pts.,  833-1265  cm-1,  sorted  by  increasing  wavenumber,  NORMALIZED  0  to  1*  -S  -G  -d  /usr/HIPS/LIBRARIES 
HyGAS_Raus s - Sprouse_6  5 


==>  do_crlib [_ERDEC_Emery ]  <== 

crlib  -t  "Wavenumber  vs.  optical  cross  section  of  chemicals  stored  in  the  IRSPEC  database  V.1.0  received  from  Edgewood 
Arsenal  and  compiled  by  Dr.  Silvio  L.  Emery"  -S  -G  -d  /usr/HIPS/LIBRARIES  HyG AS_ERDEC_Eme ry 


Table  8.  Sample  HIPS  command  lines  for  adding  spectra  to  HIPS  spectral  libraries  for  HyGAS. 
Only  the  first  ten  lines  of  each  sample  file  are  listed.  These  command  lines  can  be  executed  as  UNIX  shell 
scripts  (ex.  sh  do_add2iib_sprouse_80)  from  the  appropriate  HyGAS  directory. 


==>  do_add21ib_Sprouse_80  <== 


add_spectrum  -f 
add_spec tr um  -f 
add_spectrum  -f 
add_spectrum  -f 
add_spectrum  -f 
add_spec  t r um  -f 
add_spectrum  -f 
add_spectrum  -f 
add_spectrum  -f 
add_spectrum  -f 


2  -t  "METHYL  CHLORIDE,  800  ppm'  -L  HyGAS_Sprouse_80  -n  spec_01  . . /unsorted_4150/spec_01 .prn 
2  -t  "ETHYL  CHLORIDE,  800  ppm'  -L  HyGAS_Sprouse_80  ~n  spec_02  . . /unsorted_4150/spec_02 .prn 
2  -t  "ISOPROPYL  CHLORIDE,  49.54  ppm"  -L  HyGAS_Sprouse_80  -n  spec„03  . . /unsorted_4150/spec_03 .prn 
2  -t  "DIMETHYL  ETHER,  200  ppm"  -L  HyGAS_Sprouse_80  -n  spec_04  . . /unsorted_4150/spec_04 .prn 
2  -t  "DIMETHYL  DISULFIDE,"  -L  HyGAS_Sprouse_80  -n  spec_05  . . /unsorted_4150/spec_05 .prn 
2  -t  " PINACOLYL  ALCOHOL,  193.50  ppm"  -L  HyGAS_Sprouse_80  -n  spec_06  . . /unsorted_4150/spec_06 .prn 
2  -t  "2-CHLOROETHANOL, "  -L  HyGAS_Sprouse_80  -n  spec_07  . . /unsorted_4150/spec_07 .prn 
2  -t  "FORMAMIDE,  404  ppm,  T=60  deg  C"  -L  HyGAS_Sprouse_80  -n  spec_08r  . . /unsorted_4150/spec_08r .prn 
2  -t  "TRIMETHYL  PHOSPHITE,  19.00  ppm"  -L  HyGAS_Sprouse_80  -n  spec_09  . . /unsorted_4150/spec_09 .pm 
2  -t  "PYRIDINE,  144  ppm,  T=25  deg  C"  -L  HyGAS_Sprouse_80  -n  spec_10r  . . /unsorted_4150/spec_10r .prn 


==>  do_add21ib_Sprouse_6  5  <= 
add_spectrum  -f  2  - t  "METHYL 
add_spectrum  -f  2  - t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  - t 
add„spec tr urn  -f  2  -t 
add_spectrum  -f  2  - t 
add_spectrum  -f  2  -t 


'METHYL  CHLORIDE,  800  ppm"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_01  . . /scaled_0_to_l/spec_01 
■ETHYL  CHLORIDE,  800  ppm"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_02  . . /scaled_0_to_l/spec_02 
"ISOPROPYL  CHLORIDE,  49.54  ppm'  -L  HyGAS_Rauss-Sprouse_65  -n  spec_03  . . /scaled_0_to_l/spec_03 
"DIMETHYL  ETHER,  200  ppm"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_04  . . /scaled_0_to_l/spec_04 
"DIMETHYL  DISULFIDE,"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_05  . . /scaled_0_tO_l/spec_05 
■PINACOLYL  ALCOHOL,  193.50  ppm"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_06  . . /scaled_0_to_l/spec_06 
"2-CHLOROETHANOL,"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_07  . . /scaled_0_to_l/spec_07 

■FORMAMIDE,  404  ppm,  T=60  deg  C"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_08r  . . /scaled_0_to_l/spec_08r 
•TRIMETHYL  PHOSPHITE,  19.00  ppm'  -L  HyGAS_Rauss-Sprouse_65  -n  spec_09  . . /scaled_0_to_l /spec_09 
"PYRIDINE,  144  ppm,  T=25  deg  C"  -L  HyGAS_Rauss-Sprouse_65  -n  spec_10r  . . /scaled_0„to_l/spec_10r 


==>  do_add2 1 ib [_ERDEC_Emery ]  <== 


add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add_spectrum  -f  2  -t 
add  spectrum  -f  2  -t 


"GA  (LIQUID)  EA0001 
"GA  (VAPOR)  EA0002 • 
"GB  (LIQUID)  EA0003 
"GB  (VAPOR)  EA0004" 
"GD  (LIQUID)  EA0005 
"GD  (VAPOR)  EA0006" 
"VX  (LIQUID)  EA0007 
"CS  (LIQUID)  EA0008 
"CN  (LIQUID)  EA0009 
"HD  (LIQUID)  EA0010 


-L  HyGAS_ERDEC_Emery  -n  GA_liquid  . . /irspecOl .asc 
-L  HyGAS_ERDEC_Emery  -n  GA_vapor  . . /irspec02 .asc 
-L  HyGAS_ERDEC_Emery  -n  GB_liquid  . . /irspec03 . asc 
-L  HyGAS_ERDEC_Emery  -n  GB_vapor  . . /irspec04 . asc 
-L  HyGAS_ERDEC_Emery  -n  GD_liquid  . . /irspec05 . asc 
-L  HyGAS_ERDEC_Emery  -n  GD_vapor  . . /irspec06 . asc 
-L  HyG A S_ERDE C_Eme ry  -n  VX_liguid  . . /irspec07 . asc 
-L  HyGAS_ERDEC_Emery  -n  CS_liquid  . . /irspec08 . asc 
-L  HyGAS_ERDEC_Emery  -n  CN_liquid  . . /irspec09 . asc 
-L  HyGAS_ERDEC_Emery  -n  HD_liquid  . . /irspeclO .asc 


6.  Notes  on  resampling  libraries  for  HyGAS  sensors 

For  HyGAS,  resampling  of  the  spectral  libraries  to  the  various  HyGAS  sensors  was  done 
manually,  because  of  a  problem  with  the  add_sensor  function  in  an  interim  version  of 
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the  HIPS  code.  The  interim  version  did  not  work  correctly  with  the  specific  HyGAS 
spectral  libraries,  which  differed  from  previous  HIPS  libraries  in  ways  such  as  having 
more  spectral  channels,  being  in  cm'1  wavenumber  units  rather  than  pm  wavelength  units, 
and  not  necessarily  being  arranged  in  monotonic  increasing  spectral  order.  To  get  around 
these  difficulties,  a  temporary  library  was  created  with  spectral  resolution  characteristics 
matched  for  each  of  the  desired  HyGAS  sensors.  The  appropriate  matched  sensor  was 
then  added  to  each  temporary  library,  using  the  normal  HIPS  add_sensor  function. 

Since  there  was  no  change  in  spectral  resolution  for  this  contrived  situation  of  exactly 
matched  sensor  and  raw  library  spectrum,  the  HIPS  add_sensor  function  had  no 
problems.  The  sensor  portion  of  each  temporary  library  was  then  moved  to  the 
appropriate  HyGAS  spectral  library,  under  the  appropriate  sensor  name. 

Manual  spectral  resampling  was  done  in  one  of  two  ways.  For  the  case  of  matched 
spectral  resolutions,  differing  simply  in  wavelength  center  positions  (i.e.  for  Emery  and 
Rauss,  both  at  2  cm'1  spectral  resolution),  values  for  the  new  sensor  band  positions  were 
evaluated  using  a  spline  interpolation  function  applied  to  the  original  values.  For  the  case 
of  mismatched  spectral  resolutions  (i.e.,  resampling  Rauss  from  the  original  2  cm'  to  the 
reduced  resolution  of  8  cm'1  for  the  Kroutil  sensor),  a  simple  average  of  four  consecutive 
datapoints  was  performed.  For  this  case,  every  four  x-values  were  averaged,  giving  a 
resultant  8  cm'1  resolution  band  centered  on  the  Kroutil  sensor  band  centers.  For  use  in  a 
spectral  search  library,  the  four  corresponding  y-values  can  nearly  equivalently  be  either 
summed  (integrated)  or  averaged  without  altering  the  spectral  signature.  With  HyGAS, 
and  purely  for  ease  of  comparison  with  the  original  spectra  when  doing  HIPS  spectral 
plots,  we  chose  to  use  the  y-value  averages  rather  than  the  integral  sums. 

Numerous  batch  command  files  were  created  to  perform  the  spectral  resampling  steps 
required  for  the  various  sensors.  These  files  are  located  in  the  directories  referenced  for 
Tables  7  and  8.  These  files  can  be  executed  as  UNIX  shell  scripts.  Any  additional 
executables  required  by  these  command  files,  along  with  associated  source  code,  are 
provided  with  HyGAS.  Executables  and  source  files  are  located  in  the  -pamo/bin  and 
-pamo/src  directories,  respectively. 
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HyGAS  Simulation  Image  Cubes 

Image  cubes  with  sufficient  spectral  and  spatial  resolution  for  exploitation  of  the  standoff 
gas  scenarios  selected  for  HyGAS  were  not  available  at  the  time  of  this  effort.  Therefore, 
we  simulated  appropriate  image  cubes  for  testing,  validation  and  demonstration  of  the 
HyGAS  software.  (Use  of  these  simulation  image  cubes  with  the  various  HyGAS 
software  routines,  demonstrating  HyGAS  in  various  application  scenarios,  is  discussed  in 
detail  in  the  HyGAS  User's  Manual/T utorial,  Section  9,  Interactive  HyGAS.) 

The  HyGAS  simulation  image  cubes  were  generated  by  assigning  appropriate  spectral 
information  for  a  simulated  HyGAS  sensor  to  each  pixel  in  a  simple  two-dimensional 
spatial  scene.  Included  in  the  scene  would  be  pixels  with,  and  without,  the  spectral 
signature  of  various  gases.  [We  could  have  constructed  a  simplistic  block  diagram  image 
(as  was  done  for  an  earlier  synthetic  cube  project  for  this  customer);  however,  we  elected 
to  create  a  more  visually  appealing  and  spectrally  realistic  image  by  using  a  real  dataset 
from  a  multispectral  imaging  sensor,  and  assigning  simulated  hyperspectral  signatures 
instead  of  the  actual  multispectral  signatures  to  the  various  materials  present  in  the 
scene.] 

1.  Selecting  a  spatial  scene. 

The  spatial  scene  we  selected  was  a  Landsat  image  of  an  unpopulated  and  sparsely 
vegetated  island,  consisting  primarily  of  natural  components,  i.e.  soil,  vegetation  and 
water.  From  our  existing  library  of  local  area  Landsat  images,  we  selected  the  small 
Hawaiian  island  of  Kahoolawe  (Figure  1).  We  selected  a  relatively  cloud-free  subset  of 
this  image,  over  the  southeast  portion  of  the  island  (boxed  area  in  Figure  1). 

2.  Defining  the  materials  present  in  the  scene. 

We  created  a  multispectral  signature  model  (a  spectral  mixture  model)  of  this  scene  using 
a  spectral  Linear  Mixing  Model  (LMM)  algorithm.  We  processed  the  selected  subset 
image  of  SE  Kahoolawe,  with  six  Landsat  spectral  bands  (visible  through  near  infrared), 
using  the  HIPS  LMM  algorithm  and  the  interactive  HIPS  function  STARS.1  The 
interactive  STARS  procedure  allows  the  user  to  select  appropriate  spectral  endmembers 
representative  of  the  various  spectrally  distinct  materials  available  in  the  scene.  The 
function  then  generates  a  corresponding  model  of  the  scene  by  combining  solely  these 
endmembers  in  a  linearly  additive  fashion.  For  each  spectral  endmember,  a  fractional 
image  plane  is  generated,  representing  for  each  pixel  the  percentage  composition 
attributed  to  that  particular  endmember.  A  sum  of  all  the  fractional  image  planes  should 
total  100  percent  for  each  individual  pixel.  A  residual  image  also  is  calculated,  to  display 
the  difference  between  the  calculated  model  and  the  original  dataset. 

A  common  method  of  selecting  spectral  endmembers  is  to  start  with  two  endmembers, 
review  the  resulting  residuals  to  determine  how  best  to  revise  the  existing  endmembers  or 
add  new  endmembers,  and  then  iterate  on  this  procedure.  A  model  with  the  fewest 
endmembers  required  to  yield  a  reasonably  randomly  patterned  residual  image  is  the 
recommended  iteration  stopping  point.  The  maximum  number  of  endmembers  cannot 
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exceed  the  quantity  one  greater  than  the  number  of  available  spectral  bands;  however, 
fewer  endmembers  can  often  model  the  scene  adequately. 

Typically,  a  representative  bright  and  a  representative  dark  pixel  are  selected  as  the  first 
two  spectral  endmembers.  These  are  often  the  pixel  uniformly  brightest,  and  the  pixel 
uniformly  darkest,  in  all  available  spectral  channels.  To  a  first  order,  the  uniformly  dark 
pixel  represents  a  basic  shade  component,  which  can  be  added  in  various  percentages  to 
any  pixel,  to  change  its  overall  apparent  albedo  without  altering  its  spectral  signature. 
The  uniformly  bright  pixel,  again  to  a  first  order,  represents  any  dominant  spectral 
signature,  if  only  in  terms  of  being  the  brightest  albedo.  (In  a  case  such  as  ours,  where 
there  is  contamination  by,  for  example,  a  bright  cloud,  a  uniformly  bright  pixel  in  an 
uncontaminated  region  should  be  selected  instead  of  the  contaminated  pixel.) 


Figure  1.  Landsat  view  of  Kahoolawe,  5  January  1989.  This  is  a  false-color  composite  RGB  image  of 
Kahoolawe  (a  dry,  uninhabited  island  in  the  Hawaiian  Island  chain).  Three  Landsat  bands  (bands  4,  5  and 
2,  arbitrarily  scaled  from  12-50,  5-70  and  20-100  DNs,  respectively)  are  combined  to  form  the  false-color 
RGB  image.  For  the  simulated  HyGAS  image  cube,  we  selected  the  relatively  cloud-free  SE  area  of  the 
island  (red  box,  240  x  150  pixels).  This  area  is  composed  primarily  of  natural  elements —  soil,  vegetation 
and  water. 
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Figure  2  (left).  Fraction  planes  from  Linear 
Mixing  Model  (LMM)  representing  the  various 
spectral  endmembers  selected  to  model  the 
simulated  image  cube  for  HyGAS.  The  LMM 

function  in  HIPS  was  used  to  select  representative 
spectral  endmembers  for  this  scene.  From  top  to 
bottom  are  fraction  planes  representing  (1)  bright, 
red  soil,  (2)  dark  water/shade,  (3)  generic 
vegetation,  (4)  bright  shoreline  component,  and  (5) 
an  additional  darker  soil  or  vegetation  component. 
In  each  fraction  plane,  bright  represents  a  higher 
percentage  of  the  specific  material.  (Each  fraction 
plane  is  displayed  over  the  same  range,  0-1  display 
minimum  to  display  maximum,  representing  0%  to 
100%  fractional  composition.)  Although  not 
necessarily  representative  of  physical  reality,  these 
endmembers  represent  an  approximation  suitable 
for  creating  a  synthetic  image  cube  for  use  with 
HyGAS.  Note  the  cloud  and  its  shadow  (apparent 
in  fraction  plane  2  as  a  dark  patch  and  a  bright 
patch,  respectively),  and  how  it  appears  as  a  bright 
or  dark  anomaly  in  each  fraction  plane. 


Figure  3.  False-color  fraction  plane  image  for 
simulated  HyGAS  image  cube.  For  this 
composite,  the  first  three  fraction  plane  images 
(soil,  vegetation  and  water)  were  combined  into  a 
single  RGB  color  image.  A  fraction  plane  image 
such  as  this  is  useful  to  emphasize  the  spatial 
distribution  of  the  various  endmembers.  Since  only 
three  endmembers  can  be  combined  at  any  one  time 
in  a  simple  RGB  image,  some  false  impressions  can 
occur  because  of  the  excluded  endmembers. 
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For  our  scene,  we  selected  our  first  two  endmembers  as  (1)  a  bright  non-cloud  pixel  in  a 
region  that  appeared  in  the  false  color  RGB  image  as  reddish  soil  (red  kaolinitic  clay  soils 
are  very  common  in  the  arid  upland  plateau  regions  of  Hawaii),  and  (2)  the  uniformly 
darkest  pixel  in  the  scene.  Iteratively  selecting  endmembers,  based  on  the  appearance  of 
the  residual  image  and  an  a  priori  assumption  about  the  native  materials  appropriate  for 
the  scene  (various  soils,  sparse  vegetation  and  water),  we  selected  three  additional 
endmembers  to  model  our  specific  scene:  (3)  a  likely  vegetation  endmember,  occurring 
on  upper  slopes  and  in  valleys,  (4)  a  bright  shoreline  component  (possibly  beach  sand, 
white  water  or  sea  foam),  and  (5)  a  darker  lower  slopes  component  (most  likely,  a 
different  type  of  soil  or  a  different  type  of  vegetation).  Fractional  image  planes  for  all 
five  endmembers  are  shown  in  Figure  2. 

A  false-color  RGB  composite  image  using  the  first  three  endmembers  (soil,  vegetation 
and  water)  is  shown  in  Figure  3.  Fraction  plane  images  displayed  in  this  manner  are  often 
used  to  emphasize  the  spatial  location  or  spatial  distribution  of  one  or  more  specific 
endmembers.  While  not  necessarily  a  complete  picture  (using  only  three  colors  to 
represent  five  endmembers),  this  picture  certainly  highlights  the  regions  of  red  soil,  green 
vegetation  and  blue  water. 

3.  Defining  simulated  hyperspectral  signatures  for  selected  materials 

Once  the  spectral  endmembers  have  been  determined  from  the  multispectral  scene, 
appropriate  hyperspectral  signatures  must  be  defined  for  each  endmember.  A  recombined 
image  cube  can  then  be  constructed  using  this  new  spectral  information  and  the 
appropriate  fraction  plane  images  for  each  spectral  endmember.  The  new  image  cube  will 
contain  the  same  spatial  characteristics  as  the  Landsat  image,  but  will  have  a  simulated 
hyperspectral  sensor  response  function  appropriate  for  gas  applications,  instead  of  the 
original  broadband  response  functions  of  the  six  Landsat  bands. 

Simulated  hyperspectral  signatures  for  each  endmember  were  approximated  to  match 
both  types  of  field  data  supplied  to  us  (GFI)  for  use  with  HyGAS,  namely  (1)  the 
Rauss_Howden  field  data,  and  (2)  the  Kroutil_MIDAC  field  data. 

The  Rauss_Howden  field  data  [presumably]  came  from  an  experiment  by  E.  Howden  of 
NVEOD,  using  a  point  sensor  FTIR  spectrometer.  An  example  of  selected  background 
spectra  from  this  field  spectrometer  is  shown  in  Figure  4.  These  spectra  cover  the 
spectral  range  833-1265  cm'1  at  2  cm'1  resolution,  with  449  data  points  each. 

The  Kroutil_MIDAC  field  data  came  from  an  experiment  by  R.  Kroutil  of  ERDEC,  using 
a  MID  AC  point  sensor  FTIR  spectrometer.  An  example  of  a  background  interferogram 
from  this  spectrometer  (prior  to  performing  the  FFT)  is  shown  in  Figure  5.  These 
interferograms  contain  1024  data  points  each.  An  example  of  the  corresponding  spectrum 
(calculated  by  performing  an  FFT  on  the  example  interferogram)  is  shown  in  Figure  6. 
These  spectra  cover  the  spectral  range  0-1971  cm'1  at  8  cm'1  resolution,  with  512  data 
points  each.  [Much  of  this  spectral  range,  however,  is  outside  the  region  of  sensitivity  of 
the  sensor  (roughly  600-1600  cm'1).  For  example,  only  1 13  data  points  fall  within  the 
corresponding  spectral  region  (833-1265  cm'1)  for  the  Rauss_Howden  sensor.] 
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Figure  4.  Example  Rauss_Howden  field  sensor  datasets.  Example  background  datasets  (bgl-bg4)  from 
the  Rauss-Howden  field  sensor  (from  experiments  on  12  August  1992)  are  shown.  In  general,  these  spectra 
represent  the  overall  blackbody  radiation  emitted  by  the  background,  uncontaminated  by  any  target  gases. 
Spectral  absorptions  caused  by  the  intervening  atmosphere  in  the  optical  path  between  the  sensor  and  the 
background  are  present,  however.  All  spectra  were  supplied  (GFI)  in  SpectraCalc  format,  as  449-data  point 
spectra  over  the  spectral  range  833-1265  cm'1  wavenumbers,  at  2  cm'1  resolution  (Nyquist  sampling  of  1 
cm"1  data  point  spacing);  spectra  were  translated  into  JCAMP  and  ASCII_XY  format  using  the  GRAMS 
software  package,  for  incorporation  into  HyGAS. 

KroutiLMIDAC  Field  Data 
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Interferogram  Channel  (1-1024) 

Figure  5.  Example  KroutiLMIDAC  field  sensor  dataset  in  original  interferogram  domain.  An 

example  background  dataset  (mdtOl,  record  #1)  from  the  KroutiLMIDAC  field  sensor  (from  experiments 
in  February,  of  undetermined  year)  is  shown.  This  is  an  example  interferogram  from  an  FTIR  spectrometer, 
showing  the  strong  centerburst  intensity  response  typical  for  a  broadband  thermal  source.  The  broadband 
source  in  this  case  is  the  overall  blackbody  radiation  emitted  by  the  background,  uncontaminated  by  any 
target  gases.  These  datasets  (1024-data  point  interferograms)  were  supplied  in  binary  format,  readable  by 
the  pcida  software  supplied  (GFI)  by  Kroutil.  (The  corresponding  spectrum  generated  from  this 
interferogram  is  shown  in  Figure  6.) 
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KroutiLMIDAC  Field  Data 


Figure  6.  Example  Kroutil-MEDAC  field  sensor  dataset  in  transformed  (FFT)  spectral  domain.  The 

spectrum  shown  here  was  generated  from  the  interferogram  above  (Figure  5),  using  the  pcida  software 
supplied  (GFI)  by  Kroutil.  This  spectrum  represents  the  overall  blackbody  radiation  emitted  by  the 
background,  uncontaminated  by  any  target  gases.  Spectral  absorptions  caused  by  the  intervening 
atmosphere  in  the  optical  path  between  the  sensor  and  the  background  are  present,  however.  Interferograms 
from  this  sensor  are  converted  by  pcida  into  spectra  with  512  data  points  covering  the  spectral  range  0- 
1971  cm'1  wavenumbers,  at  8  cm'  resolution  (Nyquist  sampling  of  4  cm'1  data  point  spacing).  The  spectral 
region  below  500  cm'1  falls  outside  the  region  of  the  responsivity  of  the  sensor,  and  therefore  lies  near  zero 
intensity. 

The  datasets  for  both  experiments  were  supplied  to  us  with  no  additional  information 
regarding  the  sensor  or  the  experimental  setup.  Given  this  situation,  we  made  several 
assumptions  in  order  to  use  these  datasets  as  the  basis  to  generate  simulated  hyperspectral 
signatures  for  the  selected  spectral  endmembers. 

For  each  experimental  setup,  the  basic  measurement  setup  represents  passive  thermal 
measurements  of  the  respective  field  background.  (Passive  measurements  use  no  added 
artificial  source;  the  radiating  field  background  is  the  thermal  source.)  These 
measurements,  to  a  first  order,  are  simply  a  function  of  the  thermal  background 
temperature  (modified  by  the  emissivity  characteristics  of  the  specific  background 
material),  the  transmission  characteristics  of  the  intervening  atmosphere  between  the 
source  (the  thermal  background)  and  the  sensor,  and  the  overall  spectral  response 
function  of  the  sensor.  [No  target  gases  are  present  in  the  intervening  optical  path  for  any 
dataset  selected  as  a  field  background;  however,  the  intervening  atmospheric  constituent 
gases  (primarily  nitrogen,  oxygen,  argon  and  water  vapor,  along  with  a  host  of  minor 
constituents)  function  as  absorbing  gases  since  the  atmosphere  is  usually  slightly  cooler 
than  the  ambient  field  surface  background.] 

If  the  spectral  response  function  of  the  sensor  (i.e.,  gain  and  offset)  and  the  characteristics 
of  the  intervening  atmosphere  (i.e.,  optical  path-length)  are  not  known,  their  combined 
effect  can  be  approximated  by  estimating  the  temperature  and  the  emissivity  of  the 
background  and  comparing  these  values  to  an  actual  background  measurement.  [The 
spectral  emissivity  of  a  surface  is  defined  as  the  ratio  of  the  radiance  emitted  by  the 
surface  to  that  emitted  by  a  blackbody  of  the  same  temperature.] 
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Essentially,  the  measured  data  number  (DN)  generated  by  the  instrument  is  the  sum  of  the 
instrumental  offset  (bias)  and  the  product  of  the  instrumental  gain  and  the  radiance  at  the 
sensor  ( Ls ). 

DN  =  offset  +  (gain  x  Ls)  ( 1 ) 

The  radiance  at  the  sensor  is  the  product  of  the  emissivity  ( es )  of  the  surface  material,  the 
blackbody  radiance  ( Bs )  of  the  same  temperature,  and  the  transmission  ( zatm  )  of  the 
intervening  atmosphere. 


Ls=£sxBsxratm 


(2) 


This  translates  to 


DN  =  offset  +  (gain  xesxBsx  zatm )  (3) 

or,  combining  the  factors  representing  the  sensor  response  function  and  the  atmosphere, 

gain  x  zmm  =  (DN  -  offset )/(ej  x  Bs )  (4) 

The  instrumental  offset  (usually  a  bias  voltage  applied  to  the  A/D  converter  to  register  a 
very  small  positive  voltage  under  conditions  of  zero  signal  input)  can  be  assumed  to  be 
roughly  constant  for  all  wavelengths,  for  a  constant  bandwidth  sensor.  This  value  can  be 
estimated  from  the  instrumental  output  (DN)  from  a  spectral  region  with  little  or  no  signal 
received  at  the  sensor  (i.e.  shortwards  of  500  cm'1  in  Figure  6).  If  such  data  is  not 
available,  as  for  the  RaussJHowden  sensor,  this  offset  can  only  be  assumed  to  be 
negligible.  The  combined  atmosphere-sensor  factor  then  simplifies  to 

gain  x  xatm  =  DN/^  xBs)  (5) 

For  an  opaque  material  (absorption  only,  no  transmission),  the  spectral  emissivity  ( es )  is 
related  to  the  hemispherical  spectral  reflectance  ( ps ),  by  means  of  Kirchhoff  s  law, 

£,  =  1  -  Ps  (6) 

Substituting  Eq.  6  into  Eq.  5  yields 

gain  x  talm  =  DN/((1  -ps)xBs)  (7) 

The  DN  in  Eq.  7  is  obtained  from  the  spectral  measurements  of  the  field  background 
(Figure  4  and  Figure  6). 

The  field  background  is  assumed  to  be  grass,  which  can  be  represented  by  any  generic 
grass  for  which  the  spectral  reflectance  properties  ( ps )  are  known  (Figure  7). 

The  temperature  for  Eq.  7  is  estimated  qualitatively  by  comparing  the  raw  spectral 
measurements  of  the  field  background  to  a  Planck  blackbody  function  at  a  temperature 
appropriate  for  each  experiment.  The  RaussJHowden  dataset  was  collected  during  the 
summer;  the  Kroutil_MIDAC  dataset  was  collected  during  the  winter.  Figure  8  shows 
the  relative  spectral  emittance  for  a  blackbody  at  an  estimated  summer  temperature  of 
300K  (80°  F,  27°  C)  and  at  an  estimated  winter  temperature  of  273K  (32°  F,  0°  C). 

Plotted  on  the  same  scale  are  the  calculated  spectral  radiance  for  grass  at  each 
temperature,  and  the  corresponding  field  backgrounds  for  the  Rauss_Howden  and 
Kroutil_MIDAC  experiments — all  intensities  are  scaled  relative  to  the  peak  spectral 
radiance  for  a  blackbody  at  300K. 
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Directional  Hemispherical  Reflectance 


Wavelength  (microns) 

Figure  7.  Hemispherical  reflectance  for  grass.  The  directional  hemispherical  spectral  reflectance  for 
grass  was  provided  by  J.  Salisbury  (private  communication).  This  generic  grass  can  be  used  as  a  rough 
approximation  of  the  background  material  for  each  of  the  field  experiments,  for  purposes  of  determining  the 
emissivity  of  the  field  background  from  the  hemispherical  reflectance  using  Kirchhoff  s  law. 


Comparison  of  Spectral  Emissivity  Approximations 
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Figure  8.  Comparison  of  Spectral  Emissivity  Approximations.  The  spectral  radiance  for  a  blackbody  at 
300K  and  at  273K  are  scaled  to  the  maximum  spectral  radiance  for  the  blackbody  at  300K.  These 
temperatures,  selected  as  appropriate  for  summer  and  winter,  best  matched  a  rough  convex-hull 
approximation  of  the  raw  background  spectra  from  the  RaussJHowden  and  the  Kroutil_MIDAC 
experiments.  The  spectral  radiance  for  grass  at  these  same  temperatures  also  is  plotted,  similarly  scaled  to 
the  blackbody  at  300K.  The  raw  background  spectra  are  arbitrarily  scaled  to  match  the  appropriate 
temperature  grass  background. 

These  temperatures  were  selected  to  match  the  shape  of  the  background  spectra  near  the 
peak  spectral  response — we  estimated  a  convex-hull  approximation  for  the  raw  spectra 
over  a  small  region  near  the  peak  response,  selected  a  blackbody  temperature  appropriate 
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for  the  season,  and  then  altered  this  temperature  by  increments  of  5°  C  to  select  the 
temperature  with  the  best  match  to  this  shape.  The  initial  temperature  estimates  of  300K 
and  273K  appeared  to  be  the  closest  matches.  [Theoretically,  the  Wien  Displacement  Law 
can  be  used  to  determine  the  temperature  of  a  blackbody  source  from  the  wavelength 
corresponding  to  the  peak  of  the  spectral  radiation  curve;  however  in  practice,  the  spectral 
shape  of  the  atmospheric  transmittance  in  this  spectral  region,  along  with  the  spectral 
response  function  of  the  sensor,  can  skew  the  wavelength  position  of  the  apparent 
maximum  of  the  measured  radiation,  resulting  in  erroneous  estimates  of  temperature.] 

As  shown  in  Figure  8,  the  deviation  between  a  blackbody  emitter  and  the  actual  measured 
background  is  more  a  function  of  the  transmission  of  the  atmosphere  and  the  sensor 
response  than  of  the  emissivity  of  the  background  material.  The  ratio  of  the  measured 
DN  to  the  calculated  radiance  of  the  background  at  the  appropriate  temperature,  i.e.  the 
combined  atmosphere-sensor  factor  as  specified  by  Eq.  7,  is  shown  for  the 
Rauss_Howden  simulation  in  Figure  9,  and  for  the  Kroutil_MIDAC  simulation  in  Figure 
10.  Figure  1 1  superposes  both  factors  on  a  single  plot,  scaled  to  unity  at  the  maximum 
value. 


HyGAS  Simulation:  Assumed  Hauss— Howden  Senior  HyGAS  Simulation:  Assumed  KrauUI— MIDAC  Sensor 


Figure  9.  Combined  sensor-atmosphere  factor 
for  Rauss_Howden  simulation.  This  factor  was 
calculated  using  Equation  7,  substituting  the 
measured  DN  of  bg3  (Figure  4),  the  spectral 
reflectance  of  grass  (Figure  7),  and  an  estimated 
temperature  of  300K. 


Figure  10.  Combined  sensor-atmosphere  factor  for 
KroutilJMIDAC  simulation.  This  factor  was 
calculated  using  Equations  4  and  6,  substituting  the 
measured  DN  of  mdtoi_oi  (Figure  6),  an  estimated 
offset  value  of  75  DN,  the  spectral  reflectance  of  grass 
(Figure  7),  and  an  estimated  temperature  of  273K. 
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Figure  11.  Comparison  of  combined  sensor-atmosphere  factors  for  both  simulations.  The  plots  of 
Figures  9  and  10  are  scaled  to  unity  at  maximum  response  and  superposed  for  comparison. 
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Once  the  combined  sensor-atmosphere  factor  (and  offset,  if  appropriate)  has  been 
determined  for  each  experiment,  estimated  DN  values  (i.e.  hyperspectral  signatures)  at 
whatever  temperatures  are  desired  can  be  calculated  for  any  material  whose  spectral 
reflectance  is  known,  by  combining  Equations  3  and  6. 

Spectral  reflectance  values  for  the  five  materials  selected  to  represent  the  five 
endmembers  are  shown  in  Figure  12.  Corresponding  emissivity  values  (over  the  8-12 
pm  spectral  region  only)  are  detailed  in  Figure  13. 

Candidate  Materials  for  HyGAS  Simulation  Image  Cubes 


Wavelength  (microns) 

Figure  12.  Hemispherical  reflectance  values  for  spectral  endmembers  selected  for  use  in  HyGAS 
simulation  image  cubes.  Candidate  materials  with  existing  hemispherical  reflectance  spectra  were  selected 
to  represent  each  of  the  five  spectral  endmembers  defined  from  the  Landsat  multispectrai  image  subset  of 
Kahoolawe.  These  materials  were  selected  from  various  in-house  databases.  Three  materials  came  from 
the  Spectral  Library:  red  clay  (0150UUUSOL),  sand  (0454UUUSOL)  and  brown  soil  (0216UUUSOL), 
representing  endmember  1,  bright  red  soil;  endmember  4,  bright  shoreline  component;  and  endmember  5, 
darker  soil.  [Sample  0150UUUSOL  is  a  reddish  brown  silty  clay,  representing  a  suitable  analog  for  the 
bright  red  kaolinitic  clays  characteristic  of  the  arid  upland  plateau  regions  of  Hawaii.  Sample 
0454UUUSOL  is  a  light  brown  sandy  soil  with  high  carbonate  content,  representing  a  suitable  analog  for  a 
calcareous  beach  sand  mixed  with  kaolinitic  clay.  Sample  0216UUUSOL  is  a  gray  silty  clay  with  some 
montmorillonitic  and  kaolinitic  clays  mixed  with  some  minor  darker  organic  material;  this  represents  a 
suitable  analog  for  a  soil  near  vegetation.  Although  these  samples  were  selected  for  their  spectral 
variability  and  approximate  composition,  their  actual  composition  and  relation  to  physical  reality  have  no 
impact  on  their  function  in  the  simulated  image  cube.]  Water  and  vegetation  came  from  separate  databases 
(GFI  from  J.  Salisbury),  for  endmembers  2,  water;  and  endmember  3,  vegetation.  [The  water  spectrum  we 
used  is  the  same  dataset  available  in  the  NEF  database;  the  vegetation  spectrum  is  a  sample  of  sugar  maple. 
Again,  the  actual  composition  of  the  vegetation  and  its  relation  to  physical  reality  have  no  impact  on  the 
function  of  vegetation  in  the  simulated  image  cube.] 


For  simulated  image  cubes  for  HyGAS,  we  selected  an  estimated  base  temperature  of 
300K  for  all  materials  except  water/shade,  which  we  selected  to  be  5°  C  cooler  (295K). 
Calculated  radiance  values  using  these  temperatures  are  shown  for  all  five  endmembers 
(along  with  a  reference  blackbody  at  300K)  in  Figure  14.  [For  HyGAS  simulations,  we 
made  the  gross  assumption  that  the  combined  factor  representing  the  sensor  response  and 
the  atmosphere  is  constant  across  the  simulated  scene.]  These  same  materials  are  shown 
corrected  to  the  sensor-atmosphere  function  for  a  RaussJHowden  sensor  (Figure  15)  and 
a  Kroutil_MIDAC  sensor  (Figure  16).  Corresponding  interferograms  for  the 
KroutiLMIDAC  sensor  are  shown  at  full  scale  in  Figure  17  and  expanded  in  Figure  18. 
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Figure  13.  Emissivity  values  for  selected 
endmember  materials.  The  emissivity  values  are 
derived  from  the  spectral  reflectance  values 
shown  in  Figure  12.  Only  the  8-12  pm  spectral 
region  is  displayed. 
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Figure  14.  Calculated  radiance  for  endmember 
materials  for  HyGAS  simulations.  All  materials 
were  assumed  to  be  300K  except  water/shade  at 
295K.  A  300K  blackbody  is  plotted  for  reference. 
All  values  are  scaled  to  unity  at  the  peak  response  for 
the  300K  blackbody  emitter. 
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Figure  15.  HyGAS  endmembers  for  simulated  Figure  16.  HyGAS  endmembers  for  simulated 
Rauss_Howden  sensor.  KroutiLMIDAC  sensor,  spectral  domain. 
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Figure  18.  HyGAS  endmembers  for  simulated  KroutiI_MIDAC  sensor,  interferogram  domain 
expanded.  Differences  in  the  interferograms  for  endmembers  are  slightly  more  distinguishable  further  out 
in  the  interferogram  wings  (right)  than  near  the  centerburst  (left). 


It  should  be  noted  that  these  interferograms  (Figures  17  and  18)  are  modeled  with  a 
centerburst  response  at  channel  48,  compared  with  a  centerburst  response  at  channel  1 16 
for  the  earlier  interferogram  example  (Figure  5).  A  centerburst  at  channel  48  was 
selected  to  match  the  dataset  kft001.dat  supplied  for  testing  with  associated  linear 
discriminants.  [The  discriminant  test  operates  on  a  small  subset  of  the  full  1024-point 
interferograms;  this  subset  is  specified  in  channel  space  limits  relative  to  the  centerburst 
channel.  For  this  reason,  the  position  of  the  desired  output  centerburst  is  one  of  the  input 
specifications  required  for  the  spectrum-to-interferogram  function  in  HyGAS.] 

It  also  should  be  noted  that  these  calculated  interferograms  (Figure  17, 18),  representing 
ideal  interferograms,  are  symmetric  about  the  centerburst,  unlike  actual  measured 
interferograms  (Figure  5),  which  are  asymmetric.  This  measurement  error,  called  phase 
error  or  sampling  error,  is  usually  corrected  by  comparing  the  shape  of  a  measured 
interferogram  on  both  sides  of  the  centerburst,  and  estimating  a  corrected  interferogram 
that  is  truly  symmetric  about  the  centerburst.  This  symmetric  interferogram  is  then 
processed  to  generate  the  corresponding  spectrum.  For  HyGAS  simulations,  we  used  the 
calculated  symmetric'  interferograms,  without  any  simulated  sampling  error. 


4.  Recombining  hyperspectral  endmember  signatures  into  a  synthetic  image  cube 


Once  hyperspectral  signatures  were  generated  for  all  five  endmembers  for  both  the 
Rauss_Howden  and  the  Kroutil_MIDAC  sensors  (Figure  15  and  Figure  17),  we 
proceeded  to  create  a  simulated  HyGAS  hyperspectral  image  cube  for  each  sensor  by 
recombining  the  appropriate  set  of  hyperspectral  endmembers  according  to  the  previously 
defined  fraction  planes  (Figure  2).  A  special  command-line-only  HIPS  function  ( simcub ) 
was  created  to  generate  these  image  cubes,  given  an  appropriate  hyperspectral 
endmember  file  for  each  sensor  and  a  fraction  plane  image  cube.  [Command  files, 
executable  as  UNIX  shell  scripts,  to  generate  these  cubes  were  supplied  with  HyGAS — 
do_simcub .  rauss  and  do_simcub.kroutil  in  the  directory 
-HyGAS  /  s  imulat  ions/ cubes /xs  tuff.] 


To  simulate  potential  variations  due  to  instrumental  noise,  next  we  estimated  a  noise  level 
from  the  DN  variation  between  two  comparable  background  measurements  for  each 
sensor.  For  the  Rauss_Howden  sensor,  we  estimated  the  average  (-peak  to  peak) 
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difference  in  spectral  domain  at  ±0.025  (floating  point  DN  values);  for  the 
Kroutil_MIDAC  sensor,  we  estimated  this  difference  in  interferogram  domain  at  ±4 
(integer  DN  values).  We  then  used  the  create  cube  function  in  HIPS  to  generate  a 
random  number  pattern  "noise"  image  cube,  evaluated  within  this  DN  range,  which  we 
added  to  each  "nojnoise"  image  cube  to  estimate  a  revised  image  cube  for  each  sensor. 

5.  Adding  gas  signatures  into  the  synthetic  image  cubes 

We  created  an  arbitrary  spatial  overlay  map  for  the  SE_Kahoolawe  scene  to  simulate  the 
presence  of  various  gas  clouds.  Four  large  clouds  were  created  with  four  contour  levels 
each,  and  four  smaller  clouds  with  one  contour  level  each.  We  used  the  HIPS  display 
function  create  mask  to  sketch  in  the  proposed  spatial  distributions  for  the  gases  as 
shown  in  Figure  19. 


Figure  19.  Template  for  gas  clouds  for  simulated  HyGAS  image  cube.  Simulated  gas  cloud  contours 
were  created  for  inserting  gases  into  the  simulated  HyGAS  image  cube.  Four  clouds  with  four  contour 
levels,  and  four  smaller  clouds  with  a  single  contour  level,  were  created.  These  contours  represent  a 
simulated  dispersion  with  a  prevailing  wind  pattern  (light  NE  trade  winds).  Table  1  identifies  the  gases 
inserted  into  the  clouds  referenced  by  number  here. 

For  each  sensor,  we  selected  two  gases  to  insert  into  the  scene,  with  each  gas  represented 
under  both  absorbing  (gas  colder  than  the  background)  and  emitting  (gas  hotter  than  the 
background)  conditions.  Three  common  gases  were  used  in  the  RaussJHowden  and  the 
Kroutil_MIDAC  experiments —  sulfur  hexafluoride  (SFe),  acetone  and  methyl-ethyl- 
ketone  (2-butanone,  or  MEK).  All  three  gases  were  released  during  the  RaussJHowden 
field  experiment;  only  SF6  was  released  in  the  Kroutil_MIDAC  field  experiment.  [Sulfur 
hexafluoride  is  a  non-toxic  gas,  commonly  used  as  a  long  wave  IR  marker  gas,  with  a 
single  very  strong  and  narrow  absorption  at  -10.6  pm  (940  cm"1).  Acetone  and  MEK  are 
common  solvents  with  relatively  low  vapor  pressure,  also  with  fairly  strong  and  narrow 
vapor  phase  absorptions,  at  -8-9  pm  (1150-1250  cm'1).  Acetone  also  is  used  in  copious 
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amounts  in  the  chemical  production  of  cocaine,  and  has  related  potential  in  demonstrating 
one  commercial  application  of  HyGAS.] 

5.1  Acetone  and  Sulfur  Hexafluoride 

Since  only  acetone  and  sulfur  hexafluoride  were  included  in  the  Sprouse  spectral 
database,  we  selected  these  two  gases  for  insertion  into  the  scene  for  both  simulations. 
Table  1  indicates  which  gas  was  inserted  into  each  cloud  identified  by  number  in  Figure 
19.  Spectral  absorbance  values  for  acetone  and  sulfur  hexafluoride  (from  the  Sprouse 
database)  are  shown  in  Figure  20. 

Table  1.  Identification  of  gases  inserted  into  each  cloud  numbered  in  Figure  19. 


Large 

cloud 

Small 

cloud 

Gas  Inserted 

Absorption/Emission 

1 

5 

Acetone 

Absorption 

2 

6 

Sulfur  hexafluoride 

Absorption 

3 

7 

Acetone 

Emission 

4 

8 

Sulfur  hexafluoride 

Emission 

Absorbance  of  Acetone  and  Sulfur  Hexafluoride 


Figure  20.  Absorbance  for  Acetone  and  Sulfur  Hexafluoride.  These  laboratory  spectra,  measured  by 
Sprouse  Scientific  Systems  under  contract  to  Rauss,  were  supplied  to  us  (GFI)  for  use  with  HyGAS.  Rauss 
needed  measurements  suitable  for  quantitative  applications,  i.e.,  gases  measured  at  known  concentration,  at 
2cm"1  resolution,  and  at  roughly  ambient  temperature  and  pressure.  These  measurements  were  presumably 
done  using  a  1 -meter  length  gas-cell.  Header  information  included  with  the  data  cite  a  temperature  of  25°  C 
and  concentrations  of  146  ppm  for  acetone  and  5.01  ppm  for  sulfur  hexafluoride. 

For  HyGAS,  we  elected  to  treat  the  Sprouse  absorption  spectra  for  acetone  and  sulfur 
hexafluoride  as  appropriate  relative  absorptions  for  these  two  gases.  We  also  elected  to 
ignore  any  quantitative  implications  related  to  temperature,  pressure,  concentration  and 
optical  path-length,  and  elected  simply  to  consider  the  total  absorbance,  i.e.,  the  combined 
product  of  the  absorption  coefficient,  the  concentration  and  the  optical  path-length  (Beer's 
Law).  In  effect,  we  elected  to  assign  arbitrary  amounts  of  radiance  contributions  because 
of  the  gas,  and  treat  the  gas  as  simply  an  additive  (emission)  or  subtractive  (absorption) 
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effect  on  the  radiance  of  the  underlying  background  reference.  [This  approximation  is  the 
basis  of  the  classic  approach  of  measuring  the  background  reference  with  the  unknown 
target  gas  in  the  intervening  optical  path,  along  with  a  basic  background  reference  with  no 
target  gas  in  the  intervening  optical  path,  and  subtracting  the  measurements  of  the 
background-alone  from  the  background-with-gas  to  calibrate  the  measurements  of  the 
unknown  target  gas.  With  adequate  background  removal,  a  measurement  containing  a 
target  gas  calibrated  in  this  way  approximates  the  laboratory  absorbance  spectrum  of  the 
target  gas  alone.] 

To  approximate  estimated  radiance  contributions  from  gases  for  realistic  scenarios,  and  to 
compare  with  the  supplied  field  data  from  the  Rauss_Howden  experiment  and  the 
Kroutil_MIDAC  experiment,  we  estimated  a  lower  level  of  detectability  at  ~1  percent 
radiance  contribution  because  of  a  target  gas  (Kroutil,  private  communication).  We 
calculated  this  percent  radiance  contribution  for  acetone  and  sulfur  hexafluoride  by 
rescaling  the  Sprouse  database  absorbance  values  for  each  gas  (Figure  20)  to  unity  at  the 
peak  for  SF6,  and  multiplying  these  values  by  the  radiance  DN  values  predicted  for  each 
simulated  dataset.  We  defined  the  result  as  100  percent  radiance  contribution  for  each 
gas.  We  then  multiplied  this  result  by  the  desired  radiance  contribution  (1  percent)  we 
wanted  to  be  attributed  to  the  gas,  and  then  added  (for  emission,  or  subtracted  for 
absorption)  this  result  to  the  original  DN  of  the  background  pixel  for  each  pixel 
contaminated  by  a  gas  cloud  (Figure  19). 

We  inserted  acetone  and  sulfur  hexafluoride  into  each  simulated  cube  at  a  level 
representing  1  percent  radiance  contribution  because  of  gas,  for  the  outermost  contour 
level  of  clouds  1—4,  and  for  all  of  clouds  5-8,  as  specified  in  Table  1.  For  testing 
purposes,  we  specified  increasing  percentages  (10  percent,  50  percent,  even  an  impossibly 
high  100  percent,  purely  for  tests,  assuming  the  percentages  would  be  reduced  to 
something  more  reasonable  at  a  later  date)  for  the  interior  contours  of  clouds  1-4.  [We 
used  the  custom  HyGAS  routine  simmask  to  insert  the  appropriate  spectral  contributions 
defined  for  these  gases.]  Figures  21  and  22  show  representative  examples  of  the  1 
percent,  10  percent  and  50  percent  absorption  contributions  for  both  acetone  and  sulfur 
hexafluoride  for  the  Rauss_Howden  and  the  Kroutil_MIDAC  simulated  image  cubes. 

[The  effect  of  the  impossibly  high  100  percent  contribution  can  be  inferred  from  these 
figures.]  Figure  23  shows  a  section  (channels  100-200)  of  the  corresponding 
interferograms  calculated  from  these  spectra  for  the  Kroutil_MIDAC  simulation. 
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Figure  21.  Effect  of  gas  in  absorption  on 
Rauss  Howden  simulation. 


I  lyOAS  Simulation  Target  Gases:  Kroutil  MIDAC  Sensor 
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Figure  22.  Effect  of  gas  in  absorption  on 
Kroutil_MIDAC  simulation,  spectral  domain. 


HyGAS  Simulotfor  Target  Gases:  Kroutil  MIDAC  Sensor 


hLeiferociair  Chcnne 

Figure  23.  Effect  of  gas  in  absorption  on  KroutiI_MIDAC  simulation,  interferogram  domain. 

The  effect  of  a  narrow  spectral  feature  on  the  outer  regions  of  an  interferogram  is  clearly 
apparent  in  the  strong  sinusoidal  pattern  shown  for  the  50  percent  radiance  contribution 
for  SFg  (Figure  23).  This  is  in  marked  contrast  to  the  identical  interferogram  region  for 
the  same  background  with  acetone  (Figure  24)  and  the  various  individual  backgrounds 
with  no  gas  at  all  (Figure  25).  This  is  the  basis  of  the  linear  discriminant  approach 
developed  by  Kroutil  to  detect  specific  gases  using  specific  sections  of  the  original 
interferogram. 
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HyGAS  Simulation  Target  Gases:  Kroutil  MIDAC  Sensor 


Interferogram  Channel 

Figure  24.  Interferogram  subset  for  background  with  acetone.  The  same  interferogram  channels  are 
displayed  as  in  Figure  23,  but  the  y-axis  is  expanded. 

Candidate  Materials  for  HyGAS  Simulation:  Kroutil  MIDAC  Sensor 


100.0  120.0  140.0  160.0  180.0  200.0 

Interferogram  Channel 

Figure  25.  Interferogram  subset  for  backgrounds  with  no  gas.  The  same  interferogram  channels  are 
displayed  as  in  Figure  23,  but  the  y-axis  is  expanded.  The  y-axis  scale  here  is  the  same  as  for  Figure  24. 

5.2  Methyl-Ethyl-Ketone  (MEK) 

For  the  Kroutil_MIDAC  simulation,  we  investigated  the  possibility  of  substituting  MEK, 
for  which  we  had  a  set  of  discriminant  coefficients,  for  SF6,  for  which  no  discriminants 
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were  available.  MEK  was  not  included  in  the  Sprouse  database,  so  we  attempted  to  get  a 
suitable  MEK  spectrum  from  the  dataset  (kf  toooi .  dat)  provided  to  us  for  testing  with 
the  MEK  discriminants. 


This  dataset  [presumably]  represents  various  laboratory  gas-cell  measurements,  both  with 
and  without  MEK  in  the  optical  path.  We  selected  two  measurements  from  this  dataset 
(records  96  and  97)  to  represent  measurements  of  a  background  with  no  MEK,  and  a 
background  with  MEK  (Figure  26).  A  ratio  of  these  two  measurements  generated  a 
cursory  approximation  of  the  absorbance  of  MEK,  but  with  significant  differences  from 
what  we  expected  for  a  laboratory  measured  absorbance.  We  elected  to  treat  these 
differences  as  an  inadequate  background  removal.  We  then  estimated  a  correction  factor 
as  a  smooth  curve  baseline  approximation,  fit  to  the  points  outside  the  spectral 
absorptions  attributed  to  MEK  (Figure  27).  Removing  this  baseline  approximation 
yielded  the  estimate  of  the  absorbance  spectrum  of  MEK  shown  in  Figure  28. 
[Alternatively,  we  could  have  assumed  the  background  without  MEK  and  the  background 
with  MEK  to  have  been  measured  at  different  temperatures,  and  estimated  a  correction 
based  on  the  estimated  absolute  temperatures;  however,  for  a  first  approximation,  our 
baseline  estimate  correction  appeared  sufficient.] 


Selected  Measurements  with  and  without  MEK 

(spectra  from  interferogroms  in  kft0001.dat) 


WaventJmber  (inverse  cm) 

Figure  26.  Selected  measurements  with  and 
without  MEK  from  file  kftoooi.dat.  These 
spectra  were  computed  from  interferograms  096 
and  097.  Record  096  represents  measurements  of 
the  background  source  with  no  MEK;  record  097 
represents  a  measurement  with  MEK. 


Ratio  of  Measurements  with  and  without  MEK 


Wavenumber  (inverse  cm) 

Figure  27.  Estimate  of  baseline  correction  to 
initial  ratio  of  background  without  MEK  to 
background  with  MEK.  We  elected  to  estimate  a 
baseline  correction  by  fitting  a  smooth  curve  to  the 
points  marked,  which  we  assumed  to  be  outside  the 
absorptions  because  of  MEK. 
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Estimated  MEK  Absorbance 


Wavenumber  (inverse  cm) 

Figure  28.  Estimated  MEK  absorbance  spectrum.  This  estimated  absorbance  spectrum  for  MEK  was 
generated  by  a  ratio  of  spectra  computed  from  interferograms  96  and  97  in  the  test  datafile  kf  toooi .  flat. 
The  remaining  baseline  was  estimated  by  fitting  a  smooth  curve  to  the  points  presumed  to  be  unaffected  by 
the  absorptions  because  of  MEK,  and  removing  this  estimated  baseline  from  the  calculated  ratio.  Although 
noisy  in  the  regions  of  low  original  signal,  we  expect  data  within  the  spectral  region  800-1500  cm  to  be  a 
fair  estimate  of  the  absorbance  of  MEK,  albeit  at  an  arbitrary  y-axis  scale. 

Unfortunately,  the  provided  MEK  discriminants  did  not  perform  very  well  on  any 
interferograms  computed  using  this  estimated  MEK  spectrum.  We  attributed  this  failure 
to  the  strong  apparent  but  erroneous  absorptions  below  800  cm'1,  caused  by  the  low  signal 
levels  of  the  raw  spectrum  in  this  region.  We  tried  setting  the  spectrum  for  MEK  to  zero 
below  800  cm'1,  but  this  failed  to  improve  the  performance  of  the  discriminants.  We  next 
tried  to  approximate  the  spectrum  in  this  region,  and  also  above  1600  cm'1,  by  comparing 
it  with  another  absorbance  spectrum  for  MEK  (Figure  29). 

The  comparison  MEK  spectrum  came  from  the  NIST/EPA  Gas  Phase  Infrared 
Database.2  This  database  is  supplied  as  a  spectral  search  library,  and  as  such  is  arbitrarily 
scaled  (on  the  absorption  axis)  to  unity  at  the  maximum  absorption  for  each  individual 
gas  species.  This  precludes  the  quantitative  use  of  such  a  database,  or  any  comparison  of 
the  relative  intensity  of  absorptions  for  different  species,  although  it  preserves  the 
qualities  required  for  spectral  identification,  namely  the  uniqueness  of  the  positions  of  the 
spectral  absorptions  and  their  relative  intensity  within  any  single  species.  This  database  is 
also  supplied  at  12  cm"1  resolution,  which  would  have  to  be  resampled  to  a  higher 
resolution  (violating  the  Nyquist  rule)  to  be  used  with  the  Kroutil_MIDAC  simulation.  In 
regions  with  no  significant  spectral  absorptions  approaching  the  Nyquist  resolution  limit, 
however,  resampling  to  a  higher  resolution  is  reasonable.  To  correct  the  problem, 
spectral  regions  in  the  estimated  MEK  spectrum  from  kft000l.dat  (Figure  28),  we 
resampled  the  NIST/EPA  spectrum  of  MEK  to  8  cm'1  resolution,  and  substituted  this 
resampled  NIST/EPA  spectrum  for  the  kf  toooi .  dat  spectrum  shortwards  of  850  cm"1 
and  longwards  of  1790  cm'1.  A  comparison  of  the  kf  toooi .  dat  and  NIST/EPA  spectra 
for  MEK,  along  with  the  combined  version,  is  shown  in  Figure  29. 

Unfortunately,  the  supplied  MEK  discriminants  failed  to  perform  adequately  with  any  of 
the  estimated  MEK  spectra.  Based  on  these  failures,  we  abandoned  any  approaches 
intending  to  use  MEK  for  image  cube  simulations. 


2  NIST  Standard  Reference  Database  35,  Version  1.0. 
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Comparison  of  MEK  Absorption  Estimates 


Wavenumber  (inverse  cm) 

Figure  29.  Comparison  of  MEK  Absorbance  Estimates.  Estimates  of  the  absorbance  for  MEK  obtained 
from  the  Kroutil  datafile  kft000l.dat  were  combined  with  NIST/EPA  absorption  values  in  the  outer 
spectral  regions  to  create  a  combined  revised  estimate  of  the  absorbance  for  MEK. 

5.3  Mixtures  of  Gases 

We  created  one  test  cube  for  the  Rauss_Howden  simulation  that  contained  various 
mixtures  of  gases.  For  this  cube,  we  selected  four  gases  representing  various  pesticides 
and  precursors  (or  simulators)  of  chemical  nerve  gas  agents.  We  combined  absorbances 
in  equal  percentages  (50:50  or  25:25:25:25)  as  linearly  additive  spectral  mixtures. 

Mixture  1  was  a  50:50  percentage  mix  of  chloroethanol  and  thiodiglycol,  two  precursors 
to  the  development  of  the  chemical  blister  agent,  mustard  gas.  Absorption  spectra  (from 
the  Sprouse  database)  for  these  two  gases  are  shown  in  Figure  30  (full  scale),  and  in 
Figure  31  (expanded  over  the  800-1300  cm'1  spectral  region). 

Mixture  2  was  a  50:50  mix  of  dimethyl-methyl-phosphonate  (DMMP)  and  diethyl- 
methyl-phosphonite  (DEMP),  two  simulation  agents  for  chemical  nerve  agents. 
Absorption  spectra  (from  the  Sprouse  database)  for  these  two  gases  are  shown  in  Figure 
32  (full  scale),  and  in  Figure  33  (expanded  over  the  800-1300  cm'1  spectral  region). 

Mixture  3  was  a  25:25:25:25  mix  of  all  four  of  these  agents.  Absorption  spectra  for  all 
four  gases  are  superposed  in  Figure  34  (expanded  over  the  800-1300  cm'1  spectral  region 
only). 

Table  2  indicates  into  which  cloud  of  Figure  19  each  mixture  was  inserted.  Each  mixture 
was  inserted  at  the  same  total  radiance  contributions  as  defined  earlier  for  acetone  and 
SFg  (1  percent,  10  percent,  50  percent  and  100  percent  of  the  total  mixture  ratio). 
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HyGAS  Simulation:  Gas  Mixture  Components 
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Figure  30.  Gases  selected  for  Mixture  1,  full 
spectral  range. 


HyGAS  Simulation:  Gas  Mixture  Components 
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Figure  31.  Gases  selected  for  Mixture  1,  800- 
1200  cm'1  only. 


HyGAS  Simulation:  Gas  Mixture  Components  HyGAS  Simulation:  Gas  Mixture  Components 


Figure  32.  Gases  selected  for  Mixture  2,  full  Figure  33.  Gases  selected  for  Mixture  2,  800- 

spectral  range.  1200  cm'  only  . 
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HyGAS  Simulation:  Gas  Mixture  Components 


Wavenumber  (inverse  cm) 

Figure  34.  Gases  selected  for  Mixture  3,  800-1200  cm’1  only.  All  four  gases  were  combined  at  equal 
radiance  contribution. 

Table  2.  Identification  of  gas  mixtures  inserted  into  each  cloud  numbered  in  Figure  19. 


Large 

cloud 

Small 

cloud 

Gas  Mixture  Inserted 

Absorption  or  Emission 

1 

5 

Mixture  1 

(chloroethanol  +  thiodiglycol) 

Emission 

2 

6 

Mixture  2 

(dimethyl-methyl-phosphonate  + 
diethyl-methyl-phosphonite) 

Emission 

3 

7 

Mixture  3 

(all  four  of  the  above) 

Emission 

4 

8 

none 

n/a 

6.  Summary 

A  summary  of  the  simulation  image  cubes  delivered  with  HyGAS  is  given  in  Table  3. 


Table  3.  Summary  of  simulation  image  cubes  delivered  with  HyGAS. 


Image  cube  name 

Gases  Inserted 

Notes 

SE_Kahoolawe_event_locator 

overlay  (value  =  1)  marks  center  of 
each  "gas  event" 

SE__KahooIawe_gas_masks 

overlay  (values  =  1-4,  5-8,  9-12, 
13-16)  for  contour  levels  of  clouds 

SE_Kahoolawe_Kroutil 

Acetone,  SF6 

interferograms 

SE_Kahoolawe_Rauss 

Acetone,  SFg 

spectra 

SE_Kahoolawe_Rauss_2 

none 

spectra 

SE_Kahoolawe_Rauss_3 

Mixtures 

spectra 
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Neural  Net  Notes 

The  implementation  of  a  backpropagation  neural  net  as  the  HyGAS  module  for  spectral 
identification  was  based  on  research  initially  carried  out  by  Patrick  Rauss  at  the  U.S. 

Army  Research  Laboratory.  Rauss  implemented  a  backpropagation  network  that 
successfully  identified  65  spectra  from  Sprouse  Scientific  in  the  presence  of  large 
amounts  (e.g.  200  percent)  of  random  noise.  Rauss's  work  is  fully  described  in  a  paper  he 
submitted  to  the  1994  Army  Science  Conference.  [A  preprint  of  this  paper1  is  attached  at 
the  end  of  this  appendix.] 

1.  Choice  of  Software 

Rauss  implemented  his  neural  net  with  software  from  the  company  HNC,  Inc.  His  choice 
was  made  several  years  ago  when  there  were  fewer  players  in  the  market,  and  was  largely 
based  on  a  desire  to  use  hardware  that  was  specifically  manufactured  for  neural  net 
computations  (HNC's  Balboa  accelerator  board).  In  an  effort  to  be  fully  compatible  with 
Rauss,  and  thus  ease  the  transfer  of  knowledge,  we  initially  sought  to  use  identical  HNC 
software  for  HyGAS.  Despite  much  effort  on  our  part,  and  by  Rauss  (amid  claims  to  the 
contrary  by  HNC,  after  we  had  purchased  the  HNC  software  NeuroSoft ),  we  discovered 
that  the  code  written  by  Rauss  would  not  run  without  purchase  of  the  $10,000  HNC 
Balboa  board.  Furthermore,  after  much  frustration,  it  was  revealed  that  Rauss's  software 
might  not  even  run  on  the  newest  Balboa  board,  having  been  coded  for  an  earlier 
discontinued  version.  This,  plus  the  fact  that  the  HNC  software  was  very  user  unfriendly 
and  complicated,  and  that  the  Rauss  code  was  in  a  research  state  and  essentially 
undocumented,  resulted  in  a  decision  to  look  for  other  software.  Since  backpropagation 
nets  are  the  most  basic  of  neural  nets,  we  expected  a  commercially  available  neural  net 
software  package  would  work  in  place  of  the  HNC  software. 

After  inquiring  on  the  internet,  several  options  were  considered.  One  highly 
recommended  package  was  Matlab  by  The  MathWorks,  Inc.  Our  limited  experience  with 
basic  Matlab  suggested  this  package  was  very  easy  to  use  and  well  supported.  After 
receiving  a  trial  copy  and  a  manual  from  the  MathWorks,  we  determined  that  Matlab' s 
Neural  Net  Toolbox  could  be  used  to  identically  replicate  the  network  that  Rauss  had 
devised.  Initially,  our  intention  was  to  acquire  the  trained  network  (i.e.  weights)  from 
Rauss  and  use  them  in  Matlab.  However,  Rauss  was  unable  to  export  his  weights  out  of 
the  HNC  format.  As  such,  we  undertook  the  added  time-intensive  task  of  training  the 
new  network  in  addition  to  our  original  plan  of  interfacing  a  neural  net  search  function 
into  HyGAS. 

2.  Introduction  to  Backpropagating  Networks 

Backpropagation  nets  are  named  after  back-propagation  of  error,  a  widely  used  learning 
method  for  neural  networks.  Backpropagation  nets  are  very  easy  to  apply  and  useful  in 
many  practical  applications.  However,  it  should  be  noted  that  they  truly  are  a  proverbial 
black  box  in  the  sense  that  you  never  really  know  what  the  net  is  actually  keying  on. 


1  Rauss,  P.  J.,  Noise  Tolerance  of  Model-Based  Neural  ATR,  Army  Science  Conference,  1994. 
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Architecturally,  backpropagating  neural  nets,  like  most  networks,  contain  several  layers: 
an  input  layer,  an  output  layer,  and  one  or  more  hidden  layers.  Each  layer  contains  a 
number  of  nodes,  where  weighted  combinations  of  inputs  from  the  previous  layer  are 
processed  through  a  transfer  function  (typically:  linear,  step  function,  or  log  sigmoidal) 
and  passed  to  the  next  layer.  The  outputs  from  one  layer  are  then  used  as  the  inputs  to  the 
next  layer.  As  such,  the  structure  of  a  backpropagation  net  is  completely  defined  by  its 
size  (number  of  layers  and  number  of  nodes  in  each  layer),  and  the  type  of  transfer 
function  used  at  each  node.  [See  Matlab  manual  pages  2-1  through  2-12.  This  is  an 
excellent  description,  both  visually  and  mathematically,  of  neural  network  architecture.] 

Backpropagation  neural  nets  ar e  feed-forward,  supervised  learning,  networks.  In  order  to 
function,  they  must  be  trained  with  known  inputs  and  known  solutions.  The  network 
weights  are  adjusted  during  this  supervised  training  mode,  using  a  gradient  descent 
algorithm,  which  minimizes  the  errors  between  the  outputs  of  the  network  and  the  known 
solutions.  The  derivatives  of  the  error  between  the  output  of  the  network,  and  the  known 
solutions,  are  backpropagated  through  the  network  as  small  changes  to  each  weighting 
function,  hence  the  name  backpropagation  neural  net.  After  changes  to  the  weights  are 
made,  the  net  is  run  again.  This  process  iterates  until  the  solution  converges  to  an 
expectable  error  (e.g.,  sum  of  the  square  of  the  error  SSE  =  .01).  Once  trained, 
backpropagation  networks  can  very  rapidly  be  applied  to  new  inputs,  in  a  feed-forward 
mode,  with  the  new  inputs  beginning  feed-forward  to  the  network  and  its  trained  weights 
to  produce  the  output  solution. 

Note,  the  solution  to  training  a  neural  network  (i.e.  the  trained  net  weights ),  finds  a 
minimum  in  the  error  surface.  Like  most  iterative  techniques,  this  minimum  is  only 
required  to  be  a  local  minimum,  and  may  not  be  the  global  minimum;  thus,  it  may  not 
necessarily  be  the  best  solution.  Also,  recall  that  the  net  is  minimizing  the  error  relative 
to  all  weights.  With  the  Sprouse_65  spectra,  using  200  hidden  nodes,  and  449  channels, 
this  translates  to  (see  Script  2  below)  (449*200)  +  (200*65)  weights  for  a  total  of  102,080 
parameters. 

3.  Training  a  net. 

In  general  terms,  4  steps  are  required  to  train  and  apply  a  neural  net. 

Step  1 .  Read  in  data 
Step  2.  Set  up  neural  net 
Step  3.  Train  neural  net 
Step  4.  Apply  neural  net 

Definition  of  the  net  all  occurs  in  Step  2.  Because  Matlab  is  a  higher  order  language,  the 
basic  steps  outlined  above  are  very  easily  translated  into  a  Matlab  script  as  shown  in 
Script  1  below. 

3.1.  Example  1:  Training  on  four  spectra  with  60  percent  noise. 
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In  order  to  train,  the  size  and  shape  of  the  net  must  be  determined.  Rauss  had  used  a 
three-layer  neural  net  with  log  sigmoidal  transfer  functions,  65  spectra  (outputs),  and  a 
200  node  hidden  layer.  For  a  four  spectra  net,  the  number  of  inputs  is  449  (the  number  of 
channels  in  the  Rauss_Howden  sensor)  and  the  number  of  outputs  is  four  (four  different 
spectra).  The  only  remaining  choice  is  the  number  of  nodes  in  the  hidden  layer. 

We  chose  20  nodes  for  the  hidden  layer.  Based  on  Rauss's  results  and  the  examples  in  the 
Matlab  Neural  Net  manual,  this  choice  seemed  reasonable.  (It  is  an  order  of  magnitude 
less  hidden  nodes  than  what  was  used  for  the  actual  HyGAS  training  for  the  65  spectra). 
Both  10  and  40  nodes  produced  similar  results,  although  the  40  nodes  seemed  to  take 
longer  to  train.  However,  given  that  no  careful  studies  were  done,  the  choice  of  20  nodes 
is  rather  arbitrary. 

Recall  that  Rauss  was  using  software  from  HNC.  While  the  exact  form  of  the  Rauss 
network  was  used  in  our  application  using  Matlab,  the  Matlab  implementation  for 
training  a  backpropagation  net  is  much  more  sophisticated  than  HNC's.  First,  in  Matlab 
it  is  possible  to  send  the  neural  net  inputs  in  a  parallel  or  batch  mode  (see  Matlab  Neural 
Net  manual  pages  2-13  through  2-15),  i.e.  several  sets  of  inputs  can  be  sent  through  the 
net  at  once.  Second,  Matlab  includes  both  weights  (multiplicative  constants  for  input) 
and  biases  (additive  constants)  for  each  node,  while  HNC  only  had  weights.  Finally, 
Matlab  includes  a  learning  rate  option  in  its  back  propagation  network.  As  described  in 
Chapter  5  (see  text  and  demos)  of  the  Matlab  Neural  Net  manual,  the  learning  rate 
determines  the  size  of  the  changes  to  the  weights  that  are  made  at  each  iteration. 

For  training,  Rauss  would  send  a  different  set  of  spectra  at  each  iteration,  and  simply 
iterate  until  the  solution  converged  to  a  certain  error  level.  For  HyGAS  training,  we 
chose  to  use  a  more  sophisticated  implementation  using  Matlab’s  adaptive  learning  rate 
(, trainbpx )  training.  To  begin,  the  weights  are  set  to  small  random  numbers  (as  was  done 
in  the  Matlab  example  App  cri,  below,  and  as  was  done  by  Rauss)  and  the  biases  are  set 
to  zero.  The  initial  learning  rate  was  set  to  the  Matlab  value  of  o .  oi  with  an  error  goal  of 
0 .  ooi.  For  the  four  spectra  example,  the  two  sets  of  inputs  were  trained  in  parallel:  four 
spectra  with  60  percent  noise,  and  four  spectra  with  0  percent  noise  (as  recommended  by 
Matlab  Neural  Net  manual,  example  App  cri).  The  network  iterated  on  these  inputs,  and 
displayed  results  every  20  iterations  (disp_freq=20).  A  maximum  of  5,000  iterations 
also  was  set  (max_epoch=5000). 

Using  the  adaptive  training  rate,  at  each  iteration,  the  training  rate  is  altered  based  on  the 
results  of  the  previous  iteration.  As  described  in  the  trainbpx  reference  ( Matlab  Neural 
Net  manual,  page  13-135)  and  "Improving  Backpropagation"  (p.  5-25  through  5-32),  the 
learning  rate  increases  by  a  default  of  l .  05  times  if  the  error  ratio  between  iterations  is 
less  than  the  default  value  of  l .  04.  If  it  is  greater  than  the  error  ratio,  the  learning  rate 
decreases  by  a  default  value  of  o .  7.  An  additional  parameter,  the  momentum  which 
governs  the  sensitivity  to  small  variations  in  the  error  surface,  also  was  left  at  the  default 
setting  of  o .  9. 

To  summarize,  the  example  four  spectra  net  was  trained  using  adaptive  learning  rate 
backpropagating  neural  net  (trainbpx).  The  default  values  were  used  and  the  error  goal 
was  set  to  SSE=.001.  The  Matlab  script  given  in  Script  1  shows  all  the  Matlab  calls. 


Appendix  D-3 


HyGAS  Final  Report 


Neural  Net  Notes 


In  addition  to  the  example  given  in  Script  I,  an  example  of  character  recognition  which  is 
VERY  similar  in  format  is  given  in  the  Matlab  Neural  Net  manual:  App  cri  (pp.  1 1-40 
through  1 1-46).  This  can  be  run  by  typing  appcri  at  the  Matlab  prompt  (>>). 

The  Matlab  script  m-f  ile  (Script  1)  is  similar  to  the  character  recognition  example  (App 
cri)  in  the  Matlab  Neural  Net  manual  (pp.  1 1-40  through  1 1-46).  Refer  to  that  example 
and  the  neural  net  toolbox  manual  for  explanations.  All  calls  are  documented  in  the 
neural  net  toolbox  manual  except  the  following:  rand,  zeros,  max,  bar,  all  of  which  are 
documented  in  the  Matlab  Reference  Guide.  Note  for  an  introduction  to  basic  Matlab, 
see  Matlab  User's  Guide. 


Script  1.  Matlab  script  for  training  four  spectra.  This  script  (m-file)  is  similar  to  the  Matlab  example 

App_CRl . 

%  Training  a  net  with  60%  random  noise  on  4  spectra 
%  Acetone,  Sulfur  Hexafloride,  Methanol,  Ethyl  Acetate 
% 

%  Step  1.  read  the  4  449  point  spectra  (second  column)  into  cols 
%  of  the  matrix  "input" 

load  ace.asc; 
load  shf.asc; 
load  met.asc; 
load  eat.asc; 
input ( : , 1 ) =ace (  :  ,  2 )  ; 
input ( : , 2 ) =shf ( : , 2 ) ; 
input ( : , 3 ) =met ( :  ,  2 )  ; 
input ( : , 4 ) =eat ( : , 2 ) ; 

%  make  array  "P"  with  60%  noise  added 
%  and  normalize  each  col  to  range  from  zero  to  one 

P=input  +  0 . 6*rand{449, 4) ; 
for  I  =  1:4, 

P(:,I)=(P(:,I) /max ( P ( : , I ) )  *  2)  -1; 

end 

%  make  array  "P2"  with  0%  noise  added 
%  and  normalize -each  col  to  range  from  zero  to  one 

P2 = input ; 
for  I  =  1:4, 

P2 ( : , I ) = ( P2 ( : , I ) / max ( P2 (:,!))  *  2)  -1; 

end  I 

%  change  array  "P”  to  include  both  60%  noise  and  0%  noise 
P=[P  P2] ; 

%  set  up  correct  answer  array  "T" 

T=  [  10001000 

01000100 
00100010 
00010001]; 

%  Step  2.  set  up  backpropagation  network  with  449  inputs, 

%  20  node  hidden  layer,  4  outputs,  and  log  sigmodial  transfer  functions 

R=449 ;  Sl=20 ;  S2=4; 

[Wl,bl,W2,b2]=initff (P, Si, ' logsig’ , S2 , ’ logsig’ ) ;  I 

%  set  offsets  to  zero,  weights  to  small  random  numbers, 

%  and  set  up  training  parameters 

Wl=rand ( SI , R) * . 001;  W2=rand(S2 , SI) * . 001; 
bl=zeros  (Si, 1) ;  b2=zeros (S2 , 1} ; 
disp__f  req=20 ;  max_epoch=5000 ;  err_goal=0 . 001 ; 
tp  =[disp_freq  max_epoch  err_goal] ; 

%  Step  3 .  train  net 
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[Wl,bl,W2 ,b2 , epochs ,TR] =trainbpx (Wl , bl , ' logsig' , W2,b2 , ' logsig' ,  P,  T#  tp) ; 

%  Step  4.  Apply  trained  net.  Set  up  array  "P"  with  150%  noise 
%  and  normalize  between  zero  and  one 

P=input  +  1 . 5*rand(449 , 4 } ; 
for  I  =  1:4, 

P{ :  ,  I)  =  (P( : , I) /max(P( : , I) )  *  2)  -1; 

end 

%  run  search  with  150%  array  “P"  on  trained  net  (Wl, W2 ,bl,b2) 
simuff (P,Wl,bl, 'logsig' ,W2,b2, 'logsig') ? 

%  plot  bar  graph  {see  Matlab  manual) 

I  bar  (ans)  ; _  _  _ 


In  addition.  Figure  1  shows  the  progress  of  training.  This  is  default  Matlab  graphics  that 
plots  as  you  train.  The  top  plot  charts  the  progress  of  the  SSE  at  each  iteration  until  it 
reaches  the  error  goal  (dotted  line).  The  bottom  plot  shows  the  progress  of  the  training 
rate.  It  increases  dramatically  as  the  solution  progresses.  This  adaptive  learning  rate 
allows  for  much  more  rapid  converge  of  training.  The  net  typically  converges  after  about 
160  iterations. 


T raining  for  1 66  Epochs 


Figure  1.  Neural  Net  training  progress.  This  is  a  default  Matlab  plot  showing  the  progress  of  training 
until  the  error  goal  (dotted  line)  is  reached.  The  top  graph  plots  the  sum  of  the  square  of  the  error  (SSE). 
The  bottom  graph  plots  the  learning  rate. 

The  four  spectra  selected  for  this  example  are  acetone,  sulfur  hexafluoride,  methanol  and 
ethyl  acetate.  Examples  of  the  spectra  are  shown  in  Figure  2  with  60  percent  noise,  and 
in  Figure  3,  with  150  percent  noise.  Figure  4  shows  an  example  of  passing  methanol  with 
150  percent  noise  to  the  four-spectra  net  trained  with  60  percent  noise.  The  net  identifies 
methanol  with  a  probability  of  90  percent. 
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acetone 


Sprauss  Data  +  60%  Noise  (normalized) 


Methanol 


Ethyl  Acetate 


Figure  2.  Spectra  selected  for  the  four  spectra  net  with  60  percent  random  noise  added.  The  four 
spectra  selected  for  this  test  are  acetone,  sulfur  hexafluoride,  methanol  and  ethyl  acetate. 


aoetone  Sprauss  Data  +  1 50%  Noise  (normalized) 


800  900  1000  1100  1200  1300  1400 

Sulfur  Hexafloride 


Ethyl  Acetate 


Figure  3.  Spectra  selected  for  the  four  spectra  net  with  150%  random  noise  added.  The  four  spectra 
selected  for  this  test  are  acetone,  sulfur  hexafluoride,  methanol  and  ethyl  acetate. 
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tab  spectrum 


1 50%  noise 


1400 


network  trained  at  60%  noise 


Figure  4.  Sample  results  of  passing  a  methanol  spectrum  with  150  percent  noise  through  the  four- 
spectra  net  trained  with  60  percent  noise. 


3.2.  Example  2:  Training  Sprouse  data  (65  spectra )  with  100  percent  noise. 

Training  for  65  spectra,  is  structurally  very  similar  to  training  with  four  spectra,  but  the 
implementation  is  a  bit  more  complicated.  The  only  structural  change  is  that  there  are 
now  200  nodes  in  the  hidden  layer  instead  of  20.  The  number  200  was  chosen  by  Rauss 
as  approximately  half  way  between  the  number  of  inputs  (449)  and  the  number  of  outputs 
(65). 

The  strategy  we  employed  to  train  the  Sprouse  data  was  based  on  Matlab  examples  (e.g. 
appcrl )  and  trail  and  error.  Unfortunately,  a  careful  study  to  determine  the  optimal 
method  of  training  was  beyond  the  scope  of  the  HyGAS  project.  The  Matlab  script  for 
training  65  spectra  is  given  in  its  entirety  in  Script  2. 

After  reading  the  data  into  an  array  and  normalizing  them  from  zero  to  one,  the  first  step 
in  training  was  to  train  the  neural  net  to  correctly  identify  the  65  spectra  without  any 
noise.  However,  before  training  all  65,  some  experiments  were  carried  out  with  10  and 
15  spectra.  It  was  discovered  that  with  15  spectra  the  default  learning  rate  (.01)  was  too 
high.  Many  iterations  where  required  before  the  error  when  down.  When  the  learning 
rate  was  changed  to  .005  the  net  trained  in  360  iterations.  Using  a  learning  rate  of  .001 
was  too  small  and  the  net  did  not  converge. 

Step  1 .  Training  with  all  65  spectra  began  with  the  default  Matlab  parameters  where 
used.  However,  as  with  the  15  spectra,  the  initial  learning  rate  was  too  high 
(.01);  as  such  it  was  changed  to  .006.  In  addition,  with  65  spectra  the  net  got 
stuck  in  a  local  minimum  in  the  error  surface  at  an  error  that  was  too  high  (SSE 
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on  the  order  of  1  or  10).  Changes  in  the  momentum  (0.8),  and  to  the  learning 
rate  increase  rate  (1.02)  and  decrease  rate  (0.5)  were  made  based  on  the  nets 
performance  and  trial  and  error.  Using  these  values  the  net  was  successfully 
trained  with  0  percent  noise  to  recognize  the  65  different  spectra. 

Step  2.  Starting  with  the  weights  determined  for  0  percent  noise,  training  the  net  with 
60  percent  noise  was  then  attempted.  Here  500  sets  of  65  spectra  with  60 
percent  random  noise  (different  random  noise  in  each  of  500  presentations)  were 
presented  to  the  net.  Training  parameters  were  set  tostop  at  SSE=0.001,  a 
learning  rate  of  1.0,  and  learning  rate  increase  1.5  and  decrease  of  0.5.  These 
later  numbers  were  chosen  in  attempt  to  make  the  net  convergence  faster.  By 
presenting  the  net  with  different  sets  of  random  noise  it  should  "learn"  which 
features  are  real.  After  completion  of  this  training,  the  net  worked  pretty  well 
when  data  with  60  percent  or  less  noise  was  sent  to  it.  However,  it  did  not 
perform  well  with  large  noise  amounts,  and  did  not  always  get  the  right  answer 
even  with  60  percent  noise. 

After  re-evaluating  the  results  and  the  progress  of  the  training  at  each  of  the  500 
iterations,  it  was  clear  that  the  problem  was  that  the  net  was  presented  with  a 
single  set  of  65  spectra  with  one  set  of  random  values.  The  net  iterated  until  a 
small  error  was  reached  and,  in  effect,  memorized  the  given  noise  pattern.  At 
the  presentation  of  the  next  set  of  spectra  with  a  different  random  noise,  it 
adjusted  the  weights  again  to  match  this  new  set  of  noise  "perfectly."  Over  the 
500  iterations  the  length  of  time  to  converge  decreased;  however,  the  net  was 
still  learning  a  single  set  of  random  noise.  It  was  clear  that  spectra  with 
different  sets  of  random  noise  must  be  presented  to  the  net  all  at  one.  This  is 
what  is  done  in  the  Matlab  example  appcrl. 

Step  3.  Starting  with  the  weights  after  0  percent  noise  and  500  iterations  of  60  percent 
noise,  the  net  was  then  trained  with  100  percent  noise  inputs.  This  time,  10  sets 
of  65  spectra,  each  with  different  random  noise,  were  presented  to  the  net  at  a 
single  time.  This  process  was  repeated  500  times. 

At  the  conclusion  of  this  training  the  net  performed  very  well.  Statistical  comparisons  to 
Rauss  results  for  a  net  trained  at  100  percent  noise  were  carried  out.  For  example,  the 
trained  net  was  presented  1000  sets  of  65  spectra  with  100  percent  noise  and  99.89 
percent  of  them  were  correctly  identified  (a  correct  id  is  defined  as  a  probability  of  greater 
than  0.6  by  Rauss).  When  the  same  net,  trained  to  100  percent,  was  given  1000  sets  of  65 
spectra  with  200  percent  noise  73  percent  of  the  spectra  were  correctly  identified.  Rauss 
results  were  correct  only  63  percent  of  the  time.  As  such,  the  SETS  trained  network  was 
out-performing  Rauss's  and  was  therefore  deemed  to  be  successfully  trained  to  100 
percent  noise. 

Total  training  time  for  Steps  1-3  was  greater  than  one  week  of  CPU  time  on  a  Sparc  10- 
40  machine. 

Script  2.  Matlab  script  for  training  Spouse  Library  (65  spectra). 

I  percent  load  Sprouse  Data 
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format  long  e 

load  /scooter/home/ jess /nets /work/Sprouse_449/lab_data/spec_01.prn 

(load  statements  for  other  spectra  removed  for  space) 

load  / scooter /home /jess /nets /work/ Sprouse_4 49/ lab_da ta / spec_7 2 r . prn 

waves =lab_data/spec__01  ( : ,  1)  ; 
spr=  [  lab_data/spec__01  { : ,  2 )  ... 

1 ab_da  t a / spec_0  2 ( : ,  2)  ... 

(statements  for  other  spectra  removed  for  space) 


lab_data/spec_J71  ( :  ,  2 )  ... 

lab__data/spec_„72r  ( : ,  2)  ]; 

clear  lab_data/spec_01 

(statements  for  other  spectra  removed  for  space) 
clear  lab_data/spec_72r 
%  train  with  0%  noise 
for  I  =  1:65, 

sprZERO ( : , I) = (spr ( : , I) /max(spr (:,!))); 

end 

for  I  =  1:65, 

P( : , I) = (sprZERO ( : , I) /max ( sprZERO { : , I) >  *  2)  -1; 
end 

T=zeros (65,65); 
for  1=  1:65, 

T (I , I ) =1 ; 
end 

R=449 ;  $1=200;  S2=65; 

Wl=rands (SI , R) * .001; 

W2=rands (S2 , SI) * . 001 ; 
bl=zeros (SI , 1) ; 
b2=zeros (S2 , 1) ; 

disp_f req=20 ;  max_epoch=10000;  err_goal=0 . 001 ; 

lr= . 006 ; 

lr_inc=l . 02 ; 

lr_dec=0 . 5; 

momentum=0 . 8 ; 

err_ratio=l . 01 ; 

tp  =[disp_freq  max„epoch  err_goal  lr  lr_inc  lr_dec  momentum  err_ratioJ; 
[Wl,bl,W2,b2, epochs, TR] =  trainbpx (Wl, bl , ' logsig ' , W2 , b2 , * logsig ' ,  P,  T,  tp) ; 
save  Hygas65_449 

%  train  500  times  with  60%  noise 
tp= [ 2  500  .001  1.  1.5  .5]  ; 

for  J=  1:500, 

test=  sprZERO  + . 60*rand(449 , 65) ; 
for  I  =  1:65, 

test ( : , I) = (test ( : , I) /max (test ( : , I) )  *  2)  -1; 

end 

P=test; 

[Wl , bl , W2 , b2 , epochs, TR] =  trainbpx (Wl , bl , ‘ logsig 1 ,W2 , b2 logsig ' ,  P,  T,  tp) ; 
end 

save  Hygas65_449_60 

%  train  500  times  with  10  sets  of  100%  noise 
for  K=1 : 500 ; 
test=  sprZERO; 
for  I  =  1:65, 

test ( : , X) = (test ( : , I) /max (test ( : , I) )  *  2)  -1; 
end 

T2=T ;  P2=test; 
for  J  =  1:10, 

test=  sprZERO  +rand (449,65) ; 
for  I  =  1:65, 

test  ( :  ,  I)  =  (test  (  :  ,  I) /max(test  (  :  ,  I)  )  *  2)  -1; _ _ 
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end 

T2= [T2 ,  T]; 

P2=[P2, test] ; 
end 

[Wl,bl,W2 ,b2 , epochs, TR] =  trainbpx (W1 , bl , ' logsig’ ,W2, b2, ' logsig' ,  P2 ,  T2 ,  tp)  ; 
save  Hygas„4 4 9_1 0 0_5 
end 

simuf f ( test , W1 , bl ,  ' logsig 1 ,  W2 ,  b2 , ' logsig ' ) ; 

bar  (ans)  ; _ _ _ _ _ _ _ 


4.  Interface  with  HyGAS 

The  trained  neural  net  can  be  used  to  identify  spectra  within  HyGAS,  both  in  interactive 
HyGAS  and  as  a  stand  alone  HyGAS  function.  Spectra  can  be  input  from  HIPS  data  sets. 
If  you  wish  to  run  a  new  spectrum  through  the  neural  net,  it  is  recommended  that  the 
HIPS  function  ASCII  to  hips  be  used  to  transform  ASCII  data  into  HIPS  format.  Using 
the  stand  alone  Neural  Net  function  (on  the  HyGAS  menu),  the  spectrum,  or  any  subset  of 
any  HIPS  file,  can  be  sent  to  Matlab  for  identification  against  the  Sprouse  library. 

Spectra  also  can  be  sent  from  interactive  HyGAS. 

The  Matlab  script  called  when  the  HyGAS  neural  net  routine  is  used  is  listed  in  Script  3. 
With  the  exception  of  the  one  line  which  calls  the  trained  neural  net  ( simujf  ),  all  other 
calls  are  to  routines  described  in  the  Matlab  Reference  Guide.  For  an  introduction  to 
basic  Matlab  see  the  Matlab  User's  Guide.  On-line  help  is  available  by  typing  help 
(recommend  setting  "more  on"  prior  to  help). 


Script  3:  Matlab  script  for  executing  HyGAS  Neural  Net  Identification  Routine. 


echo  off 

%  make  sure  to  NOT  echo  this  script  to  screen 

%  Load  trained  neural  net 
clear 

%  load  trained_NN  (saved  weights  and  biases,  spectra 
load  • /usr/HIPS/DATA/hygas/neural_net/trained„NN 1 

and  names  of  spectra) 

%  Load  unknown  spectrum 
clc,  disp('  '),  disp(’  ') 
f name  =  [ ' unknown ' ] ; 
eval(['load  ' , fname] ) 

spec=eval (fname) ; , spec=spec ( : , 2 ) ; , spec=spec-min(spec) 

spec=spec/max (spec) *2-1; 

%  Run  unknown  through  neural  net 

ident=simuf f (spec, Wl, bl , 'logsig' ,W2,b2, 'logsig' ) ;  ident=ident*100 ? 

[a ( : , 2 ) , a ( : , 1) ] =sort ( ident) ; ,  a=flipud(a) ; 

%  renormalize  spectrum  between  zero  and  one 
spec=spec+l; ,  spec=spec/2; 

%  Clear  or  open  new  figure  and  set  size 

figure (' Position' ,  [20  100  600  700],  'NumberTitle ' ,  ' 

Identifications ' ) ; 

off',  'Name',  'Neural  Net 

repeat= ' Y ' ; 

while  repeat (1) == *Y' ; 

%  number  of  id's  with  prob.  ?.05% 
num=size ( f ind( ident> . 05)  ,  1)  ; 

clc, disp ( '  ' ) ,  disp ( '  ' ) 

disp ( [ ' The  Neural  Net  has  identified  ' , int2str (num) , ' 

disp ( ’  '),  disp ( '  ') 

disp ( ' The  best  identification  will  be  displayed...') 

possible  spectral  matches' ] ) ; 

%  Plot  unknown  spectrum 
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subplot (311) ,plot (waves, spec) ,  title ( 'Unknown  Spectrum* ,  1  color ’  ,  'y ' ) 
xlabel (' Wavenumber ') ,  ylabel (' Arbitrary ') ,  zoom  on 

n=l; 

while  n  <=  num; 

%  Plot  id'ed  spectrum 

subplot (312) , plot (waves,  sprZERO { : , a (n, 1) ) ,  * c ' ) , axis ( [ 800  1300  0  1]) 
xlabel ( ' Wavenumber ' ) ,  ylabel ( ' Arbitrary ' ) 

title { [ 'Spectrum  #* ,  int2str (a (n, 1) ) , ' :  ’ , names (a (n, 1) , : ) ,  '  Probability: 

num2str (a (n, 2) ) , ' %  1 ] , 1  color ' ,  1 c ' ) 

%  Plot  probability  bar  graph 
subplot (313) ,bar (ident , 1 r ' ) ,hold  on 

[ x, y ] =bar (a (n, 1) , ident (a (n, 1) ) ) ; , patch (x,y, ‘c ' ) ; bar (ident,  *r 1 ) 
xlabel ( 1  Spectrum  Number  * ) ,  ylabel ( *  Probability  (percent) ’ ) 
title ( 'Neural  Net  Identifications  from  Sprouse  Library *,' color r ’) , 
hold  off 

%  Ask  if  want  to  see  next  spectrum 

in=upper ( input ( '  Display  next  best  identification?  { [Y]  or  N)  *,  *s')); 

if  size(in)==0,  in='Y';,  end 
if  in{l) == 'N* ,  break,  end 
clc; 

%  Check  to  see  if  there  are  any  more  possible  spectra 
if  n==num 
disp ( '  ' ) ; 

disp('  No  more  significant  identifications'); 

end 
n=n+l; 
end 

%  Ask  if  want  to  see  all  possible  id's 
disp ( 1  ' ) ,  disp  ( '  ' ) 

in=upper ( input ( '  Would  you  like  to  see  all  the  possible  identifications  at  once  ? 
( [Y]  or  N)  • ,  's'  ) ) ; 

if  size{in)==0(  in='Y';,  end 
if  in ( 1 ) == ' Y ' 

%  Plot  id'ed  spectrum 

subplot (312) ,  plot (waves, [sprZERO ( : , a (1 :num, 1) ) ] ) , axis ( [800  1300  0  1]) 
xlabel ( ’Wavenumber* ) ,  ylabel ( ’Arbitrary' ) 

title ( [ ‘All  Significant  Spectral  Identifications '],' color ’ ,  'w' ) 

%  Plot  probability  bar  graph 
clrs=  '  ymcrgbymcrgb,'  ; , 
subplot (313 ) ,  bar (ident ), hold  on 
for  1=1: num , 

[x,y] =bar (a (1, 1) , ident (a(I, 1) ) ) ; 
patch (x,y, clrs (I) ) ; 
bar (a (1,1) , ident (a (I, 1) ) , clrs (I) ) ; 
end 

hold  off 

xlabel ( ' Spectrum  Number' ) ,  ylabel ( ’ Probability  (percent) ' ) 
title ( 'Neural  Net  Identifications  from  Sprouse  Library'),  hold  off 
end 

repeat =upper ( input ( '  Repeat  with  same  unknown??  (Y  or  [N] )  ’,  's')); 

if  size (repeat) ==0,  repeat='N';,  end 

end; 

disp { '  ' ) ,  disp ( ’  ' ) 

disp ( '  Hit  any  key  to  end  ■ ) ; 

pause; 

close 

exit  _ _ _ _ _ 
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Gas  Coefficient  datafiles  ( GFI) 

The  gas  coefficient  datafiles  (GFI)  supplied  by  Kroutil  must  be  converted  into  a  HIPS  gas 
coefficient  library  for  use  with  HyGAS.  Kroutil  gas  coefficient  datafiles  consist  of  two 
files  per  gas: 

•  discrim.inc 

•  filter.inc 

These  files  must  be  edited,  renamed,  and  put  into  a  HIPS  gas  coefficient  library  directory 
structure  to  be  used  by  HyGAS.  The  HIPS  gas  coefficient  library  directory  structure  is  a 
set  of  two  ASCII  files,  discrim.dat  and  matrix .  dat,  for  each  gas  (Table  1).  Example 
listings  of  these  files  for  the  gas  acetone  are  given  in  Tables  2  and  3.  Tables  4  and  5  list 
the  contents  of  the  corresponding  original  Kroutil  datafiles  (discrim2 .  inc  and 
f  Uter2 .  inc),  from  which  the  HIPS  format  files  were  derived.  The  HIPS  format  files 
are  essentially  identical  to  the  Kroutil  format  files  except  for  the  various  text  and 
delimiters  included  in  the  Kroutil  format  files.  These  files  can  be  edited  with  any  basic 
text  editor  (the  unix  editor  vi,  for  example)  to  eliminate  the  extraneous  characters  and 
convert  them  into  the  appropriate  format  for  a  HIPS  gas  coefficient  library. 

Table  1.  HIPS  gas  coefficient  library  directory  structure  for  HyGAS.  Two  gases,  acetone  and  MEK, 
are  included  in  this  library.  For  each  gas,  two  files  are  required,  discrim .  dat  and  matrix .  dat.  These  files 
are  derived  from  the  original  Kroutil  (GFI)  gas  coefficient  files. 


$ DOD_DIR/ LIBRARY /gas_Coeff : 
acetone/  mek/ 

$DOD_DIR/LIBRARY/gas_coe£f /acetone: 
discrim.dat  matrix.dat 

$DOD_DIR/LlBRARY/gas_coef f /mek: 
discrim.dat  matrix.dat _ 
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Table  4.  Example  listing  of  coefficient  file  discrim2  .±nc,  as  supplied  by  Kroutil  (GFI),  for  acetone.  This  is  the  original  Kroutil  file  discrim2 .  inc. 
equivalent  HyGAS  directory  path  for  this  file  is  /d2/hygas/Kroutil/new/mtrx/discrim2 .  inc . 
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Table  5  Example  listing  of  coefficient  file  f  iiter2 .  inc,  as  supplied  by  Kroutil  (GFI),  for  acetone.  This  is  the  original  Kroutil  file  f  ilter2 .  inc. 
equivalent  HyGAS  directory  path  for  this  file  is  /d2/hygas/Kroutil/new/mtrx/f  iiter2 .  inc . 
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},  -0. 1015457 E- 02,  0 . 17873 16E-01 ,  -0 . 5627253E-01 ,  0 . 2446596E-01 , 

{  0 . 4138545E-01 ,  -0 . 4920137E-01 ,  0 . 8277126E-01 ,  -0 . 5583 157E-01 , 

-88,  -20,  -14,  -13,  -8 ,  -7,  0,  0 . 5483646E-01,  0 . 7511845E-01,  0 . 7321317E-01 ,  0 . 5719284E-01 , 
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Introduction: 

Model-based  automatic  target  recognition  (ATR)  is  the  focus  of  much  attention.  One  approach  is 
the  use  of  artificial  neural  networks  (NNs)  that  have  been  trained  with  model  data.  This  paper  will  present 
the  results  of  a  research  effort  entitled  Neural  Networks  Applied  to  Infrared  Signal  Processing.  This 
project  evolved  from  an  effort  to  automatically  identify  infrared  chemical  spectra  collected  in  the  field 
with  a  Fourier  Transform  Infrared  Spectrometer  (FTIR).  A  number  of  common  spectroscopic  analysis 
techniques  were  investigated  and  found  unable  to  tolerate  the  level  of  noise  common  to  this  sensor  sys¬ 
tem.  i  Due  to  the  poor  performance  of  the  common  techniques  a  new  approach  was  required.  A  literature 
search  was  undertaken,  as  well  as  discussions  with  other  Army  organizations,  to  identify  methods  to  pull 
signals  out  of  a  high  noise  environment.  A  number  of  promising  methods  were  identified  and  subse¬ 
quently  ruled  out  due  to  constraints  on  the  program  hardware.  Several  sources  mentioned  NN  s  tolerance 
to  noise,  suggesting  these  may  be  a  promising  approach. 

Overview: 

The  neural  network  chosen  to  investigate  was  the  simple  back  propagation  of  error,  or  backprop 
(BP),  network.  The  BP  network  is  a  supervised  network.  The  architecture  used  in  this  work  consists  of 


fully  interconnected  in  each  layer. 


Figure  2:  Processing  Element 
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three  layers;  input,  hidden,  and  output.  The  input  layer  consists  of  571  nodes.  These  nodes  are  fully  con¬ 
nected  to  200  processing  elements  (PEs)  in  the  hidden  layer.  The  hidden-layer  PEs  are  fully  connected  to 
66  output-layer  PEs  (see  fig.  1).  The  NN  is  given  an  input  and  told  what  the  required  response  or  truth  is. 
The  network  compares  its  response  to  the  desired  response,  computes  an  error,  and  propagates  the  error 
correction  back  through  the  layers.  The  PEs  take  the  sum  of  the  weighted  inputs  and  apply  a  nonlinear 
thresholding  function.  The  function  used  in  this  work  is  the  logistic  sigmoid  function,  as  shown  in  figure 
2.  The  output  of  the  PE  depends  upon  the  activation  of  the  sigmoid  function. 

The  current  environment  used  for  this  work  consists  of  a  VME  BALBOA  i860  neural  accelerator 
board  and  NeuroSoft  software  from  HNC,  Inc.,  hosted  on  a  SUN  630MP.  The  program  is  written  in  C 
with  function  calls  to  NeuroSoft. 

The  goal  of  this  effort  is  to  identify  an  unknown  infrared  spectrum  as  one  of  a  library  of  66  spectra 
or  as  an  unknown/unidentifiable.  The  structure  of  this  NN  with  66  output  nodes,  one  for  each  class,  has 
been  shown  to  give  Bayesian  a  posteriori  probabilities,  once  trained,  if  the  training  requires  the  nodes  to 

respond  in  a  binary  fashion:  true  =  1  and  false  =  0.^  The  neural  network  is  trained  on  a  library  of  66 
spectra  that  were  collected  in  a  laboratory  with  a  laboratory-grade  FTIR  configured  to  collect  at  the  same 
resolution  as  the  fieldable  sensor.  The  library  data  has  no  noticeable  noise  in  the  signal. 

The  problem  with  training  an  NN  in  this  fashion  is  its  ability  to  memorize.  In  this  environment  the 
training  of  an  NN  continues  until  some  Mean  Squared  Error  (MSE)  threshold  is  crossed.  That  is,  over  the 
last  complete  pass  through  the  training  file  did  the  MSE  drop  below  the  set  threshold;  if  so,  training  is 
complete.  Equation  (1)  shows  the  MSE  calculation  used  by  this  package. 

N  M 

=  S<',-zi>2  <*> 

k =  1  i=l 

where  k  indexes  the  iterations,  one  iteration  being  the  presentation  of  one  spectrum;  N  is  the  number 
of  consecutive  iterations  of  the  network  since  the  last  error  calculation,  here  set  to  66  (check  MSE  after 
each  training  pass  through  the  library);  i  indexes  the  PE’s;  M  is  the  total  number  of  PE  s;  t  is  the  truth 
value  for  the  ith  PE,  and  z  is  the  output  of  the  ith  PE. 

The  difficulty  lies  in  answering  the  question:  What  is  a  suitable  MSE  threshold?  The  MSE  thresh¬ 
old  could  be  set  to  zero,  which  would  require  the  network  to  train  until  it  exactly  reproduced  the  required 
outputs.  However,  in  this  case  that  would  mean  the  NN  memorized  the  original  spectral  inputs,  and  exact 
reproductions  of  these  spectra  would  be  required  to  get  a  proper  response  from  the  network.  The  problem 
is  the  actual  sensor  will  have  quite  a  bit  of  noise  in  it,  which  would  confuse  a  network  trained  to  this  level. 

Alternatively  one  could  train  only  to  a  relatively  high  MSE:  however,  this  network  would  have  poor 
recognition  qualities  in  general  or  would  identify  some  of  the  spectra  well  and  others  not  at  all.  The 
objective  is  to  find  an  optimal  MSE  which  identifies  all  the  spectra  but  is  not  trained  so  far  that  it  needs  an 
exact  reproduction  of  the  spectra  it  was  trained  on  before  it  registers  a  match. 

In  order  to  avoid  having  to  collect  the  library  with  the  fieldable  sensor  or  to  extensively  preprocess 
the  collected  data,  the  approach  was  taken  to  add  random  noise  to  every  input  point  to  ‘fuzz  up’  the 
model  data.  The  questions  are  how  much  noise  can  be  added  before  the  network  will  no  longer  learn  and 
how  well  will  an  NN  trained  in  this  manner  perform  on  sensor  data.  Figure  3  shows  an  example  of  the 
original  library  spectra  and  the  resultant  noisy  spectra.  This  approach  also  forces  the  network  to  become 
very  robust  by  compelling  it  to  focus  only  on  the  most  robust  features  of  the  data  set.  The  network  will 
not  focus  only  on  peak  location  or  relative  intensities  between  specific  peaks  because  these  values  will 
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continuously  vary  during  training. 

Approach: 

The  following  procedure  was  used  to  train  the  networks  discussed  in  this  paper.  A  network  was  ini¬ 
tialized  with  small  near-  zero  random  weights.  A  noise  level  was  chosen.  This  level  is  a  percentage  of  the 
height  of  the  preprocessed  library  spectrum.  The  preprocessing  simply  normalizes  each  spectrum 
between  0  and  1.  The  noise  is  added  to  the  spectrum  just  before  presentation  to  the  NN.  The  resulting 
noisy  spectrum  is  then  normalized  again  to  either  a  0  to  1  or  - 1  to  1  range  and  fed  to  the  NN  for  training. 
This  paper  primarily  focuses  on  the  - 1  to  1  range  as  this  range  appeared  to  train  faster  in  preliminary 
work. 

The  noisy  spectra  are  renormalized  in  order  to  make  the  baseline  of  the  spectra  vary  somewhat.  Of 
course,  the  amount  of  noise  added  will  influence  how  much  the  baseline  varies.  There  was  some  concern 
that  the  NN  may  only  recognize  spectra  with  a  baseline  on  the  zero  or  - 1  line  and  thus  not  recognize  sen¬ 
sor  data  that  can  have  a  varying  baseline  due  to  background  subtraction  methods;  figure  4  shows  an 
example  of  very  bad  baseline  from  sensor  data. 

Once  training  is  begun  the  MSE  is  monitored  every  complete  pass  through  the  training  library. 
When  the  MSE  drops  below  0. 1  the  network  state  is  saved.  Training  is  then  continued  from  the  current 
state  to  the  next  threshold.  The  thresholds  used  in  this  effort  are  0.1, 0.01, 0.001,  and  0.0001.  Four  thresh¬ 
olds  were  chosen  to  allow  an  investigation  into  what  would  be  a  good  MSE  level  for  training.  Too  much 


Figure  4:  Sensor  data  from  field  collection. 
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training,  even  if  there  is  a  slight  improvement  of  performance,  can  be  a  waste  of  resources. 

Networks  have  been  completely  trained  on  the  following  noise  levels:  20%  (2),  50%,  60%,  70%, 
80%,  90%,  100%  (7),  140%,  150%  (2),  and  200%.  The  (2)  indicates  that  two  networks  were  trained  start¬ 
ing  in  different  random  areas  of  weight  space.  This  was  done  to  ensure  that  a  trained  network  did  not  just 
“happen”  to  start  in  a  weight  space  area  near  a  solution.  The  seven  different  networks  for  100%  noise  are 
the  result  of  an  error  in  the  computer  code,  which  defaulted  to  100%  noise  for  any  value  entered  greater 
than  100%.  This  error  was  corrected  and  networks  trained  recently  to  140%,  150%,  and  200%  noise  lev¬ 
els. 

In  order  to  determine  the  noise  tolerance  of  the  trained  networks  a  statistics  collection  routine  was 
written  that  uses  the  library  spectra  trained  with,  as  a  test.  The  library  spectra  are  used  at  this  stage 
because  a  truth  file  exists  that  is  used  to  determine  a  correct  hit.  Starting  with  0%  noise  added,  the  library 
is  run  through  the  trained  network  with  learning  turned  off.  A  threshold  of  60%  is  set  for  a  confidence- 
level  check.  The  output  nodes  are  sorted  and  the  largest  response  is  checked.  If  it  does  not  cross  the 
threshold,  the  presentation  is  classified  as  NoCross.  If  it  crosses  or  equals  the  threshold  it  is  checked  for 

Network  trained  with  100%  noise 
Threshold  level  is  0.6 
Ambiguity  level  is  0.4 


Noise  level  is  0  Percent. 
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Figure  5:  Sample  output  of  statistics  collection  program. 
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correctness  against  the  truth  file.  If  it  is  incorrect  the  False  counter  is  incremented  and  the  node  number  of 
-the  false  spectra  is  logged  in  an  incorrect  hit  table  (see  fig.  5  for  an  example).  If  the  largest  responding 
node  is  the  correct  identification,  the  Correct  counter  is  incremented  and  the  node  number  of  the  spectra 
is  incremented  (fig.  5).  One  further  check  is  made  at  this  point.  An  ambiguity  level  is  set  when  the  pro¬ 
gram  starts.  For  most  of  this  work  it  is  0.4.  This  means  that  if  a  correct  hit  is  logged  a  comparison  is  made 
with  the  next  highest  responding  node.  If  they  are  less  than  0.4  apart  the  multi-cross  counter  is  incre¬ 
mented. 

Figure  5  shows  the  output  of  the  statistics  collection  program  for  a  network  trained  with  100%  ran¬ 
dom  noise  and  trained  to  0. 1  MSE.  The  noise  level  indicated  in  figure  5  is  the  noise  level  presented  to  the 
trained  network,  here  0% .  The  next  block  of  numbers  are  the  correct  hit  registers  for  each  spectra,  show¬ 
ing  here  2000  presentations  for  each  spectra.  The  left-most  column  is  a  row  index.  The  next  number  to 
the  right  is  the  number  of  hits,  and  the  following  number  to  the  right  is  the  spectrum  number.  The  next 
block  down  are  the  false  alarms  of  each  spectra.  Below  are  the  number  of  multi  crossings,  no  crossings  of 
the  threshold,  total  number  of  false  alarms,  total  number  of  correct  hits,  total  number  of  spectra  presented, 
and  the  false  alarm  and  correct  hit  rates  in  percentages. 

Once  the  program  has  presented  the  desired  number  of  spectra,  the  results  are  written  to  disk  and 
the  program  then  increments  the  noise  level  10%  and  repeats  for  the  same  number  of  total  passes.  The 
program  continues  until  it  completes  collecting  results  on  the  upper  noise  level  bound  that  was  set.  For 
this  work  the  upper  bound  is  200%  noise. 

Results: 

The  first  observation  to  be  made  upon  examination  of  the  statistics  files  collected  thus  far  shows 
that  training  to  an  MSE  of  0. 1  is  insufficient  to  ensure  complete  training.  At  this  level,  at  least  three  of  the 
spectra  are  not  correctly  classified  (see  fig.  5).  In  most  cases  the  spectra  not  identified  correctly  are  not 
identified  at  all,  as  in  the  results  shown  in  figure  5.  In  a  few  cases  one  of  the  spectra  is  consistently  classi¬ 
fied  incorrectly  as  another  spectrum. 

The  second  observation  noted  is  shown  in  figure  6.  Up  to  the  noise  trained  on  there  is  little  or  no  dif¬ 
ference  between  the  performance  of  networks  trained  to  an  MSE  of  0.01, 0.001,  or  0.0001.  This  observa- 
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tion  holds  for  networks  trained  with  up  to  100%  noise.  On  networks  trained  with  noise  levels  between 
100%  and  150%  noise,  networks  trained  to  an  MSE  of  0.0001  incorrectly  identify  one  spectrum  at  low 
noise  levels.  This  can  be  seen  in  figure  7  as  the  slight  dip  at  the  low  noise  end  of  the  150%  curve  and  also 
in  figure  8.  As  the  noise  level  during  testing  increases,  the  network  begins  to  correctly  identify  this  spec¬ 
trum  more  regularly  until  the  performance  of  the  entire  network  begins  to  degrade  at  higher  noise  levels. 
The  exact  training  noise  level  where  this  begins  has  not  been  identified  yet,  since  statistics  have  not  yet 
been  collected  for  all  the  nets  trained  on  noise  levels  between  100%  and  150%,  although  figure  8  shows  a 
similar  effect  for  the  network  trained  with  140%  noise.  This  activity  may  indicate  some  over  training  of 
the  network,  so  it  appears  that  at  least  for  high  noise  training  levels  an  MSE  of  0.001  is  more  useful. 

With  the  60%  confidence  threshold  before  classifying  a  node  active,  the  false  alarm  rate  remains 
very  low  even  as  the  classification  performance  decays  (see  figs.  8  and  9).  Figure  8  also  shows  the  incor¬ 
rect  classification  of  a  single  spectrum  by  nets  trained  with  high  noise  and  trained  to  an  MSE  of  0.0001 
when  processing  low  noise  spectra.  By  comparing  figures  8  and  9  we  see  that  by  using  nets  trained  to  an 
MSE  of  0.001,  in  general,  there  is  very  little  increase  in  false  alarms,  except  for  nets  trained  with  very  lit¬ 
tle  noise. 

Preliminary  results  of  collecting  statistics  with  no  thresholds  set  (i.e.  finding  the  node  with  the  larg¬ 
est  magnitude)  show  that  some  networks  generalize  well  out  to  nearly  three  times  the  noise  level  they 
were  trained  on  (fig.  10).  Also  note  that  the  false  alarm  rate  also  increases  much  faster  because  every  pre¬ 
sentation  requires  an  identification  (fig.  1 1).  At  this  time  no  information  is  available  to  indicate  how  close 
the  next  highest  node  is  under  these  conditions.  Nor  is  information  collected  on  the  intensity  of  the  best 
node.  It  is  likely  that  both  of  these  values  would  show  low  node  response  and  little  difference  between 
best  and  second  best  hits  at  very  high  noise  test  levels. 

Conclusions: 

The  results  to  date  show  that  by  training  with  noise  added  to  the  model  signal  a  neural  network  can 
be  generated  which  has  high  noise  tolerance.  For  this  data  set,  a  reasonable  noise  level  to  train  with 
appears  to  lie  between  100%  and  140%  and  training  to  an  MSE  of  0.001  appears  adequate. 


Percent  noise  added  for  test 


Figure  8:  Trained  to  MSE  of  0.0001 -False  alarms. 


Figure  9:  Trained  to  MSE  of  0.001 -False  alarms. 
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Preliminary  indications  of  the  NNs  performance  on  sensor  data  are  encouraging,  but  only  a  small 
number,  <20,  of  sensor  spectra  have  been  run  through  a  network  trained  with  50%  noise.  The  results  were 
mixed  as  expected  because  some  of  the  sensor  data  had  very  skewed  baselines. 

Finally,  by  requiring  a  minimum  confidence  level  for  identification  one  can  suppress  many  false 
alarms  down  to  1  to  2%  and  maintain  very  high  classification  rates. 

Future  Efforts: 

A  more  thorough  test  of  performance  with  sensor  data  will  be  conducted  shortly.  The  data  is  cur¬ 
rently  undergoing  background  subtraction  and  will  be  formatted  for  input  into  the  trained  networks.  A 
small  rewrite  of  the  statistics  software  may  be  needed  to  properly  record  the  results  to  allow  a  determina¬ 
tion  of  the  best  noise  level  to  train  on  for  sensor  data  identification. 

Also  further  statistics  will  be  collected  out  to  higher  noise  levels  to  determine  if  the  networks  all 
begin  exhibit  a  linear  roll-off  at  high  noise  levels  as  seen  for  the  50%  network  in  figure  10.  Optimal  con¬ 
fidence  thresholds  will  also  be  determined 

Shortly,  these  techniques  will  be  applied  to  Forward  Looking  Infrared  (FLIR)  imagery  in  order  to 
develop  a  FLIR  ATR  system  that  exhibits  the  high  classification  rates  displayed  under  this  effort.  Appli¬ 
cation  of  this  technique  to  FLIR  ATR  should  not  only  allow  training  with  a  model  or  simulated  data,  but 
should  also  improve  the  performance  of  the  ATR  under  dynamic  battlefield  conditions. 


Percent  noise  tested  with 

Figure  10:  Correct  hit  rate-No  Threshold. 


Figure  1 1 :  False  alarm  rate-No  Threshold. 
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